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ABSTRACT 


Accurate measurement of particle size distribution of rocket motor exhausts is 
essential for predicting the combustion efficiency and infrared plume signature. This 
thesis presents an automated method for extracting particle size distribution from 
scanning electron microscope (SEM) images. The SEM images were taken off a 
filter paper placed at the end of a collection probe inserted into the exhaust plume. 
The automated SEM extraction system consists of an IBM AT-based computer 
system fitted with a 512 x 480 pixel frame grabber. Photographic images taken off 
the SEM are acquired via a vidicon camera. A C language program was written to 
control the hardware and automate the extraction process. A threshold is first 
applied to the digitized image and the resulting binary image is subjected to object 
segmentation. Each object is then sized and the distribution from one or more 
images can be plotted. The main bulk of this thesis is to document the software 
specially written to undertake this set of tasks. Results obtained were compared with 
that from a Malvern MasterSizer particle sizer and found to be favorable. Particles 


as small as 1/8 um have been successfully sized. 
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I. INTRODUCTION 

The rationale for analyzing scanning electron microscope (SEM) images of 
combustion exhaust stems from the on-going research to investigate particle behavior 
in exhaust nozzles and plumes of solid propellant rocket motors. Netzer and Powers 
describe various methods being tried in the quest to obtain accurate particle size 
distribution [Ref. 1]. Accuracy is essential because the various computer 
codes available for predicting performance and infrared signatures are highly 
sensitive to the particle size distribution. 

Traditionally, SEM images are analyzed visually by a trained observer, with the 
aid of linear scales (stage or eyepiece micrometer), graticules (British Standard BS 
3406), or other visual cues. Full use is made of the human observer’s ability and 
talent to differentiate shapes and features. Humphries [Ref. 2] has a concise 
description of the various methods. Disadvantages include low productivity, fatigue, 
observer subjectivity, and a high probability of error. The generally small sample 
population which a human observer is able to analyze also results in sampling errors 
due to poor statistical averaging. 

The standard procedure in particle size analysis (PSA) is to allot each particle 
to its appropriate size class and, from the resultant data, calculate (or determine 
graphically) parameters that adequately describe the size distribution. These 
parameters include particle projection lengths, perimeter, form factor, area and 


volume. Automatic PSA offers a significant advantage over manual microscopy in 


that it allows large numbers of particles to be measured consistently and accurately. 
However, automatic PSA is not without limitations. Resolution and contrast 
performance of the imager is generally poorer than the human eye. Also, the 
instruments for such systems tend to be rather elaborate and expensive. 

The current effort is a spin-off from work done with holograms, first by 
Redman [Ref.3] and Orguc [Ref.4] using the FORTRAN language, and 
subsequently by Kaeser [Ref. 5] and Hockgraver [Ref. 6] in the C language. The 
C algorithms used in this thesis have been optimized for speed, yet have reduced 
false feature identification. In the earlier work with holograms, significant problems 
were encountered with speckle in the reconstructed hologram images. This limited 
the particle size determination to particles larger than 10 microns [Ref. 7]. 
The current work attempts to analyze particles down to the resolution limits of the 
imaging system. For the SEM images analyzed, particle size distributions down to 
one-eighth of a micron have been achieved. 

This thesis is divided into six chapters. Chapter II describes the SEM images, 
what they represent and how they are obtained. Chapter III provides an overview 
of all the hardware and software required to support this work. Chapter IV goes into 
the detailed description of the program modules making up the Scanning Electron 
Microscope Extraction (SEMEX) program which forms the bulk of the current effort. 
Chapter V describes the procedures required to acquire a SEM image from a 
photograph, to process it to a form suitable for extraction of particles, and, finally, 


to extract and compile the results for plotting or further statistical analysis. Chapter 


VI describes the experimental results beginning with calibration, speed performance 
comparison with earlier programs, and finally summarizes the results from SEM 
images obtained from an actual motor burn. The latter are correlated with results 
obtained from the Malvern MasterSizer particle sizer [Ref. 8]. The last 
chapter lists the conclusions arrived at and recommendations for future efforts in this 


area. 


Il. SEM IMAGES OF COMBUSTION EXHAUST 
This chapter introduces the combustion mechanisms involved in metallized 
solid-fuels and solid-propellants, and the need for accurate measurement of the 
plume particle size distribution. (One such method using SEM images of the 


collected exhaust is described. 


A. COMBUSTION PHENOMENA 

The higher energetic performance of metallized solid-fuels (1.e., fuels containing 
boron, magnesium, titanium, aluminum, etc.) has sparked considerable interest in 
solid-fuel ramjets. While the exact combustion phenomena for metallized fuels has 
not been fully understood, Gany and Netzer [Ref. 9] ee one possible 
scenario, shown in Figure 1. During combustion, the hydrocarbon fuel vaporizes at 


the surface exposing the metal fuel particles. As there is little or no oxygen at the 
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Figure 1. Combustion and flow characteristics of a metallized solid-fuel 
[From Ref. 9:p. 424]. 


fuel surface, these metal particles can heat up but will not ignite. As more of the 
surrounding hydrocarbon fuel vaporizes, the metal particles tend to coalesce before 
being ejected into the main flow which is rich in oxygen. A particle may collide with 
other particles and agglomerate, or it may contact oxygen and ignite intensely. 
Agglomeration of metal particles may be one reason for poor combustion efficiency. 
Another may be due to the oxide coating on the particle, thus slowing the rate of 
chemical reaction. Agglomeration and oxidation also result in a large variation in 
particle sizes. 

Metals are also added to solid-propellants for rocket applications to provide 
combustion stability and/or increased performance. Aluminum is the most often 
used performance enhancer. The particles burn to produce aluminum oxide of 
varying sizes. When these particles pass through the exhaust nozzle, they can change 
in size distribution due to breakup and collisions. This two-phase flow results in a 
loss in delivered thrust. When these particles exit the nozzle into the plume, they 
significantly affect the plume radiation. The infrared signature of the exhaust plume 
is a function of the temperature profile of the plume and this, in turn, is dependent 
On its composition of hot gases and metal particles. From high-speed infrared 
imaging, it is noticed that the center of the plume appears hottest. Some of this is 
due to afterburning of the fuel-rich exhaust gases with the ambient oxygen. It is also 
suspected that the larger particles cannot turn as rapidly as the gas within the exhaust 
nozzle, resulting in these particles being located nearer to the plume centerline. 


They also have a greater heat capacity and, hence, cool off slowest after exiting from 


the exhaust. Accurate prediction of plume radiation requires accurate knowledge of 
the particle size distribution at the nozzle exit. To validate the codes, axial and 
radial variations in the plume particle size distribution need to be accurately 
measured. 

The need for accurate determination of particle size and their distribution is 
recognized and accurate methods are available to determine modal distributions. 
However, a complication lies in the multi-modal distributions caused by the wide 
range of particle sizes expected. Particle sizes may vary from sub-micron sizes to 
several tens of microns. For particles larger than about 2 um, various methods, such 
as forward laser-diffraction measurements [Ref. 10] have been successful. For sizes 
down to 0.5 um, the Malvern MasterSizer particle sizer has been used; utilizing Mie 
corrections to the diffraction equations [Ref. 8]. Measurements using these 
techniques have not been successful for particles in gas flows smaller than 0.5 um 
because of the requirement for short focal length lenses and the resulting restricted 
measurement volumes. Analysis of the collected exhaust particles using the SEM 
avoids these problems and shows promise as it allows for very clean images even at 
high magnification. Resolutions down to one-eighth of a micron have been possible, 


as is shown in this work. 


B. SOURCE OF SEM IMAGES 
The SEM images are extracted from combustion exhaust taken from a small 
(5 cm) experimental rocket motor developed by the Naval Postgraduate School. 


Different propellant compositions are being tested. The firings and subsequent 


extraction of the SEM images were carried out by CAPT Lyle Kellman, a thesis 
student working under Prof. David Netzer. A detailed description of his work is 
found in his thesis [Ref.11]. Each firing produces two sets of Malvern 
MasterSizer data, one generated directly during the burn and another from dissolving 
the filter paper removed from the collection probe. Part of the filter paper from the 
collection probe is also analyzed under the SEM. The resulting SEM images are the 


inputs to SEMEX and the focus of this thesis study. 


1. Experimental Setup 

The combustion exhaust from the rocket motor is sampled with the aid of 
a collection probe. The probe is placed at a desired location behind the exhaust 
nozzle and is protected from the intense heat of the exhaust by a plume deflector. 
The deflector prevents the plume from impinging on the collection probe tip until 
the rocket motor achieves steady-state, and is then lowered for one second. The 
exhaust enters the collection probe through a small entry nozzle and expands inside, 
depositing the combustion products on a filter paper. At about the same time, the 
Malvern MasterSizer is activated to scan the captured exhaust plume and to store its 


results for later analysis. A schematic of the setup is shown in Figure 2. 


2. Preparation of SEM Specimens 
After firing, the filter paper is removed from the rear of the probe. Two 
half-inch diameter samples are:cut out from the edge and center of the filter paper. 


These samples are gold-plated in preparation for scanning under the SEM. 
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Figure 2. Experimental setup to collect exhaust particles from a solid- 
propellant rocket motor [Ref. 11]. 


The rest of the filter paper is dissolved in acetone. The metal particles are 
allowed to settle to the bottom and the supernatant liquid is siphoned off. This 
process is repeated several times until the paper is fully dissolved and removed. In 
the final step, acetone is added again and the mixture of acetone and solid particles 
is agitated to get a homogeneous mixture. A drop of this mixture is extracted and 
analyzed with the Malvern MasterSizer. This produces a second set of Malvern 
results. 

One source of error comes from dirt or other debris being introduced on 


the filter paper during combustion and subsequent handling. Sampling errors can 


occur if the half-inch samples are not taken from representative portions of the filter 


paper where particles are uniformly distributed. 


3. Photographing the SEM Images 

The prepared specimens are placed inside a Hitachi Model S450 SEM and 
the SEM chamber evacuated to a high vacuum. The specimens can be moved about 
in search of representative portions of the specimens. This is carried out with the 
SEM on low magnification. The area being scanned is displayed on a high- 
persistence CRT. Once particles are found, the magnification is increased to 1,000. 
The scanning electron micrographs are then recorded on Polaroid photographic film 
for subsequent analysis by SEMEX. The resulting images have between 5 and 1,000 
particles, depending on the type of propellant used, the time of exposure, the 
distance of the collection probe from the exit nozzle, the radial position of the probe, 
and the location of the filter paper from which the sample was taken. To increase 
the number of particles on a single photograph, multiple exposures were taken from 
different fields. It was found that up to four exposures could be taken without 
significant degradation to the image. The SEM was set for maximum contrast and 
minimum brightness so that the background stays relatively dark even after four 
exposures. 

A source of error comes about from the choice of the area to be 
photographed. This is a subjective task. Many of the areas may contain only a few 
particles while a few areas may contain a large splattering of particles. Each 


photographed area (70 um by 80 um) only represents about 1/22,500th of the sample 


area and typically only 30 to 40 different areas can be photographed in a 4-hour 
session. Thus, the size distribution of the particles can be heavily biased by the SEM 


operator if the particles are not evenly distributed. 


C. CURRENT EFFORTS 

Current work carried out by the Naval Postgraduate School for the Air Force 
Phillips Laboratory include analysis of the exhaust plume by high-speed video and 
infrared cameras, programmed data acquisition of temperature and pressure using 
transducers, and analysis of the exhaust plume composition by the Malvern 
MasterSizer and SEMEX. The experimental data will be used both as inputs to and 
for validation of plume prediction computer codes. The work is primarily carried out 
by the Aeronautics Department with the Electrical and Computer Department 
assisting in the area of SEM image digitization and feature extraction. 

In this thesis, the work done encompasses the development of the SEMEX 
program and the setting up of a reliable methodology for the digitization and feature 
extraction of the SEM images. Calibration of the system has been carried out and 
comparisons have been made with the Malvern MasterSizer. The performance of 
several SEMEX modules (in particular, CLIP, TAG and SIZE) have also been 
compared with relevant portions of the HOLOGRAM C program developed by 
Kaeser [Ref. 5] and improved upon by Hockgraver [Ref. 6]. The results have been 


favorable. 
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III. OVERVIEW OF THE AUTOMATED SIZING SYSTEM 
The automated particle sizing system consists of an IBM AT microcomputer 
fitted with dedicated hardware and run by a program written in the C language. 


Figure 3 shows the equipment involved in the extraction of particle size data from 
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Figure 3. Diagram showing the equipment setup used for extracting 
particulate size data from SEM photographic images [After 
Ret. 12]: 


SEM photographic images. SEM photographic images are acquired by the vidicon 
camera and digitized by the frame grabber. Processing of the images is carried out 
by the IBM AT and the processed image is displayed on the video monitor. 


Digitized images in the frame memory can be saved to disk for later retrieval. User 
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dialog is carried out through the use of the computer monitor and the keyboard. A 
hardcopy of the digitized image can be printed on the video copy processor. The 
details of the hardware configuration and the software support are described in the 


rest of this chapter. 


A. HARDWARE CONFIGURATION 
The system configuration consists of the following hardware: 


- IBM AT with Intel Inboard 386/AT accelerator board and EGA monitor. The 
Inboard 386/AT disables the Intel 80286 microprocessor inside the IBM AT 
and executes programs with the 32-bit Intel 80386 microprocessor running at 
16 MHz. This enhancement gives a speedup of 17.6 over a standard IBM PC. 
The computer also has a 80287 math coprocessor running at 10 MHz to speed 
up floating point operations. The EGA monitor functions in text mode and 
displays the dialog boxes. The user makes his or her choices by typing in 
commands on the keyboard. 


e Imaging Technology PCVISIONplus frame grabber board. This is a video 
digitizer and frame memory capable of digitizing the standard RS-170 
composite video input from the vidicon camera. Figure 4 shows a block 
diagram of the frame grabber. Two camera inputs are available. Each camera 
input, in turn, has three input channels (red, blue, and green) to handle pseudo- 
color images. However, only the green channel is being used for monochrome 
images because it carries the synchronization signals. The built-in gain and 
offset circuits can be adjusted through software control to achieve optimum 
brightness and contrast. The look-up tables (LUTs) are special memories used 
to transform pixel values without altering the contents of frame memory. 
Although there are 32 LUTs (8 input LUTs, and 8 output LUTs for each of the 
red, blue and green channels), only one input LUT and one output (green) 
LUT are required. The digitized image, stored in a special on-board high- 
speed memory (called the frame memory), can be manipulated by specially 
written software functions and saved to a disk file for later analysis. The frame 
memory can handle up to two images, each 512 x 512 pixels wide. However, 
the RS-170 video signal can only provide 512 pixels by 480 lines. Therefore, 
the 481st to 512th rows of the frame memory are not used. The digitized 
image has to be converted back to analog RS-170 video signals before the 
image can be displayed on a composite video monitor. [Ref. 12] 
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Figure 4. Block diagram of the PCVISIONplus frame grabber board 
[From Ref. 12]. 


e Panasonic Model WV-1410 CCTV camera and light table. The Closed-Circuit 
Television (CCTV) camera uses a 525-line vidicon imaging element with 
25 mm f/1.4 lens and is able to operate down to 5 lux incandescent 
illumination. The internal electronics scan and convert the image into 
composite monochrome video which is routed to the frame grabber. The light 
table (not shown) consists of an adjustable mount for the camera and four 
75 W incandescent lamps fitted with diffusers. Setting the height at 
approximately 13 inches from the SEM photograph enables the 3.5" by 4.5" 
photographic image to be captured onto the 512 x 480 pixel video frame. The 
camera and illumination setup Is critical for obtaining a good digitized image. 
A piece of flat transparent glass, used to hold the photograph flat, completes 
this setup. Dust particles should be excluded. The placement of the lamps is 
crucial for obtaining even illumination. A detailed procedure is provided in 
Chapter V. 


* Panasonic TR-196M monochrome composite video monitor. This is used to 
display the digitized SEM image and any video operations carried out on it. 
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on it. The horizontal size of this particular monitor cannot be adjusted to 
display the full frame. Hence, the left and right edges are not visible and can 
only be seen using the video copy processor. This is a serious shortcoming as 
edge artifacts caused by poor light placement cannot be detected readily. 


e Tektronix HC01 video copy processor. The video copy processor accepts 
composite video from the output channel of the frame grabber board and is 
used for making a hardcopy of the video image displayed on the composite 
video monitor. The video copy processor uses a thermal process to etch the 
video image onto 4-inch wide heat-sensitive paper. 


IOMEGA Bernoulli Box II dual removable cartridge disk (RCD) system. Each 
RCD consists of two flexible disks enclosed in a plastic housing providing 
20 MB of formatted disk storage. As each digitized image takes between 
100 kB to 250 kB in a compressed form, this allows storage of approximately 
100 images in each RCD. One drive is dedicated to storing image and data 
files. The data files contain the tabulated particle sizes (.dat files), the 
histogram results (.his files) and the session activities (.ses files). The other 
drive contains the executable SEMEX file and the C source codes for all the 
modules, together with the library files and the software development system. 
A 44 MB version of the RCD is also available. 


- HP LaserJet series II printer for producing a hardcopy of the results and the 
histogram plots (optional). 


SOFTWARE SUPPORT 
The software used to support this work include:- 


Microsoft C Optimizing Compiler, version 5.1, with Large Memory Model Run- 
Time Library, CodeView Symbolic Debugger and Overlay Linker. This 
provides the C programming language and the software development system for 
the SEMEX program. The compiler converts the C source files into object 
files which are linked together with the libraries by the linker. The large 
memory model library (LLIBCE.LIB) is recommended by Imaging Technology 
to support their ITEX PCplus software. This model allows for multiple 64 kB 
segments for both code and data to contain the large image and data arrays. 
CodeView aids in the debugging process by allowing single-step program 
operation while observing key variables. Microsoft also has a MAKE facility 
that automates the compilation and linking process through the use of a special 
batch file shown in Appendix A. 
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- Star Guidance Window BOSS windowing package. The Window BOSS is an 
extensive library of C functions for the creation, management and manipulation 
of text windows, which are essentially dialog boxes for program-user interaction. 
It allows for multiple windows to be created using dynamic memory allocation 
and manages these windows. Window attributes like border styles, border 
colors, and window foreground and background colors can be set. Context- 
sensitive help screens can be easily implemented. Forms and input functions 
are also available to aid the programmer in developing user-friendly dialog 
boxes. The Window BOSS is a shareware product. [Ref. 13] 


- Imaging Technology ITEX PCplus Large Memory Model (/TEXPCML.LIB) 
Library. The ITEX PCplus is a library of image processing and graphics 
functions specially written to support the PCVISIONplus Frame Grabber. Only 
a few of the functions are actually used in SEMEX and these include setting 
up of the frame grabber’s control registers and LUTs, reading and writing a 
single pixel value, reading and writing the frame memory from and to disk, and 
thresholding an image. Details can be found in the ITEX PCplus 
Programmer’s Manual [Ref. 14]. 


* Mathworks MATLAB 386. This mathematical programming package is used 
to produce the histogram plots subsequent to the ANALYZE stage. A script 
file was written to automatically read in a SEMEX data file and plot the 
results. It can also plot the Malvern data. This file appears at the end of 
Appendix B. 


* SEMEX (Scanning Electron Microscope EXtraction) Program is a set of 
program modules, written in the C programming language as part of this thesis 
work. This program comprises thirty-seven locally-written functions and 
handles the six stages required in extracting particle size information. The first 
involves setting up the system to ensure correct illumination, camera input gain 
and offset, and calibration. The second stage deals with the digitization of 
SEM images, cropping the digitized image to the desired area of interest, and 
complementing the image, so as to obtain a digitized image of the area of 
interest with particles dark against a light background. The third stage employs 
thresholding to form a binary image. The fourth and fifth stages involve 
particle tagging and particle sizing respectively. Here, particles are identified 
from the background and their x-chord and y-chord lengths and areas are 
measured. The final stage merges the data from several images and analyzes 
the resulting data. The analysis involves extracting three-dimensional features 
from the two-dimensional images. Presently, this last stage has been written to 
generate results similar to the Malvern MasterSizer. The latter puts out a 
histogram plot of percentage of particle volume against a logarithm scale of the 
particle diameters. SEMEX generates the histogram data and makes use of 
MATLAB to generate the histogram plots. The source listings for the various 


15 


modules are found in Appendix B and detailed descriptions are given in 
Chapter IV. 
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IV. DETAILED DESCRIPTION OF SEMEX 

SEMEX is the name given to the executable file made up of the seven € 
language program modules described in the following sections. SEMEX stands for 
Scanning Electron Microscope EXtraction. SEMEX was written with the aim of 
providing a user-friendly environment which is flexible, yet highly efficient for 
digitizing SEM photographic images and extracting particle size distribution data. 
Windows (dialog boxes) guide the user through the whole process of setup, 
acquisition, processing (more specifically clipping, tagging, and sizing), and analysis 
as shown in Figure 5. Generally, SETUP and ANALYZE need only be done once 
at the beginning and at the end of a session, respectively. ACQUIRE, CLIP, TAG, 
and S/ZE need to be done as many times as there are images to extract. 

Formatted inputs disallow illegal user responses while extensive error trapping 
prevents unintentional user inputs from resulting in fatal problems. A record of each 
session's activities is automatically created in an ASCII text file format which can be 
reviewed or printed out as a permanent record by any text editor program. The 
sessions file uses a filename formed from the current date (1.e., the first three letters 
of the month and a two-digit day) with the extension of .ses. In this way, each day's 
activities are recorded in a separate file. If more than one session is started in a 
single day, the second session will be appended to the end of the first. This prevents 
any record from being written over. Alternatively, the user can type in a filename 


during setup. 
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SEMEX | 
Invoked by typing 
This Is the entry program which calls ail SEMEX at the DOS 
the other modules. prompt. 


SETUP 
This needs only be 


This module handles the setting up of the run once for every 
camera, lights and program defaults. sesslon. 





ACQUIRE 


This module brings an SEM image into 
frame memory from disk or via a camera. 












CUE 


This module converts a gray scale Image 
into a binary image by a threshoiding 
process. 


Repeat these 
steps for the 
rest of the 
images. 


TAG 


This module performs object segmentation by 
assigning adjoining pixels of a feature with 
a unique fid value. 





SIZE 






This module uses the fid values to find 
the area, X-chord, and Y-chord of each 
feature which are then saved to a file. 


ANALYZE (with aid of MATLAB) 
This need only be 
This module merges selected data files and run after all images 
finds the particle volume distribution which from a rocket burn 
can then be plotted by MATLAB. have been processed. 





Figure 5. Typical sequence of activities in SEMEX. 
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Block diagram of SEMEX modules (rounded rectangles) and their 


functions. 


Figure 6. 
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Figure 6 gives an overview of the various modules and functions, and the 
filenames they are contained in. Detailed descriptions of each of the program 


modules, the algorithms used and any trade-offs made follow. 


A. PROGRAM MODULES 


1. Module MAIN 
The MAIN module consists of three functions main(), fginit() and 


session name(). Figure 7 gives the algorithms, written in pseudo-code, for the 


main() 
{ 
initialize frame grabber by calling fginit(); 
blank out image filename; 
call session_name() to create sessions filename from current date; 
display popup menu of SEMEX options; 
accept user response and call respective functions; 
return to display popup menu again; 
update sessions file with elapsed time; 
display sign-off message; 


) 


fginit() session_name() 
{ 

set hardware definition; call dos_getdate() to get current date; 
set frame dimensions; call dos gettime() to get current time; 
turn frame grabber on; convert month and day to filename string; 
initialize registers and LUTs; append file extension ”.ses" to string; 
select camera input, open session file using this string; 
clear frame memory; record current date and time; 
select and display frame memory; close session file; 
select input and output LUTs; return to main(); 
return to main(); 





Figure 7. Pseudo-codes for SEMEX main program main() (top), frame 
grabber initialization fginit() (bottom left), and sessions file 
naming function session name() (bottom right). 


functions in this module. All C programs must begin execution with the function 
main(). It ties together all the other modules and functions in a program, and this 


is no different for SEMEX. In addition, it initializes all the global variables to their 
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default values. These global variables are used by the various program modules for 
range checking and for customizing the SEMEX program to suit user needs. 

The function feinit() initializes the frame grabber by loading the correct 
addresses and configuration data. Once done, the frame memory is cleared and the 
input and output LUTs selected to prepare the frame grabber to receive digitized 
images captured with the vidicon camera. Function session_name() handles the 
formation of the sessions filename using the current month and day from the DOS 
system (obtained by calling dos getdate() function), and appends a .ses file extension 
to it. However, the user can specify a different filename when inside the module 
SETUP. The sessions file records all the session activities for that particular day, 


regardless of the number of sessions. Figure 8 shows a typical sessions file. 


Opening Session on 26 Feb 1991 at 18:09. 
CLIP: Image from Filename: sem02 . img 
Auto Threshold: 187 User Threshold: 187 
TAG: 273 features tagged in 14.2 seconds 
SIZE: Vertical scale used: 8.400000 pixels/unit length 
Min Length spec: 1 Max Length spec: 100 
Sized 273 Features within specifications 
Sizing took 10.4 seconds 
Conversion constants: Cx=0.141691 Cy=0.119048 Ca=0.016868 
AREA_M X-chord Y-chord 
5.803 12.044 4.167 
0.017 0.142 0.119 


ANALYZE: Merging data files 
Extracting data from SEM02.DAT 
274 extracted. Running Total is 273 
Total volume of minuscule particles is 3.726466 or 5.252 
Printing results to Feb26.his 
Printing MAT-file sem02.mat 
SEMEX was on for 2.4 minutes. 





Figure 8. A typical session recorded in the sessions file. 
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2. Module ACQUIRE 
The ACQUIRE module is actually a single function acquire() whose 
pseudo-code is shown in Figure 9. It accepts one of two sources of image input. The 


first makes use of a previously digitized image stored on disk. The second allows for 


acquire() 
{ 


open dialog box; 
if disk image desired 
call SEMIO function getim(); 
else 
display live video; 
if default selected 
use default gain; 
else 
interactively adjust gain; 
if default selected 
use default offset; 
else 
interactively adjust offset; 
snap a single video frame; 
transfer synchronization to frame grabber; 
deselect camera to prevent interference; 
crop left and right margins: 
complement digitized image if desired by calling ITEX PCplus function complement(); 
if image is to be saved 
append gain, offset and margins to comments; 
call SEMIO function putim(); 
update sessions file; 
close dialog box; 
return to main(); 





Figure 9. Pseudo-code for image acquisition function acquire(). 


acquisition of a live image from the vidicon camera. If SETUP has not been run, 
acquire() will allow the gain and offset of the camera input to be interactively 
adjusted to obtain a high contrast image. After positioning the object image in the 
field of view of the camera, a single frame can be acquired (a process called 
snapping a frame by Imaging Technology). The camera input is then disconnected 
by software to reduce the jitter (caused by synchronization conflicts between the 


vidicon camera and the frame grabber) which can corrupt the digitized image. The 
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acquired image is 512 pixels wide by 480 pixels high with each pixel represented by 
8 bits, thus allowing for 256 shades of gray. 

Once acquired, the digitized image can have its left and right margins 
cropped. This will eliminate any edge effects and textual information that could 
interfere with the extraction of particle data. Next, the image can be complemented 
to produce a final image where particle features are dark and the background light. 
This is carried out by an ITEX PCplus function complement() which performs a one's 
complement (one bits become zeros and vice versa) on every pixel in the image. 

Finally, the image can be saved to disk or left in the frame grabber's 
memory (frame memory) for further processing. If the image save option is 
exercised, program control is transferred to the putim() function in module SEMIO 
where a filename is requested. An extension of .umg is automatically appended. The 
user will be asked to supply a line of comments which will be stored in the image 
header, together with information on the gain and offset of the camera input, 
margins set and the vertical scaling factor. This information will be displayed 
together with the digitized image whenever it is recalled by the other modules. The 
stored margins are used by the modules to define the area of interest in the digitized 
image. This can result in an improvement in performance, as the whole frame need 
not be processed. The vertical scaling factor provides the S/ZE module with the 


conversion factor between pixels and some unit of measure (e.g., microns). 
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3. Module CLIP 


The CLIP module consists of four functions namely, clipmain(), clip(), 


autoclip() and findthd(). Their pseudo-codes are shown in Figure 10. 


clipmain() 
{ 
open dialog box; 
open session file, 
if disk image desired 
call SEMIO function getim(); 
else 
use image in frame memory; 
call autoclip() to threshold the image; 
call clip() to modify the threshold; 
update sessions file; 
close sessions file; 
display options and get user response; 
if image to be restored 
linearize output LUT; 
else 
map output LUT into frame memory; 
if image to be saved to disk 
call SEMIO function putim(); 
close dialog box; 
return to main(); 


} 


findthd() 
{ 
within each region 
determine pixel value; 
if pixel value is more than mid-gray 


compare it to existing threshold, 


take the minimum of the two and 
assign it as the new threshold; 
return to autoclip(); 


Figure 10. 





clip() 
{ 


start repetition 
display clipped image; 
interactively get threshold; 
end repetition if user satisfied; 
supply threshold to calling program; 


autoclip() 

{ 
define regions to be sampled; 
assume threshold is peak white initially; 
for each region 

call findthd() to find new threshold, 

use the minimum threshold for all regions; 
return to clipmain(); 


Pseudo-codes for clipmain(), interactive clipping clip() (top 


right), automatic clipping autoclip() (center right), and 
background threshold determination findthd(). 


Function clipmain() is called by main() and uses the other functions in this 
module to convert the digitized gray scale image to a binary two-tone image by a 
process called thresholding. The gray scale image can be retrieved from the disk, or 


the existing image in the frame memory, acquired previously, can be used. The 
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IEEE Standard Glossary of Image Processing and Pattern Recognition Terminology 
(IEEE Std. 610.4-1990) defines thresholding as 
The process of producing a binary image from a gray scale image by assigning 
each output pixel the value 1 if its corresponding input pixel is at or above a 
specified gray level (the threshold) and the value 0 if the input pixel is below 
that threshold. 


In mathematical terms, for a threshold T, the process can be represented 


as 
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where p(i, j) represents the gray scale value of the pixel at column i and row j. The 
gray scale value 0 represents peak black while the value 255 represents peak white. 

The purpose of thresholding is thus to produce a binary image where the 
features are black against a white background. This is critical to the subsequent 
Stages of processing. To ease the user workload, autoclip() and findthd() are called 
to automatically determine an estimate of the optimal threshold level. The algorithm 
assumes that the particle-to-background contrast is good and that the particles are 
dark against a light background. Given these two conditions, the algorithm examines 
six Small regions around the edges of the image and one more in the center to 
determine the background gray scale levels. The location and size of the regions 
have been determined empirically, and is based on the location of the lighting. Pixels 
with gray scale levels lower than mid-gray (pixel value of 128) will be skipped over 
as they are assumed to be part of a particle and not the background. The darkest 


of the background levels obtained, corresponding to the lowest gray scale value for 


25 


all the regions, is used as the threshold. Sampling over several regions helps to 
account for any non-uniformity in illumination. This algorithm tries to ensure that 
background pixels do not accidentally become converted to particle pixels. While the 
algorithm is not robust, it was found to work satisfactorily for images satisfying the 
above two conditions. The processing time required to determine the automatic 
threshold level is less than one second. Once the threshold level is determined, a 
ITEX PCplus library function, threshold(), is called to perform the actual thresholding 
process. 

As the automatic thresholding algorithm assumes a prion knowledge of the 
image, Manual intervention is allowed using the clip() function. This allows the 
current threshold to be interactively changed to fine tune the image. Changes made 
to the image are constantly updated to the video monitor. 

The resulting two-tone image can either be saved to a disk file for later 
processing, carried forward to the next stage of processing, or reverted to its original 
gray scale form to start over. Clipped files have a default file extension of .iml. 
Typically, the image is not saved but is kept in frame memory. A SETUP option 
allows the subsequent stages to be activated automatically to reduce user intervention 


and increase throughput. 


4. Module TAG 
The TAG module consists of the seven functions shown in Figure 6 earlier. 
Collectively, these functions determine whether a particle is made up of an isolated 


pixel or a group of adjacent pixels. The general term for this process of screening 
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the image for objects (features or particles) is called image or object segmentation. 
The IEEE Std. 610.4-1990 defines image segmentation as 


The process of dividing an image into regions for the purpose of object 
extraction. 


tagmain() tag() 
{ { 
open a dialog box; determine maximum feature size limits; 
open sessions file; Start timer; 
if disk image desired tag the first row of pixels by calling tagrow0(); 
call SEMIO function getim(); tag subsequent rows by calling tagrows(); 
else determine total features tagged; 


use image in frame memory; check search window size by calling checkmerge(); 
tag image by calling tag(); check for joined features by calling tagmerge(); 
close sessions file; display final feature count; 
if tagged image to be saved display time elapsed; 

call SEMIO function putim(); update sessions file; 
close dialog box; update frame memory to make process permanent; 
return to main(); return to tagmain(); 





Figure 11. Pseudo-code for TAG algorithm. Left column shows tagmain() 
while right column shows tag(). 


Image segmentation or tagging, as it 1s subsequently referred to in this 
thesis, is performed in two stages. Figure 11 shows tagmain() calling tag() to scan the 
clipped image previously stored in the frame memory (or retrieved from a user- 
specified .im1 disk file). Each non-background pixel is assigned with a feature 
identification number (fid). Isolated pixels are assigned unique numbers while 
adjacent pixels share the same number. Essentially, the fid number identifies pixels 
belonging to a single feature or particle. 

The actual tagging is done by two functions tagrow0() and tagrows() which 
operate on the first row and the subsequent rows in the frame, respectively. 
Figure 12 shows the pseudo-code for the two functions. Essentially, the functions 
Scan a row from left to right, skipping white pixels which are assumed to be part of 


the background. When a black pixel is detected, the pixel above is examined. If that 
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tagrow0() 
{ 
for each pixel in the first row of the image 
skip to next pixel if present pixel is background; 
if left pixel occupied 
adopt its fid value; 
else must be new feature 
assign new fid value; 
return to tag(); 


} 


tagrows() 
} 
for each row of pixels in the image 
for each pixel in the row 
skip to next pixel if present pixel is background; 
if left or above pixel occupied 
adopt its fid value; 
else must be new feature 
asSign new fid value; 
if subsequent adjacent pixels connected 
assign same fid value to these; 
return to tag(); 





Figure 12. Pseudo-codes for tagrow0() and tagrows() algorithms used for 
differentiating features from background in the first and 
subsequent rows, respectively. 


pixel is not white, it already would have been assigned an fid value and the current 
black pixel is rewritten with the fid of the pixel above. In addition, the preceding 
pixel in the row is also checked and, if it has been assigned an fid value, this is also 
changed. This effectively identifies adjacent pixels as being part of the same feature, 
as long as the features are convex polygons (that is, there are no jagged edges caused 
by clusters of particles). 

Because only 8 bits are allowed for storing a gray scale value, this allows 
for only 256 values. Allocating value 0 to peak black and value 255 to peak white 
leaves only 254 values. This means that only 254 particles can be uniquely identified. 
To overcome this limitation, the frame area is windowed into one or more 
rectangular regions. A simple function checkmerge() (see Figure 13) uses the initial 


number of features counted by tagrow0() and tagrows() to determine the dimensions 
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checkmerge() 
{ 
determine safe feature size within a group; 
if feature size is larger than safe size 
warn user; 
prevent larger features from being tagged; 
return to tag(); 


} 


tagmerge() 
{ 
for each row in the image 
asSign row search limits; 
for each pixel column in the image 
get pixel value of current pixel; 
get pixel value of pixel above; 
if either is background or 
if they are part of the same feature 
no action required, skip on to next column; 
else 
adjoining pixels have different fids and need to be merged; 
assign column search limits; 
for each row within the search limits 
for each pixel column within the search limits 
if pixel is part of the current pixel 
re-assign the value above to it; 
return to tag(); 





Figure 13. Pseudo-codes for checkmerge() and tagmerge() used for checking 
the number of features in a merging window and for carrying out 
the merging, respectively. 


of this sizing window. The formula used is given by 


480 (2) 


int (==) + 2 


SAFESIZE = int 





where SAFESIZE is the allowable size of the feature. The value 480 represents the 
number of rows in the frame and 255 is the number of unique fid values available 
plus one. The pseudo-function int() takes the integer portion of the quotient. For 
example, if the feature count is 254, the denominator will yield 2 and SAFESIZE 
would evaluate to 240 which is half the frame height. The sizing window is twice this 
value which implies that the whole frame is used. Features having a vertical length 


greater than SAFESIZE will be partitioned into two. If this happens, the program 
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will warn the user to rerun the TAG module after specifying that these larger 
features be excluded. Alternatively, the user may attempt to use a lower threshold 
level to reduce the number of features if the image is suspected to be noisy or the 
background is heavily textured. For any image, the maximum feature size should not 
be greater than a quarter of the image frame area (1.e., 256 x 240 pixels). Features 
within different regions may have the same fid number because, in addition, they are 
assigned a different group identification (gid) number to distinguish between them. 

After the first pass described in the preceding paragraphs, another pass is 
taken through the whole image to merge joined features. This is carried out by 
tagmerge() and is necessary because the first pass assumes that all features are convex 
polygons. However, if some of the particles are overlapping or have irregular cross- 
sections, then they would be assigned different fid numbers although they are joined. 
Consequently, these non-convex polygons have to be merged together and assigned 
a common fid number. The number of features merged will be highlighted to the 
user together with the number of features identified or tagged. The latter can be 
used for verification of the TAG algorithm with the SIZE module described in the 
next sub-section. Both should yield the same feature count. 

A limitation of the present algorithm is that border regions are ignored. 
This means that particles truncated by the edge of the image are treated as if they 
were complete particles. This may distort the actual distribution of particle sizes. 
This generally affects a small number of particles and is not expected to give rise to 


large errors in the result. However, should this error be deemed significant, the 
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usual practice is to define an area of interest somewhat smaller than the image such 
that a border exists. The border must be as wide as the largest particle. Particles 
lying on this border, and not truncated by the edge of the frame, will be tagged while 
those outside the border will be reset to the background color. In this way during 
SIZE, truncated particles will not be sized and hence do not affect the distribution. 
While this would yield more accurate results in terms of classifying particles, the 
smaller area of interest would yield a smaller sample population and more images 
would have to be used to compensate for this. In addition, it would not be practical 
to incorporate a border when particles could be as large as a quarter of the image. 

The tagged image can be saved to disk. It will have a default file 


extension of .im2. 


5. Module SIZE 

The SIZE module consists of four major functions namely sizemain(), 
size(), pixelsize(), and outdata(). Collectively, they perform the function of 
determining the area and physical dimensions of the particles. The main function 
sizemain() allows for a stored tagged image to be retrieved, calls size() to coordinate 
the sizing process, and, finally, calls outdata() to tabulate and save the results. The 
pseudo-codes for sizemain() and size() are shown in Figure 14. 

The function size() dynamically allocates storage arrays for particle areas, 
horizontal particle lengths (x-chords), and vertical particle lengths (y-chords). These 
arrays are indexed by the fid value of the feature whose dimensions are being stored. 


The function also keeps track of the processing time taken to size the image. Actual 
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sizemain() 
{ 
open a dialog box; size() 
open sessions file; { 
if disk image desired determine scale factor; 
call SEMIO function getim(); determine feature size limits; 
else check for sufficient memory; 
use image in frame memory; start timer; 
Size image by calling size(); size the image with pixelsize(); 
if results are to be saved display time elapsed; 
call outdata(); return to sizemain(); 
close sessions file; 
if sized image to be saved 
call SEMIO function putim(); 
close dialog box; 
return to main(); 





Figure 14. Pseudo-codes for SIZE module functions sizemain() (left) and 
size() (right). 


sizing begins when size() calls function pixelsize() to carry out the algorithm given in 
Figure 15. 

Starting with the top left corner of the image, each row of pixels is scanned 
for the fid numbers previously assigned by tag(). Three registers are used for 
tracking the current feature area, feature x-chord, and y-chord. These are initialized 
to zero whenever a feature is first found. When pixels are found with an fid value 
corresponding to the current feature being sized, these pixels are reset to peak black 
to prevent them from being recounted. At the same time, the area and chord 
registers are updated. When the end of a feature is detected (no subsequent rows 
having the same fid value as the current feature being sized), these registers are 
stored into the respective storage arrays for later retrieval by outdata(). Statistics like 
the largest and smallest features are also determined, together with any features 
being rejected for failing to satisfy the size limits specified during setup. These size 


restrictions can be used to filter out background noise or undesirably large features. 
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pixelsize() 
{ 
for each row in the image 
assign row search limits; 
for each pixel column in the image 
get pixel value of current pixel, 
if pixel value is background or has already been sized 
skip on to next column; 
else 
pixel is part of a feature yet to be sized; 
assign column search limits; 
for each row within the search limits 
for each pixel columm within the search limits 
if pixel is part of the current feature 
increment feature area; 
increment horizontal feature length; 
tag the pixel as sized; 
if horizontal feature length is zero 
end of feature is reached; 
go size next feature; 
else 
update maximum feature length; 
reset horizontal feature length; 
increment vertical feature length; 
save vertical feature length; 
save feature area; 
reset vertical feature length and area; 
determine min and max feature dimensions sized so far; 
determine smallest and largest feature sized so far; 
check if any size limits exceeded; 
display any rejects together with number of features sized; 
update sessions file; 
return to sizemain(); 





Figure 15. Pseudo-code for function pixelsize() which scans each pixel in the 
image and identifies it as part of a feature. 


The above procedure takes place inside a sizing window determined during 
setup (or modified by check _merge() as described in the earlier sub-section). The 
window must be as large as the largest pixel. This sizing window is moved from left 
to right, top to bottom over the whole frame. If the sizing window is large, a 
considerable amount of computation time is required to scan the area in the sizing 
window for pixels belonging to a particular feature. A performance enhancement has 
been incorporated by detecting the end of the feature and aborting the search for 
more pixels belonging to this feature. For small features, this results in a significant 


improvement in performance, as is shown in the next chapter. 
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Up to this point, all processing has been in terms of pixel dimensions to 
avoid floating point operations. For the vidicon camera used, a pixel is not square 
but has an X:Y aspect ratio of 1.2:1. To convert the pixel dimensions to microns 
requires the vertical scale factor which was determined during setup. The conversion 
is done in outdata(), shown in Figure 16. It uses the stored pixel values of the x- 
chords and y-chords and the total pixel area of each feature, and converts these to 
microns and square microns respectively. The results are tabulated on the computer 


monitor and saved to a .dat data file for subsequent analysis. 


outdata() 
{ 


open a dialog box; 

calculate conversion factors; 

request for data filename; 

open data file; 

display feature horizontal and vertical dimensions; 
display feature area; 


write the same results to data file; 

close data file; 

display min and max feature sizes for the whole image; 

display smallest and largest feature area for the whole image; 
update sessions file; 

close dialog box; 

return to sizemain(); 





Figure 16. Pseudo-code for function outdata() which tabulates the feature 
sizes and writes it to a data file. 
6. Module ANALYZE 
Many methods exist for analyzing particles which are normally distributed 
[Ref. 15]. However, few exist for multi-modal analysis. For this application, 
the algorithm is designed to collect all the data files specified by the user, to extract 
the area data from them, and then to calculate the equivalent diameter and volume 


of asphere. The end result of ANALYZE is a histogram of the sample population. 
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These steps are accomplished by the functions analyze(), merge_data(), 
extract_data(), and histo_vol(). The equivalent spherical volumes are used because 
a large particle has a much greater impact on the plume characteristics than many 
small ones. Referring to Figure 17, analyze() opens a dialog box and the sessions file 
and calls the other three functions. The function merge data() builds up a list of data 
filenames by prompting the user with the name of each data file in the current 
directory and asks whether each should be used. The function extract_data() then 
begins extracting the area from each of the data files. The user can choose between 
using the calculated area (AREA C) or the measured area (AREA M). The 
function histo_vol() takes the areas and calculates the equivalent spherical diameters 
and volumes. Although the particles may be irregular in shape, the use of equivalent 
spheres facilitates the plotting of histograms against a single size dimension (particle 
diameter). The function then outputs the data into a .his histogram data file. This 
file can be read out or printed with any ASCII text editor program. 

An option allows the user to save the same histogram data into a 
MATLAB-compatible .mat file. This makes use of a function called savemat() 
provided by MATLAB and has been locally modified for SEMEX. The histogram 
data can be plotted using a MATLAB script file called SEM.M. This file reads in the 
.mat file and plots a histogram of percent of total volume against a logarithmic scale 
of particle diameter. The listing for SEM.M is found in Appendix B. To enable 
comparison with the Malvern MasterSizer [Ref. 8], the same upper and lower limits 


for each bin is used. The only difference is that SEMEX has 38 bins as against only 
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analyze() 
{ 
open a dialog box; 
open the sessions file; 
call merge data() to merge data from different data files; 
call histo_vol() to calculate volume and histogram the result; 
close session file; 
close dialog box; 
return to main(); 


} 


merge data() 
{ 
use dos findfirst() to get first data file and its creation date; 
display file and date created and ask user whether to include this 
if user response is positive 
update session file with datafile name; 
allocate memory for the list; 
add filename to list; 
find the rest of the data files using _dos_findnext(); 
repeat the above steps until no more data files found; 
allocate memory for data array; 
call extract_data() to extract area from data file; 


} 


extract_data() 
{ 
for each of the selected data files in the list 
open the data file; 
read area into data array; 
repeat until end of file; 
close data file; 
erase the list; 


} 


histo_vol() 
{ 
determine bin limits; 
allocate memory for volume and diameter arrays; 
for each particle in the data array 
calculate the equivalent diameter assuming a circle, given area; 
calculate the equivalent volume assuming a sphere, given area; 
accumulate total volume; 
Sieve the volumes and collate into the correct bins; 
update sessions file with particles outside the defined sizes; 
print out results to histogram file; 
if desired, a MATLAB .mat file can be created for SEM.M to plot the histogram; 
de-allocate memory thus erasing all the arrays; 





Figure 17, Pseudo-codes for the ANALYZE module showing, from top to 
bottom, analyze (), merge_data(), extract data(), and histo_vol(). 


31 for the Malvern MasterSizer. The latter can only size down to 0.5 um while 
SEMEX has a resolution down to 0.125 um. SEM.M allows the results from the 
Malvern MasterSizer to be simultaneously plotted for comparison. These plots can 


be seen in the next chapter. 
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7. Module SETUP 
The module SETUP consists of two functions check equipment() and 
setup(). The first helps the user to set up the camera and light fixtures (refer to 


Figure 18), while the second is used to customize the SEMEX program. 


check_equipment() 
{ 

open dialog box; 

if changes required to equipment setup 
acquire live video; 
interactively set video input gain; 
interactively set offset; 
digitize a single video frame; 
synchronize to frame grabber for stability; 
deselect camera to prevent interference; 
call measure_line() to measure 5 micron reference line; 
repeat measurements as desired; 
complement image to get dark features on light background; 
update frame memory to record changes; 
call ITEX PCplus threshold() to threshold image at maximum level; 
call clip() to interactively threshold the image; 
repeat whole process if desired; 
open sessions file; 
update setup parameters; 
close sessions file; 

close dialog box; 

return to setup(); 


} 
Figure 18. Pseudo-code for the function check _equipment() which aids in 
the physical setup of the camera and the lights. 


To set up the camera and light table, the function check _equipment() turns 
on the frame grabber and begins acquiring live video through the camera. The user 
can manually focus the vidicon camera and set its aperture. At the same time, the 
gain and offset of the camera input can be interactively set, and the photographic 
image aligned. In order to fit the 4:3 frame aspect ratio of the camera, the SEM 
photograph has to be rotated onto its side such that the textual information is to the 
right. Having done this, a single frame of the image is digitized. 

Once acquired, measure line() is called to allow features on the digitized 


image to be repeatedly measured with the aid of a graphics cursor. The graphics 
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cursor, drawn and erased by put _cursor() and unput_cursor(), respectively, is moved 
around on the video monitor by using the arrow keys. The function chkkey() 
translates these keystrokes into x and y values which put line() and unput line() use 
for drawing and erasing lines. The lines are measured by calc line() and these 
measured lengths are recorded in the sessions file. Since SEM images have a 5 um 
reference line on the image, measurement of this line allows for size calibration of 
the system. The function measure line() assumes that this reference line is vertical 
and proceeds to calculate the vertical scaling factor to be used for converting pixels 
to dimensional lengths. Figure 19 shows the pseudo-code for this function. The 


graphics cursor and line manipulation functions are shown in Figure 20. 


measure line() 
{ 
display graphic cursor using put_cursor(); 
await pressing of cursor keys; 
remove graphic cursor using unput_cursor(); 
decode key pressed using chkkey(); 
if valid cursor key pressed 
display graphic cursor at new location; 
else 
repeat above sequence; 
await key pressed to get second position; 
remove graphic cursor using unput_cursor(); 
decode key using chkkey(); 
if valid cursor key pressed 
put graphic cursor at new location using put_cursor(); 
draw line to new location using put_line(); 
calculate length of line using calc_line(); 
if not zero length 
calculate scaling for a 5 micron reference line; 
open sessions file; 
update sessions file; 
remove line using unput_line(); 
remove graphic cursor using unput_cursor(); 
close sessions file; 
else 
go back to wait for second position; 
return to check _equipment(); 





Figure 19. Pseudo-code for the function measure_line() used to determine 
the vertical scaling factor of a 5 um line in the SEM image. 
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put_cursor() 


{ 
determine top- and left-most pixels, given the center; 
save pixel values under the graphic cursor, 
if center pixel value is greater than mid-gray 
draw black graphic cursor; 
else 
draw white graphic cursor; 
return to measure_line(); 
} 
unput_cursor() 
{ 
determine top- and left-most pixels, given the center; 
restore original pixel values thus erasing the cursor; 
return to measure _line(); 
} 
put_line() 
í 
if center pixel value is greater than mid-gray 
use black pixels to draw line; 
else 
use white pixels to draw line; 
if horizontal line is longer than vertical 
save horizontal pixel values under the line; 
draw a horizontal line; 
else 
save vertical pixel values under the line; 
draw a vertical line, 
return to measure_line(); 
) 
unput_line() 
{ 
if horizontal line is longer than vertical 
restore horizontal pixel values under the line; 
else 
restore vertical pixel values under the line; 
return to measure _line(); 
) 
chkkey ( ) 
{ 
test for valid cursor keys; 
if either arrow keys is pressed 
increment/decrement x or y respectively; 
if Home, End, PgUp or PgDn key is pressed 
step x or y respectively (step size is 10 by default); 
return to measure line(); 
} 
calc line() 
{ 
if horizontal line is longer than vertical 
return length of horizontal line to measure_line(); 
else 
return length of vertical line to measure _line(); 
) 
Figure 20. Pseudo-codes for functions put cursor() and unput cursor() 


which manipulate the graphics cursor while functions put_line() 
and unput line() deal with lines. 
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Next, the digitized image can be complemented (to produce an image 
where particles are dark against a light background) and clipped. Clipping is carried 
out within SETUP by calling clip(), described earlier. This causes pixels with the 
same or lower values to be displayed as black, thus allowing the user to visually 
gauge the uniformity of the background. Figure 21 shows two images with varying 


uniformity in illumination. 
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Figure 21. Samples of clipped images showing the non-uniformity of the 
illumination. Darkened areas have the same or lower pixel 
values. 


The left image is poorer than the right due to the poorer location of the 
lights. It is important that the illumination be evenly distributed so that the 
background intensities do not differ significantly over the image. One way of 
ascertaining this is to use a blank sheet of paper in place of the SEM image. The 
digitized image is then clipped with the threshold at maximum (threshold value 255). 
This should yield an almost black image initially. Decreasing the threshold will cause 
increasing portions of the image to turn white. The background should turn from 


completely black to completely white in about 50 levels of gray. If it takes more than 
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this or if the blacks appear blotchy, the illumination is uneven and needs to be re- 
adjusted. 

If the initial clipped image is not almost completely black, it indicates that 
the full dynamic range of the contrast has not been exploited. The offset may need 
to be changed and then this whole process repeated until the results are satisfactory. 

The second function setup() in the SETUP module allows the user to 
configure the SEMEX system, thereby streamlining the SEM extraction process. 
Figure 22 shows the pseudo-code for setup(). Various default parameters are 


setup() 
{ 


check camera and lighting by calling check_equipment(); 
open dialog box; 
display current default settings; 


accept user changes; 

display new default settings; 
close dialog box; 

return to main(); 





Figure 22. Pseudo-code for the function setup() in the module SETUP. 
displayed and can be modified to enhance the accuracy of the extraction. The 
processes are automated to reduce user workload. Several flags can be defined 
which cause program flow to be altered, thus changing the appearance of the 
program and the amount of user intervention. In addition, a sessions file can be 
specified to record Session activities and to store intermediate results, thus relieving 
the user from the need to write notes. A permanent record of each session can thus 
be kept. Figure 23 shows the SETUP dialog box with all the default parameters 


which the user can alter to tailor the SEMEX program. 
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SETUP DEFAULTS = 


GAIN LEVEL O Use Default[Y, 
OFFSEReLEVEL o 50 Use Default[Y, 
LEFT MARGIN aw) Use Default[Y, 
RIGHT MARGIN gees iy) Use Default[Y, 
Y~SCALE FACTOR : _ 1.000 Use Default[Y, 
Max Feature Size : 100 

Min Feature Size : l Use Defaults[Y,N]: Y 


Max Feature Count : 2000 Auto-Allocate Memory[Y,N]: Y 


N] 
N] 
N] 
N] 
N] 


ALL: Enable HELP screens [Y,N] 

SEMEX: CLIP, TAG and SIZE without asking 
ACQUIRE: Complement Image without asking 
CLIP: Load RAW Image without asking[Y,N] 
TAG: Load CLIPPED Image without asking[Y 
SIZE: Load TAGGED Image without asking[Y 
Session Filename: Febl9.ses 





Figure 23. SETUP dialog box showing default parameters which the user can 
change to customize SEMEX. 
8. Module SEMIO 

The SEMIO module consists of several globally-used functions getim(), 
putim(), and chkext(), which are called by most of the other modules. It performs 
such functions as reading a disk image file into frame memory, writing out an image 
stored in frame memory to disk, and checking the filename extension. Image files 
are stored in a compressed format to reduce the disk storage space required for each 
image. However, this has a small penalty on the time required to read and write an 
image file (one or two seconds). Figure 24 shows the pseudo-codes for the three 
functions. The functions make use of two ITEX PCplus functions readim() and 


saveim() to actually perform the image read and save operations, respectively. 
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getim() 
{ 
call chgext() to append default file extension to filename; 
prompt user with filename; 
get user response; 
get image and comments from disk using ITEX PCplus library function readim(); 
if error detected 
suggest trying again with new filename, 
display image on monitor; 
display comments in dialog box; 
extract margins and scaling factor from the comment line; 
return to calling program; 


} 


putim( ) 
{ 
call chgext() to append default file extension to filename; 
prompt user with filename, 
get user response and comments; 
append gain,offset, margin and scaling factor to comment string; 
save image and comments using ITEX PCplus library function saveim(); 
if error detected 
suggest trying again with new filename; 
return to calling program(); 


} 


chgext() 
{ 
determine the start of the file extension; 
discard old file extension; 
append default file extension; 
return to calling program; 





Figure 24. Pseudo-codes for the SEMIO module. Functions getim() and 
putim() handle image transfers to and from disk while chgext() 
is used to set the default file extension. 


The functions also handle the insertion of image parameters into the image 
header. Information like the gain, offset, margins and the scaling factor used in 
creating the digitized image are stored, together with any comments the user may 
wish to include. These parameter will be automatically retrieved when the image is 
brought in from disk. This ensures that no useful data is lost and it also reduces the 
setup time. Mundane tasks such as checking the existence of a filename and 


availability of disk storage space are also carried out by the functions in SEMIO. 
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V. RUNNING SEMEX 


A. SEMEX SETUP PROCEDURE 

SEM photographic images first have to be digitized and stored in the frame 
grabber’s memory before any processing or extraction of information can be started. 
The camera and lighting setup is critical for obtaining a good digitized image. Using 
a 28 mm focal-length lens, the vidicon camera has to be mounted at a height of 
approximately 13 inches from the SEM photograph on the light table. A transparent 
glass plate is placed over the photograph to keep it flat. After ensuring that the 
video cables have been properly routed, the IBM AT is powered up and SEMEX 
Started by typing SEMEX at the DOS prompt. A windowed menu like the one 
shown in Figure 25 will appear. Pressing a number will cause a corresponding 
selection to be highlighted. Alternatively, the arrow keys can be used to make a 
selection. 


Once SEMEX has been started up, the first thing to do is to run SETUP by 


atado ao dopado tato llo atan to toto oo o ota toto 
TC tt 


consists of five steps: 
1. Initial alignment and focusing. 
2. Setting camera input gain and offset. 
3. Determining the system scaling factor. 


4. Adjusting the illumination. 
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S EMEX 
SEM Extraction Program 
Naval Postgraduate School 


Setup SEMEX 
Acquire Image 
Clip Image 

Tag Features 
Size Features 
Analyze Features 


Use cursor keys to select 
Press [Enter] to execute 
Press [Esc] to quit 





Figure 25. SEMEX main menu showing the six options. 


5. Verifying and changing system defaults. 


The detailed procedures are given in the following sections. 


1. Initial Alignment and Focusing 

Running SETUP causes another window to appear over the main SEMEX 
window. A prompt will appear asking whether the user wishes to set up the camera 
and lighting. Pressing any key other than ‘N’ or ‘n’ will select the default answer 
of ‘Yes’ and this turns on the camera. A digitized image of the SEM photograph 
can now be seen on the video monitor. Rotate the photograph counter-clockwise 
such that the aspect of the photograph is the same as the video monitor, with the 
textual information on the right. Adjust the height and focus of the camera so that 


the maximum area of the photograph can be seen without borders. The Sum 
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reference line and part of the textual information on the photograph should also be 
visible on the right edge of the video monitor. Figure 26 shows the orientation of 


a typical image. 





Figure 26. Digitized SEM image showing the orientation of the Sum 
reference line and other textual information. 


Ensure that the image is in focus. Once this is done for one photograph, 
the height and focus of the camera need not be adjusted again unless the setup is 
disturbed or dismantled. The photograph, the lens, and the glass plate should be 


dusted to prevent dust particles from being digitized. 


2. Setting Camera Input Gain and Offset 
The program will next ask the user to adjust the gain of the camera input; 
initially, the gain is set at maximum (value 0). The user can interactively alter this 
gain setting or, alternatively, adjust the aperture of the lens. A setting of gain 
value 0 at f{/4 was found to be satisfactory. Pressing [Enter] moves the program on 
to adjust the offset. Again, the user can interactively change the offset or it can be 
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left at the default of 60. In general, this value should give the image on the video 


monitor good contrast. 


3. Determining the System Scaling Factor 

The vertical scaling factor can be determined by measuring the length of 
the 5 um line located at the top right edge of the video monitor with the graphics 
cursor (a small cross-hair) located nearby. The cursor can be moved using the arrow 
keys for single pixel movements in either of four directions. For faster movement 
of the cursor, the [Home], [End], [PgUp], and [PgDn] keys can be used to move left, 
right, up, and down respectively, in steps of ten pixels. Pressing [Enter], when the 
cursor has arrived at one end of the 5 um line, acquires the point. Next, move the 
cursor to the other end of the line and press [Enter]. This will cause the graphics 
line to disappear and the measured length in pixels will be displayed on the computer 
monitor, together with the vertical scaling factor (see Figure 27). If desired, other 
features can be measured. However, it should be noted that the scale factor from 
the latest measurement will be stored. This factor can also be changed in the SETUP 


screen. 


4. Adjusting the Illumination 
The placement of the lamps is crucial for obtaining even illumination. 
After visually determining that there is no glare or reflection from the lamps or 
overhead lighting, an image should be captured by pressing the spacebar. This 
causes the live video mode to be stopped and a single frame acquired. Press the 


[Enter] key if the image appears satisfactory or press the spacebar again to replace 
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ACQUIRE IMAGE - SNAP MODE 


Position cursor at first point 
using ARROW keys. Press [ENTER] when done 
Coordinates: X1: 456 VIO 


Stretch line to second point 


using ARROW keys. Press [ENTER] when done 
Coordinates: X2: 456 Y 2 2 


Feature is 42.0 pixels long. 
For 5 micron vertical line, 
Vertical Scale factor is 8.4000 pixels/micron 





Press Any key when done 


Figure 27. Measuring a line during SETUP. This determines the vertical 
scaling factor. 


the old image with another digitized snapshot of the photograph. Next, complement 
the image if particles are white against a black background. The program then 


continues by clipping the image at threshold value 254 as shown in Figure 28. 


ACQUIRE IMAGE - SNAP MODE 
Reduce to gauge Lighting Uniformity 


THRESHOLD LEVEL: 254 


Use [+] and [-] keys to adjust 


Press [ENTER] when done 


res: 


Figure 28. Adjusting the threshold level in SETUP to determine the 
uniformity of the illumination. 
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This should cause the whole image to appear almost black. If this is not 
the case, the offset may have to be changed. The user should use the [-] key to alter 
the threshold level until the background starts to turn white. If the photograph has 
been evenly illuminated, the background should turn white evenly. If the 
illumination is uneven, the background will appear blotchy. Press [Enter] to quit. 
This process can be repeated to re-adjust the lamp positions and camera input offset 
for optimal illumination and image contrast. An example of good even illumination, 


achievable with the current setup, is shown in Figure 26. 


5. Changing System Defaults 
The last step of SETUP is to verify and, if necessary, to allow the user to 
change the defaults that SEMEX and its modules will use. The defaults are shown 
in Figure 23. The user can sequence through the defaults, using the up and down 
arrow keys. If no changes are necessary, the user can simply press [Esc] at any point 


to terminate SETUP and return to the main menu. 


B. IMAGE ACQUISITION PROCEDURE 

If the gain and offset have been adjusted during setup, then there will be no prompts 
to adjust them here in ACQUIRE. If the user has elected to use a disk image, a 
dialog box like the one in Figure 29 will appear asking for the image filename. 
Upon retrieving the file, any comments embedded in the image header will be 


displayed. The program will prompt to allow the left and right margins to be 
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cropped by pressing the spacebar, so as to remove any artifacts or textual information 


(see Figure 30). Press the [Enter] key when done. 


IMAGE ACQUISITION 


READ IMAGE FROM FILE 
Filename: test . img 
COMMENTS: 
SEM 00001 from burn taken on 25 Sep 90 
Gain= 0; Offset= 60; Margins= 0,460; 


VSCALE= 8.400 


Press any key to continue 








Figure 29. Acquiring a stored image into the frame memory. When the 
image is fetched from disk, the comments stored in the image 
header will be displayed. 


IMAGE ACQUISITION aes 


noterar esn sese rere tere 
ts 


Press [SPACEBAR] to crop image 
one vertical line at a time 


Right margin: 465 


When done, press [ENTER] 





Cropping the right margin of the image in ACQUIRE. This is 
to remove the textual information on the image. 


Figure 30. 


If the features appear white against a dark background, press [Enter] to 


complement the image, else type ‘N’ or 'n' to leave the image unchanged. The 
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user can then choose to save the image in a disk file or leave it in frame memory. 
Choosing the former will bring the user back to the SEMEX main menu after 
supplying the image filename and any comments. The program will automatically 
insert the image parameters (gain, offset, margins and scale factor) into the image 


header. The sessions file is also automatically updated. 


C. IMAGE PROCESSING PROCEDURE 
Image processing inside SEMEX consists of three steps: 


1. A threshold is applied to the image to convert it into a binary image where 
particles are peak black and background pixels are peak white. 


2. The image undergoes object segmentation whereby objects (i.e., particles) are 
discriminated from the background. Each object is given a unique identification 
number in the form of a gray scale value. 


3. Each object is sized by measuring all pixels with the same identification 
number and converting to microns. 


The details of these procedures are given below. 


1. Thresholding with CLIP 


The thresholding function is selected by typing 3 on the main menu which 


OOOO 
tota 


To threshold the digitized image already on the video monitor, the user simply 
presses [Enter] at the prompt ‘Read from disk file [N]?’. The program then begins 
its automatic threshold level determination. Almost immediately, the clipped image 
is displayed on the video monitor. The user can choose to vary the threshold (to fine 
tune the thresholding process) by using the [-] or [+] keys (see Figure 31). The 
threshold value is displayed, as well as recorded in the sessions file. A copy of the 
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clipped image can also be saved to disk, although this is not necessary (see 


Figure 32). The default file extension is .im1 for a clipped image. 


CLIPPING IMAGE 


Determining threshold... 
THRESHOLD LEVEL: 171 


Use [+] and [-] keys to adjust 


Press [ENTER] when done 








Figure 31. Clipping image using the automatic threshold which the user 
can subsequently change using the [+] and [-] keys. 


CLIPPING IMAGE 


1: RESTORE image to original and abort 
2: EXIT without saving 
3: SAVE the modified image 


Select option by NUMBER [3]: 


OS 
retos 


ta ere ete ate 


Figure 32. Screen for selecting the outcome of the clipped image. Option 
3 is the default which can be invoked by pressing [Enter]. 
2. Image Segmentation using TAG 
Image segmentation is performed by the TAG module. This module may 


be automatically invoked after CLIP completes, or it can be selected by pressing 4 


JZ 


first scans the clipped image and tags each pixel with an identification number (fid). 
Isolated pixels have unique numbers while adjacent pixels share the same number. 
Essentially, the fid number identifies pixels belonging to a single feature. The second 
pass searches for joined features with different fid numbers. These are merged, 
leaving only one fid number for all pixels belonging to that joined feature. Figure 33 


shows the TAG screen after it has completed tagging an image. 


== TAGGING IMAGE 
TAGGING FEATURES in Progress 
190 
Largest permissible feature is 240 pixels 
Combining joined features... 


FEATURE COUNT: 184 


Elapsed Time: 16.0 seconds 





Save image to Disk File [N]? 





Figure 33. TAG screen showing the number of features, the number 
merged and the final count. The final prompt requests 
whether to save the tagged image. 


3. Feature Sizing using SIZE 


To run the SIZE option, the user selects 5 to highlight the option 5, Size 


seest 


a oro, 
PeT E o* 


scan the tagged image in frame memory or on disk, and measure the pixel 
dimensions of each feature. If any feature fails to satisfy the size limits imposed 
during setup, these will be flagged out to the user. Figure 34 shows the information 


on the SIZE screen. 
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SIZING FEATURES 


Sized 184 Features within specifications 


Elapsed time: 8.3 seconds 


Save data [Y]? 





Figure 34. SIZE screen showing number of features sized and the time 
taken to size them. There is also a prompt to save the sized 
feature data. 


If the sizing is successful, the user can save the results into a data file 
which has the same filename as the image but with the extension .dat. The results 
are also tabulated on the computer monitor twenty features at a time (see 
Figure 35). The largest and smallest feature dimensions are also recorded in the 


sessions file. 


D. IMAGE ANALYSIS PROCEDURE 


1. Merging the Data Files 
Analysis of the images is usually carried out after a series of SEM images 
from the same firing have been sized. Each sized image will generate a data file 
containing areas and chord values. In order to obtain a Statistically significant 
sample population and thereby reduce sampling errors, the data files from all the 
images belonging to one rocket motor burn have to be merged. Module ANALYZE 


handles this merging operation and is invoked by selecting 6 at the SEMEX main 
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TABLE OF FEATURE DATA 
AREA C AREA M xX-Chord Y-Chord 
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Figure 35. Tabulated data showing the calculated equivalent elliptical 
areas (AREA C), the measured pixel areas (AREA M), the 
X-Chords and Y-Chords for 20 features. 


menu 6. Analyze Features. The module will begin by listing all the .dat data files 


in the current directory one at a time, for the user to select (see Figure 36). 


2. Calculating Particle Volume 
ANALYZE calculates the particle volume by assuming that the particles 
are spheres. It takes the measured particle areas from all the selected data files and 
determines the equivalent particle diameters and volumes. A 38-bin histogram is 
tabulated and displayed 10 bins at a time. This result is also saved into an ASCII . his 
file and can also be saved in a MATLAB-compatible .mat file. If the latter is 
chosen, the user can quit SEMEX and have the histogram plotted by running SEM.M 
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MERGING FILES 


DATA FILENAME DATE CREATED 


TESTÍ. DATE leI FEDNI 


include (Is 





Figure 36. ANALYZE screen prompting one data file at a time with its 
date of creation. Pressing ‘Y’ or ‘y’ accepts the datafile. A 
count is kept of the number of files selected. 


inside MATLAB. The histogram plot shows percentage of total particle volume 
against a logarithmic scale of particle diameter. This plot format is similar to that 
put out by the Malvern MasterSizer. For ease of comparison, both SEMEX and 


Malvern data can be plotted. Sample plots are shown in the next chapter. 
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VI. EXPERIMENTAL RESULTS 
This chapter describes the preliminary experimental results obtained. It is 
divided into three sections: system calibration and error quantification, comparison 
of the performances of SEMEX and HOLOGRAM, and, correlation of the results 


with that of the Malvern MasterSizer. 


A. CALIBRATION 


1. Determining Pixel Aspect Ratio 
In order for outdata() in SIZE to properly convert pixel counts into 
dimensional lengths and areas, there is a need to first determine the pixel aspect 
ratio. In the vidicon camera used, the aspect of a single pixel is not square. To 
determine the actual aspect ratio, a three-inch rule was placed, first horizontally, then 
vertically. The digitized lengths were then measured using the measure line() 
function in SETUP. The measurements are tabulated in Table I. 


Table I. CAMERA PIXEL SIZE 


3-inch Reference length 350 pixels 420 pixels 
Pixel Length (inch/pixel) 0.00857 0.00714 






The pixel aspect ratio can be obtained by taking the ratio of the horizontal 


length to vertical length of a single pixel. This is equal to 420 : 350 (or 1.2) since the 
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same 3-inch length was measured. In SEMEX, this constant is defined as the 


ASPECT RATIO. This works out to be the same as the frame aspect ratio. 


2. Quantifying System Errors 
To quantify the system errors, a test pattern was created by placing on a 
white background, 48 black etch-resistant circles used in printed circuit artwork 
(commonly called donut pads). Each circle has an outer diameter of 0.187 + 0.003 
inches and an inner diameter of 0.062 + 0.003 inches (see Figure 37). The variability 


of the outer diameter gives rise to an error of 1.6%. 





Figure 37. Calibration test pattern consisting of 48 donut pads placed in 
four rows. 


The gain, offset and illumination were adjusted in accordance with the 
setup procedure and precautions given in Chapter V. The measure line() function 
was used to determine the vertical outer diameter of the donut pad, in pixels. The 
vertical scale factor, VSCALE, could then be determined by taking the ratio of the 
measured pixel length (28 pixels) and the true vertical length (0.187"). For the setup 
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used, VSCALE was found to be 149.733 pixels/inch. Since measure line() only 
determines VSCALE correctly for a5 um line, VSCALE has to be manually entered 
in the SETUP screen. SEMEX calculates the horizontal scaling factor by dividing 
VSCALE by ASPECT RATIO. The area represented by a single pixel can then be 
found by taking the reciprocal of the product of the horizontal and vertical scaling 
factors. This area, multiplied by the number of pixels covering a particle, gives the 
measured area of that particle. 

Next, the ACQUIRE module was executed to digitize the test pattern into 
frame memory. CLIP was then invoked and the automatic thresholding function 
autoclip() applied. This yielded a visually-clean binary image for a threshold of 157. 
After exiting the CLIP module, TAG and SIZE were invoked. The results are 
tabulated in Figure 38. The maximum percentage error is given by the maximum 
deviation from the mean expressed as a percentage of the mean. 

The elliptical areas (labelled AREA C in column 2) were calculated by 
taking the products of the x-chords and y-chords (these correspond to the major and 
minor axes of an equivalent ellipse) and multiplying by the constant, 7/4. The 
AREA C values corresponded reasonably well with the actual pixel areas measured 
(labelled AREA M in column 3). The difference of the two mean areas amounted 
to 3.9% of the mean of AREA M. 

Although the effect of the non-square pixel aspect ratio had been 
compensated for, the results showed that the x-chord and y-chord lengths were still 


not equal. The difference between the two mean chords was 9.4% of the mean y- 
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Mean 
Max % 


Std Deviation 


Figure 38. 


ID AREA C 
i 0.029 
2 0.029 
3 0.026 
4 0.026 
5 0.027 
6 0.028 
7 0.028 
8 0.027 
9 0.027 

10 0.024 

11 0.024 

12 0027 

13 0.028 

14 oo 

15 YO 

16 0 OT 

17 0.026 

18 OMOT 

19 0.027 

20 0.027 

21 0.024 

22 0.025 

23 0.026 

24 0.024 

25 0.025 

26 0.026 

Dy 0.027 

28 B02 7 

29 OA 

30 0.027 

31 0.028 

32 0.028 

33 0.027 

34 0.027 

35 0.024 

36 0.025 

37 0.024 

38 0.024 

39 0.025 

40 0.026 

41 0.026 

42 0.028 

43 MOZ 

44 0.028 

45 0.028 

46 0.027 

47 0.026 

48 0.027 

0.027 
Error 9.43 
0.001 


Output from SIZE showing the results obtained from the 
calibration test pattern. The calculated area, AREA C, is 


AREA_M 
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X CHORD 


0.184 
0.184 
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0.184 
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0.160 
0.160 
0.176 
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DAT 
SAS 
ONIS 
0s T76 
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0.168 
0.168 
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0.184 
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0.160 
0.168 
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0.168 
On 76 
DIS 
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OLAS 
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given by 7/4 * X CHORD * Y CHORD. 
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Y CHORD 


0.200 
0.200 
0.194 
0.194 
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0.194 
0.194 
0.194 
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0.194 
0.194 
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0.194 
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0.194 
0.194 
0.187 
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0. Taz 
0.187 
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0.187 
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chord. The deviation within each chord amounted to a 8.2% error for the x-chord 
and a 3.7% error for the y-chord. Taking into account the variability of the donut 
pads (1.6%), the system has contributed to errors of 6.6% and 2.1% for the x-chords 
and y-chords respectively. 

In an attempt to identify the source of this error, a square grid was 
digitized. A rectangular graphics box was then placed over the grid lines. lt was 
found that vertical lines near the right edge were slanted while lines along the other 


three edges were correctly digitized (remained orthogonal). Figure 39 shows the 
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Figure 39. Digitized image of a square grid showing that vertical lines 
near the right edge were slanted when compared with the dark 
rectangular graphics box. 


square grid image with the artifact enhanced by the rectangular box. The maximum 
divergence was found to be three pixels. Using a plumb line and a liquid level, the 
light table and camera were checked for squareness and were found to be accurately 


aligned. It is suspected that this error is due to variations in the horizontal scan 
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velocity of the vidicon camera down the frame. If this is true, then the vidicon 
camera has introduced a scanning error of three pixels. However, this amounts to 
only 0.6% error over the whole frame length and is insufficient to account for the 
deviation in the chord values. 

Next, a 3-dimensional plot was made by sampling the pixel values around 
a single donut pad (32 x 32 pixels) in the calibration test pattern. Each grid square 
represents one pixel. The pixel values are represented by the height (z-axis) of the 


plot. Figure 40 shows that the edges of the donut pad are not vertical. This suggests 





Figure 40. 3-D plot showing the pixel values from a digitized image of a 
donut pad. Note that the plot has been inverted for clarity. 
This gives rise to high peaks for dark regions. 


that the analog-to-digital (A/D) converters in the frame grabber are not fast enough 
to digitize a sharp transition from peak white to peak black (the latter is represented 
by the flat circular region on the plot). The contrast transfer function of the camera 


and lens system also contributes to this limitation. The slope represents the amount 
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of error. From the plot, the change from peak white to peak black takes one to two 
pixels. As there are two slopes for each feature, the error could be as much as 4 
pixels. Depending on the threshold selected, the cross-sectional area would change. 
The worst case slew rate error for the case of the donut pad would be 4 pixels out 
of 28 (the measured pixel length), or 14.3%. The error is a function of the size and 
contrast of the feature. This error is very significant for small features of the order 
of one to five pixels. For the magnification used in the SEM, a 1/8 um particle 
would occupy only two pixels and have a slew rate error of up to 33.3%. As the 
particle size increases, the error drops dramatically. A 1 um particle, for example, 
would have a maximum slew rate error of only 6.3%. 

A major problem faced is the determination of the proper threshold to use. 
In chapter V, an automatic thresholding algorithm was described. The use of an 
automatic threshold can significantly improve the throughput. However, the criteria 
for threshold determination is an important consideration. Figure 41 shows a three- 
dimensional plot of a 32 x 32 pixel gray scale image. The left portion shows the 
varying gray scale levels of each pixel. When the region is clipped, the gray scale 
levels below the threshold are suppressed. However, different thresholds result in 
vastly different results as can be seen by the center and right portions of the figure. 
The right portion was obtained using the automatic thresholding algorithm. From 
the figure, it can be seen that the limited slew rate has prevented small features from 
exhibiting a high pixel value. Hence, they tend to be suppressed. In this particular 


case, the autoclip() function was able to bring out thirteen features. 
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Figure 41. Three-dimensional plot showing a 32 x 32 pixel region drawn three 
times at different thresholds. The left portion is the gray scale 
image. The center and right portions are at thresholds of 33 and 
165. 


B. PERFORMANCE COMPARISON WITH HOLOGRAM PROGRAM 
The HOLOGRAM program written by Hockgraver [Ref. 6] consists of the 
following modules: 
1. Filtering routines. 
2. Image threshold. 
3. Feature identification. 


4. Feature sizing. 
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Of these four modules, only the last three modules have similar functions in 
both HOLOGRAM and SEMEX. The first module (filtering routines) is used to 
reduce the effect of speckle introduced during image reconstruction from a hologram. 
This is not applicable to the SEM images as there is no speckle. In addition, 
SEMEX has two essential stages not carried out bb HOLOGRAM. These are the 
acquire and analyze stages). HOLOGRAM depended on the ImageActionplus 
software to perform image acquisition and on Statgraphics for analysis. This 
generally took much longer, as the ImageActionplus was not tailored to perform the 
acquisition, cropping and complementing of the images in an efficient manner. Also, 
the scaling factors had to be manually determined. Similarly, Statgraphics required 
its Own setup procedure. 

Comparison of the performance of the three modules common to SEMEX and 
HOLOGRAM was carried out by subjecting both programs to a common set of 
images. These images were previously acquired and then stored on disk. The images 
had different numbers of features ranging from 48 to 920, and are shown in 
Figure 42. The execution times are tabulated in Table II and the speedup calculated. 


Speedup is defined as 


Execution Time of Slow System (HOLOGRAM) (3) 


Speedup = 
p Es Execution Time of Fast System (SEMEX) 


and is a standard measure for the performance improvement of two systems. 
From the results, it can be seen that the speedup becomes more significant as 


the images become more complex. Almost an order of magnitude improvement has 
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Image +1 








Figure 42. 


Image #3 Image #4 


Images used for testing the performance of SEMEX against 
HOLOGRAM. Image #1 Is the calibration test pattern with 48 
features while images #2, #3 and #4 are actual SEM images with 
177, 495 and 920 features respectively. 
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Table Il. COMPARISON OF EXECUTION TIMES FOR HOLOGRAM 
AND SEMEX. 


|| Procramruncrion | 41 | 42143 | 42 
preamp |a| sef 06 | iss 


| SPEEDUP (H3/53) | 66 | 152 | 255 | 244 
SEMEX Noea) | SO | 66] 78 | 108 








1. The times for steps H1 and S1 do not include times to input filenames (approximately 
10 s each). The H1 times are based on one iteration only. In practice, 3 to 5 iterations 
may be required before arriving at a satisfactory threshold. For S1, the thresholds are 
determined automatically by autoclip(). These same thresholds are used in THRESHIT 
so that the same binary images are processed by both programs. 


2 The times for steps H2 do not include times to input filenames. For S2, no filename 
entry is required as the frame memory is used. 


3. The times for steps H3 do not include times to input filenames and answering prompts 
(approximately 25 s). For S3, no user intervention is required. 


4. The times for H4 and S4 includes all the user input and setup times. 
a: Four images (#1 to #4) were used with 48, 177, 495 and 920 features, respectively (see 
Figure 42). 
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been achieved. This is because the four nested loops used in both THRESHIT and 
SIZEIT (in the HOLOGRAM program) become highly inefficient. In SEMEX, this 
inefficiency is prevented by prematurely aborting the loops whenever the end of a 
feature 1s detected. 

Another observation was that, for the same image and threshold, the modules 
within the HOLOGRAM program produced different feature counts. For example, 
with image #2, THRESHIT counted 177 features (same as TAG and SIZE) but 
SIZEIT counted only 175 features. Only image #1 produced consistent results for 
HOLOGRAM, whereas SEMEX produced consistent counts for all the four images. 
This suggests that one or more of the HOLOGRAM modules may not have been 


correctly coded. 


C. CORRELATION WITH MALVERN MASTERSIZER 

From one of the recent test firings using 4.69% aluminum solid-propellant, two 
sets of SEM images were extracted and analyzed. The results were then compared 
with that from the Malvern MasterSizer. The two sets of results are plotted in 
Figure 43 and Figure 44. 

The results show that as the particle counts increase, the distributions produces 
a better correlation with the Malvern. Hockgraver [Ref. 6] showed that about 1,000 
particles were required to produce a Steady-state distribution for hologram images. 
From this particular set of images, 1,062 particles appeared to be insufficient. 
Unfortunately, there were no more SEM images from the same burn to extract. 


Hence, it is not possible to verify the minimum sample size requirement at this time. 
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Figure 43. MATLAB plot showing the histogram data from SEMEX and 
Malvern MasterSizer. The SEMEX data was obtained from three 
images with a total of 584 particles. 


It is suspected that an exact value cannot be determined even if sufficient SEM 
images were available. This is because the particle distribution varies considerably 
with the location of the probe tip (in both the radial and longitudinal directions), the 
portion of the filter paper from which the SEM images were taken, and the 
propellant characteristics. There is also a high possibility of debris and other 
contaminants being collected on the filter paper and this could skew the resulting 


particle distribution. 
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Figure 44. MATLAB plot showing histogram data from SEMEX and Malvern 
MasterSizer. The SEMEX data was obtained from six images with 
a total of 1062 particles. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 
This chapter lists the limitations of the present SEM extraction system and 
makes recommendations on how to further improve its performance. The 
recommendations are divided into three sections: hardware, software and 


methodology, respectively. 


A. HARDWARE LIMITATIONS AND RECOMMENDATIONS 


1. Light Table 

The light table setup can be improved by using a more uniformly 
distributed light. The use of a fluorescent ring light with the camera lens placed 
though the center will help distribute the light uniformly on the photographic image. 
Other external illumination (from windows and ceiling) should be excluded. Use of 
diffusers would further reduce any uneven distribution of the illumination. A rigid 
photograph mount or stage would also be helpful for accurate alignment of the 
images. This would allow for the use of image subtraction techniques whereby an 
image of the blank background is first taken. This image is then subtracted from the 


actual SEM images to remove any dirt or artifacts present in both images. 


2. Video Monitor 
The present video monitor is unable to display the full image frame. 


Consequently, any edge artifacts cannot be seen. The use of a video monitor that 
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has adjustable vertical and horizontal static convergence controls (V-SIZE and H- 


SIZE) would enable all the frame borders to be seen. 


3. Video Camera 

The use of a Charge-Coupled Device (CCD) camera with square detector 
elements would overcome the problem of the pixel aspect ratio not being square and 
also eliminate the errors due to the scanning of the vidicon camera. In order to 
improve the accuracy of the extraction, higher resolution CCD cameras (typically 
1024 x 1024 pixels) would be required. Unfortunately, these cameras are not RS-170 
compatible and cannot be used with the existing frame grabber. The contrast 
transfer function (CTF) of the camera and lens system may have to be further 


investigated to determine its contribution to the slew rate errors. 


4. Direct Acquisition of SEM Images 

Presently, the Hitachi S450 SEM has a high persistence video display for 
the operator and a high resolution line scan for exposing the photograph. The 
problem of providing uniform illumination during acquisition of the SEM 
photographic images may be overcome by acquiring images directly off the SEM 
video display. However, the video display resolution is significantly poorer than the 
line scan. Hence, it is not advisable to record SEM images directly off the former. 

A better arrangement would be to tap out the signal direct from the SEM 
line scan. The line scan takes approximately five seconds to cover the whole frame 
area. To use this signal would require building an interface with the correct video 


impedance matching and timing synchronization. The feasibility of this approach 
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would depend on the availability of technical information on the SEM. However, by 
mounting the video camera in place of the Polaroid camera, it may be possible to 
integrate all the frames (approximately 300, based on 5 s at 30 frames/second) to 
form a complete image. 

An ideal solution would be to use a SEM that has a RS-170 video output 
and a programmable stage. In this approach, programmed instructions can be 
inserted to make the SEM stage scan all the areas on the sample without overlap. 
The images captured can be processed in real time if a faster processor is available. 


This will reduce the sampling errors due to small sample population. 


5. Sun Workstation 

Various optimization techniques have already been employed to enhance 
the speed of the system. This has resulted in an order of magnitude increase in 
speed over that of the HOLOGRAM program. Further attempts to increase the 
speed would require disproportionate amounts of effort (Amdahl’s Law). Presently, 
the processing speed is limited by the IBM AT/386 and the frame grabber memory. 
In addition, the resolution of the current setup allows reliable sizing down to only 
1/8 um (with some degradation in accuracy). This resolution is limited by the video 
camera and the frame grabber. 

By upgrading to a Sun SparcStation 1, it is expected that another order of 
magnitude increase in speed would be possible with a two-fold increase in resolution. 
However, the frame grabber, video copy processor and monitor would also have to 


be replaced. The Sun SparcStation 1 has a 32-bit SPARC CPU running at 1.25 MIPS 
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coupled with a floating-point unit for faster floating-point computations. The higher 
integer and floating-point performance over a PC would auger well for processing the 
larger image arrays. It has a 104 MB internal SCSI hard disk and a 3.5" 1.44 MB 
floppy disk drive for easy data transfer and DOS compatibility. The latter would 
allow the SEMEX program to be ported over. A higher-resolution frame grabber 
with function calls compatible with the PCVISIONplus should be used to minimize 
software portability problems. The monitor resolution is 1152 (h) by 900 (v) pixels 
with a pixel aspect ratio of 1:1. The monitor has horizontal and vertical static 


convergence controls to adjust the frame size. 


B. SOFTWARE ENHANCEMENTS 


1. Automatic Camera Input Adjustment 

Automated setup procedure for adjusting the camera inputs can be done 
by generating a histogram of the distribution of pixel gray scales. An image with 
good contrast would have a good spread of gray scales ranging from peak black to 
peak white. An algorithm could be devised that would analyze the gray scale 
distribution and would adjust the camera offset automatically to maximize image 
contrast. This would eliminate the subjective determination of the camera input gain 
and offset. The setup time could also be shortened, as there would be less user 


intervention required. 
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2. Automatic Threshold Algorithm 
Currently, the background is sampled at fixed points which have been 
empirically-determined, based on the existing lighting arrangement. Changing the 
lighting could affect the performance of the automatic threshold algorithm. 
Consequently, the regions may not be optimal and may need to be changed. With 
faster processors available, the whole frame area can be scanned to determine the 
background threshold. This would result in a more accurate background threshold 


assessment. 


3. Improved Image Processing Algorithms 

Improved image processing, such as independent sizing of overlapping 
particles, filling in of hollow areas, and filtering off of irregularly-shaped features that 
fail certain roundness and sphericity tests, could be added. Currently, joined features 
are treated as one large feature with an irregular perimeter. However, if the 
curvature of each of the overlapped feature could be found, the outline of each of 
the overlapped features could be determined. This requires a significant amount of 
computation and, hence, was not implemented in SEMEX. 

Up to this point, the only sphericity test carried out is by calculating the 
area of an equivalent ellipse using the x-chord and y-chord lengths. This could be 
improved by measuring other parameters, such as the perimeter or radius of 
curvature. In some cases, particularly when large particles are involved, it was noted 
that the centers of the particles were not filled. This would result in smaller 


measured areas (AREA M) and could introduce false features within the hollow of 
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the particle (see Figure 45). TAG may have to be modified to test for the existence 


of unfilled areas and to fill these up automatically. 
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Figure 45. Three-dimensional plot of a particle showing that the center 
region is not a plateau. A smaller particle can be seen rising 
out of the center. This would give rise to two particle counts. 





4. Frame Border 
The addition of an area of interest may be justified if particles are likely 
to be small. This would then give rise to a border region whereby particles cutting 
the border would be sized and counted as a fraction of a particle, depending on the 
proportion inside the area of interest. Particles outside the area of interest are not 
sized. This will eliminate sizing errors due to incomplete particles. An alternative 


scheme would be to simply exclude all particles touching the edge of the frame. 
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C. IMPROVEMENTS IN METHODOLOGY 


1. Photographing SEM Images 

As far as possible, images should be taken with the same magnification 
and brightness and contrast settings. This would enable the acquisition stage to 
proceed faster. SETUP would have to be called only once and all the parameters 
would then be left unchanged for the whole set of images. 

From the five sets of images processed, it was found that up to four 
exposures could be made on one photograph without significant degradation in the 
results. This saves time in the photographing process and also in the extraction 
process. The only limitation is where overlapping of particles begins to occur, due 
to the multiple exposures. This will result in larger particles being sized. 

In chapter VI, more than 500 particles were used in the histogram plots. 
This was found to be insufficient. Hence, there should be sufficient SEM images to 


extract so that a steady-state distribution can be obtained. 


2. Running SEMEX 
Although SEMEX allows the user to customize the sequence of the 
processing, there is a preferred sequence to extracting SEM images, namely: acquire, 
clip, tag and size. This was shown in Figure 5 and the aim is to reduce disk reads 
and writes which takes between 5 and 8 seconds per image. An image which has 
been acquired and saved has its pixel values stored in the frame memory. Hence, 
CLIP can be initiated without recalling the image from disk. Similarly, after each 


clip and tag operation, the frame memory can be used by tag and size, respectively. 
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In this way, there is no need to retrieve an image from disk. Another advantage is 


that intermediate disk files need not be maintained (i.e., .im1, .im2 and .im3 files). 


D. CONCLUSIONS 

The IBM AT/386 running SEMEX has been found to speed up the extraction 
of particle sizes from SEM images. The results show that for the three stages of 
clipping, tagging and sizing, a speedup of nearly an order of magnitude was possible 
over that of the HOLOGRAM program. The speedup and throughput would be 
considerably more if all the stages of acquisition, clipping, tagging, sizing and analysis 
were compared. The errors obtained for small particle sizes on the order of half a 
micron or less are considerable. To reduce these errors, the magnification of the 
SEM images would have to be increased or a higher-resolution acquisition system 
would be required. The simplest approach would be to increase the magnification. 
This would require more photographic images and consequently more operator time 
and manpower cost. A higher-resolution system would entail an increase in 
computational cost on the order of the square of the frame length in pixels. This 
would strain the IBM AT/386 but would be comfortably handled by a Sun 
SparcStation. The capital outlay could be justified by the gain in system resolution 


and the savings in manpower cost. 
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APPENDIX A. NOTES TO THE PROGRAMMER AND THE USER 
This appendix contains a description of the files and libraries used by SEMEX. 
It describes the process of compiling the source codes into an executable file using 
the Microsoft MAKE facility. The corresponding makefile (MAKESEM) is also 
listed. In addition, pointers to programmers and users are provided to help them 


understand and maximize the potential for this program. 


A. SEMEX PROGRAM FILES 

Figure 46 lists the files required to make up the SEMEX executable file. The 
files with the .C extension are the C source files. These are ASCII text files and can 
be edited with any ASCII text editor. The full listings of the source files are found 
in Appendix B. The GLOBAL.H file contains all the prototype definitions required 
by the C language for SEMEX; global variables and constants are also defined here. 
A full listing of GLOBAL.H is also found in Appendix B. 

When compiled, each source file produces a corresponding .OBJ/ object file. 
The object files contain the machine language instructions required to perform the 
commands contained in the source codes. The process of linking resolves any 
external function calls or variables and binds the object codes and the libraries 
together to produce an executable file. After linking, the SEMEX.EXE file is 
created. This process of compiling and linking can be automated by using the MAKE 


facility provided by Microsoft. The MAKE command requires a makefile which 
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Volume in drive F is LEE 
Directory of F:\SEMEX 


GLOBAL 4082 
ACQUIRE 10809 
CETE 9063 
SEMIO 8967 
TAG 13987 
ANALYZE 20055 
SIZE 18975 
SETUP 26789 
SEMEX 10687 
TAG 4809 
SEMIO 3360 
SIZE 9204 
ACQUIRE 4348 
ANALYZE 7965 
SETUP 11023 
CLIP OBJ 2941 
SEMEX OBJ 3700 
SEMEX EXE 127481 
IMGCVT EXE 57405 





Figure 46. Listing of SEMEX C source files, object files and executable 
files using-the DOS DIR command. 


contains the file dependencies and the commands to be executed. The MAKESEM 
makefile used in developing the current version of SEMEX is shown in Figure 47. 
MAKE will determine the validity of the object files based on the dates last 
modified. If the source file has a date later than that of the corresponding object 
file, MAKE will recompile the source file to generate a new object file. If the object 
file is current, no action is taken. In this way, only files that have been modified will 
be compiled. This saves compilation time. Text between the number sign ‘#’ till 
the end of the line are treated as comments and are ignored by MAKE. 
The compile switches currently set are: 


/AL This compiles the source file using the large memory model. Far 
pointers are allocated. By default, maximum optimization is used. 
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# Make file for semex.c and all the dependencies 
# 

# To execute type: make makesem 

# 

# Options for optimizing and no Codeview 

OPT =/Ze 

LOPT =/E /ST:8192 

# Options for no-optimization and Codeview 
#OPT =/Zi /Od 

#LOPT =/CO /ST:8192 /F /PAC 


# Library paths - LARGE ITEXPCplus library and LARGE Microsoft C library 
LIB1 =\pcplus\itex\itexpcml 
LIB2 =\msc\lib\llibce 


acquire.obj: acquire.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 acquire.c 


analyze.obj: analyze.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 analyze.c 


clip.obj: clip.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 clip.c 


semex.obj: semex.c + makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 semex.c 


semio.obj: semio.c 4 makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 semio.c 


setup.obj: setup.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 setup.c 


size.obj: size.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 size.c 


tag.obj: tag.c # makesem 
cl /c /AL $(OPT) /G2 /Fs /SI 100 tag.c 


semex.exe: makesem semex.obj semio.obj size.obj setup.obj tag.obj clip.obj 
acquire.obj analyze.obj 
link /NOD $(LOPT) semex semio size setup tag clip acquire analyze,,, 
Iwin $(LIB1) $(LIB2) 
del *.Ist 


Figure 47. MAKEFILE used in creating SEMEX. 
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/G2 This compiles using the 80286 code instead of the 8088 code. 

WLS This enables extensions to ANSI C and offers additional features beyond 
that provided by ANSI C. One particular feature required by SEMEX 
is the use of casts to produce lvalues (left-hand values). 

/S1 100 This sets the width of the listing file to 100 characters. If a compilation 
error is detected, MAKE will abort leaving a listing file which the user 
can refer to. If no errors are detected after linking, all the listing files 
are deleted. 

The variables $(OPT) and $(LOPT) in MAKESEM are used for macro 
substitutions and will be replaced by their equivalent right-hand expressions specified 
at the beginning of the makefile. Two sets of equates are provided; one is for enable 
CodeView debugging while the other optimizes for speed. The former is currently 
commented out. 

After successful compilation, the linking process will be initiated. As long as 
one object file has been modified, the LINK command will be invoked. This causes 
the object files and the libraries to be linked together. The libraries must be in the 
correct sub-directories as indicated by their paths. If the link is successful, all the 


listing files will be deleted. The libraries used are 


LWIN.LIB Large Memory Model Library for WINDOW BOSS. This 
currently resides in f:\SEMEX. 


ITEXPCML.LIB Large Memory Model Library for ITEX PCplus. 


LLIBCE.LIB Large Memory Model Library for Microsoft C. 
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B. SEMEX OUTPUT FILES 

SEMEX may generate one or more output files. Some of these files are used 
by other modules in SEMEX, while others are files which the user can subsequently 
access. Figure 48 shows the files that can be generated by SEMEX and what they 
represent. SEMEX automatically appends the correct extension. The extension 
provided should not be changed, as this may confuse SEMEX. The filename can be 
entered by simply typing in the characters. The [Ins] key can be used to insert 
characters between existing characters and the [+] and [>] arrows keys can be used 
to position the cursor. 

The image files are all stored in ITEX PCplus compressed file format. This 
generally achieves a storage efficiency of 1.8. However, for highly textured images, 
the compressed file may actually take up more disk space. With the exception of the 
.mat and .met files, the rest of the files are ASCII text files. It should be noted that 
non-ASCII files should not be read by a text editor as the latter may try to format 
the files by inserting special codes. This would destroy the integrity of the files 


render them unusable. 


C. WORKING FROM DIFFERENT DIRECTORIES 

Generally, subdirectories may be used to contain image files from a particular 
burn. In this case, SEMEX should be run from that particular subdirectory from 
which the images are to be processed. To do this, simply type the following 
commands. For example, if the images are contained in g:\SEM\f14, type 

g: [Enter] 


83 


FILE EXT 


SCS 


mg 


unl 


m2 


m3 


.dat 


«his 


.mat 


met 


Figure 48. 


DESCRIPTION AND USAGE 


Sessions file used to record all the activities for the day. The 
filename is based on a 3-letter month code and a 2-digit date 
code (e.g., Mar28.ses). 


Image file generated by ACQUIRE if the save option is 
exercised. The image is generally cropped and complemented. 


Image file generated by CLIP if the save option is exercised. 
The image has been clipped and can be read by the TAG 
module. 


Image file generated by TAG if the save option is exercised. The 
image has been tagged and can be read by the SIZE module. 


Image file generated by SIZE if the save image option is 
exercised. The image should be identical to .im1. This is 
provided merely for test and verification purposes. 


This is the data file generated by SIZE if the save data option is 
exercised. The data consists of the fid, the calculated area 
(AREA C), the measured area (AREA _M), the x-chord and the 
y-chord, for every feature that meets the size specifications (i.e., 
not rejected). 


This is the histogram data file generated by ANALYZE. It 
contains 38 rows of data. Each row represents one bin and 
contains the upper and lower bin limits, the volume in um’, the 
percentage of total volume, the feature count and the percentage 
of total features counted. 


This is the MATLAB data file generated by ANALYZE and 
contains information similar to the histogram data file except that 
it is ina MATLAB format. This is the file which will be asked 
for by SEM.M. 


This is the MATLAB graphics output file which is generated 
when the MATLAB meta command is invoked. The .met files 
can be plotted by typing GRAPH when in DOS. 


List of output files generated by SEMEX. 
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cd \SEM\f14 [Enter] 

Then, to execute SEMEX from f:\SEMEX, type 

f:\SEMEX\SEMEX [Enter] 

To run MATLAB, first make sure that the MATLAB path Is set. For example, 
if MATLAB resides on drive f: in a directory called \MATLAB, then the following 
line must be added to the CONFIG.SYS file 

MATLABPATH =f:\MATLAB 
In order for the path to be set, the IBM system would have to be reset by 
simultaneously pressing [Ctrl], [Alt] and [Del]. After the system has booted up, type 
MATLAB at the DOS prompt. The histogram plotting function is invoked by typing 

sem [Enter] 
when inside MATLAB. The function will request for a histogram filename. Type 
in the filename and the .his extension. It will then ask whether to input the Malvern 
data. If the user chooses to, the function will display each bin limit and request for 
the corresponding Malvern data. After the bins are filled, MATLAB will plot out 
the histogram. The MATLAB meta function can be invoked to save the plot which 


can then be printed out using the GRAPH.BAT command in DOS. 


D. SPECIAL KEYS TO NOTE 
Besides the keys mentioned above, the following keys are important and useful 
to note. 
[Esc] This key is to terminate SEMEX when in the main menu. It can also 


be used to abort any unwanted command. For example, S/ZE usually 
proceeds automatically after TAG finishes (unless disabled during 
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SETUP). However, if the user finds that the tagging Operation has 
produced wrong results (possibly because the size parameters have been 
wrongly set), then he or she could press the [Esc] key when any prompt 
appears. This will terminate TAG and return the user to the main 
menu. 


[F1] This function key is used to obtain help information. The help is 
context-sensitive, where available. 


Mouse support is not available at this time. However, the programmer can 
refer to the WINDOW BOSS manual on how to incorporate functions supporting the 


mouse. 
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APPENDIX B. PROGRAM LISTINGS 
This appendix contains all the source code listings for SEMEX and its modules. 

The listings are documented with comments at the start of each function and at the 
end of each line. These are marked with a slash and an asterisk in the following 
way, /* Comment */. Also included is the MATLAB script file SEM.M. The 
sequence of the listings are as follows: 

GLOBAL.H 

SEMEX.C 

Se PUP C 

ACQUIRE.C 

CHIP -C 

TAG.C 

SIZE.C 

ANALYZE.C 


SEM.M 
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[JARRARARAROEREAARERRERRARA AE ARANA RAR ANI AARERA CARA EAN FANIA NARA ARAN 


* GLOBAL .H Include files for use with SEMEX programs 

He Ve Ae Ye AAA AAA AA ve A AA AAA VA AAA AAA AAA AA AA AV Y A ve ve We We ve Ye ve We ve ve We ve ve ve Wy We ve ye ve ve Wy vr 
j 
#include <math.h> 
finclude "itexpfg.h" /* ITEX PCplus function definitions */ 
#include "stdtyp.h” /* ditto */ 
#include "window.h" /* WINDOW BOSS function definitions */ 


/* Program constants for the SEM Extraction (SEMEX) program 
used in conjunction with the PCVISIONplus Frame grabber and 
w ITEXPCplus library functions */ 


/* PCVISIONplus Board Settings w/ 


#define MEMBASE OxDOO00L /* base memory start address */ 
#define REGBASE 0x300 /* base register start address */ 
#tdefine MEMORY DUAL /* memory type */ 


/* Frame Dimensions */ 


#define XSIZE 312 /* Number of pixels in X direction */ 
#define YSIZE 512 /* Number of pixels in Y direction */ 
#define DEPTH 8 /* Number of bits per pixel */ 


/* AOI (area of interest) settings */ 


#define IXS 0 /* Initial X Starting Point */ 

#define IYS 0 /* Initial Y Starting Point +/ 

define NROW 480 /* Total Number of rows in image */ 
define NCOL 512 /* Total Number of colums in image */ 
#define LASTROW 479 /* Last row in image */ 

#define LASTCOL 311 /* Last columm in image */ 

/* Threshold Limits */ 

#define LOWEST 0 /* Equates to Black for lowcut value */ 
#define HIGHEST 255 /* Equates to White for highcut value */ 
#define BLACK LEVEL 0 /* Indicates a feature */ 

#define WHITE_LEVEL 255 /* Represents background */ 

#define HIGH 254 /* Highest Tag value */ 

#define LOW 1 /* Lowest Tag Value */ 

#define GRAY 128 /* Middle Gray level */ 


/* Miscellaneous Limits AND Defaults */ 
#define MAXFLEN 
#define MAXLEN 
#define MAXCLEN 


/* Define enter key */ 


#define ENTER 


20 


200 


50 


Ox0d 


/* 
/* 
p* 


ye 


Max filename length */ 
Max comment length allowed by ITEX */ 
Comment length */ 


[ENTER] key signature */ 
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15% 


* Global Variables Declarations 


ay 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 


5 


char 
char 
char 
char 
char 
FILE 
floa 
long 
int 
int 
int 
int 
int 
int 
int 
ant 
int 
int 
floa 
int 
int 
int 
int 
int 
int 
int 
int 


filename [MAXFLEN] ; 
session[(MAXFLEN] ; 
comline[MAXLEN] ; 
comlinel([MAXCLEN] ; 
comline2[MAXCLEN]; 
"fp; 

t ASPECT_RATIO, 
TOTAL; 

OVERSIZE; 

UNDERSIZE: 

GAIN_LVL; 

OFFSET LVL; 

LT_MARGIN, 

RT_MARGIN; 

DF_GAIN ; 

DF OFFSET; 

DF_LM; 

DE TRM; 

t VSCALE; 

DF_VSCALE,; 

DF_INVERT; 

LOAD_RAW; 

LOADCLIP: 

DF SIZE: 

LOADTAG; 

DO_SEQ; 

HELP LYL; 


fa 
iB} 
/* 
/* 
Ja 
Je 
ya 
po 
rs 
/* 
{* 
/* 
/* 
/* 
q" 
/* 
Jae 
{/* 
/* 
/* 
qe 
/* 
TS 
/* 
/* 
/* 
J" 


image filename */ 

session filename */ 

entire comment line */ 

comment line 1 */ 

comment line 2 */ 

session file pointer */ 

X to Y aspect ratio of a pixel */ 
Total feature counter */ 


Default too large feature */ 
Default too small feature */ 
Initial Gain level */ 

Initial Offset level */ 

Initial Left Margin */ 

Initial Right Margin */ 

Default gain flag */ 

Default Offset flag */ 

Default left margin flag */ 
Default right margin flag */ 
Initial vertical scale factor */ 
Default vertical scale factor flag */ 
Auto Complement flag */ 


Auto load RAW image flag */ 

Auto load Clipped image flag */ 

Default size limits flag */ 

Ask to load tagged image flag */ 

Flag to sequence through whole process */ 
Help Level */ 


* External function prototypes found in SEMIO.C 


rd 


extern int getim(WINDOWPTR w,int n); 
extern int putim(WINDOWPTR w,int n); 
extern void chgext(char *fd, 


/* End ot file GLOBAL.H */ 


char *fs, 


/* read image function */ 

/* save image function */ 

char *ext); 

/* change extension function */ 
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/* ERARARARAARARAAARNAARARARARARRARRRARARRRARRARARRRRARRRARARARARARARARARRARAARARRRARARRA ARA AAA 


» FILENAME.:. SEMEXRE 
* LAST MODIFIED: 12 Mar 91 by LEE YEAW-LIP 


w = 





=== a me ee ee e m e = 





* PURPOSE : Scanning Electron Microsope Extraction Program (SEMEX) 


X This program uses extensive windowing and calls the 

* following functions: 

* acquire() - acquire images from SEM photographs using 
w vidicon camera and PCVISIONplus framegrabber 
‘ clipmain() ~ clips (threshold) the image 

" tagmain() - tags the features on the image 

4 sizemain() - sizes the features 

x analyze() - analyze the results by building histogram 
A setup() - set up equipment and change the level of 

* user control 

Y Arre hr fe de fe ye ee ee ee ee ee de Ye de de ve we e we ve Ye Ye Ye de de Yr Pe ve fe Pr Fe ve Pe de Pe de de de ve Fe we A Y e 
Y 


/ 


ffinclude "global.h" /* global defines */ 
include <time.h> /* timing prototypes */ 
#include <graph.h> /* graphics prototypes */ 
#include <dos.h> /* dos function prototypes */ 
#include <string.h> /* string handling prototypes */ 
/* 
* Global Variables Defines 
oy 
char filename (MAXFLEN] ; /* image filename */ 
char session(MAXFLEN); /* session filename */ 
char comline[MAXLEN]; /* entire comment line */ 
char comlinel(MAXCLEN] ; /* comment line 1 */ 
char comline2(MAXCLEN]); /* comment line 2 */ 
ELLE" fp /* session file pointer */ 
float ASPECT_RATIO = 1.200; /* X to Y pixel aspect */ 
long TOTAL; /* keeps track of total feature count */ 


/* The following can be redefined during run time by using setup() */ 


int GAIN LVL = 0 : /* Initial Gain - set to highest */ 

int OFFSET_LVL = 60 ; /* Initial Offset - set to midpoint */ 

int DF_GAIN = FALSE; /* Don't use default gain */ 

int DF_OFFSET = FALSE; /* Don't use default Offset */ 

int LT_MARGIN = 0 ; /* Initial Left Margin - leftmost */ 

int RT_MARGIN = 512 ; /* Initial Right Margin */ 

int DF LM = TRUE ; /* Use default left margin = 0 */ 

int DF_RM = FALSE; /* Don’t use default right margin */ 

int OVERSIZE = 100 ; /* Default too large feature */ 

int UNDERSIZE = 1 . /* Default too small feature */ 

int DF_SIZE = TRUE $: /* Use default size limits */ 

float VSCALE = 1.0 ; /* Default Vertical scale factor */ 

int DF_VSCALE = TRUE ; /* Use default scale */ 

int DO SEQ = TRUE ; /* Auto sequence through whole process */ 
int DF_INVERT = TRUE ; /* Complement Automatically */ 

int LOAD_RAW = FALSE; /* Don’t load RAW image automatically */ 
int LOADCLIP = FALSE; /* Don’t load Clipped image automatically */ 
int LOADTAG = FALSE; /* Don’t load Tagged image automatically */ 
int HELP LVL = TRUE ; /* Enable help screens */ 

unsigned blue = BLUE ; /* remap for mono */ 
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main(void) 

{ 
7* Prototype declarations */ 
extern int setup(void); 
extern int acquire(void); 
extern int clipmain(void); 
extern int tagmain(void); 
extern int sizemain(void); 
extern int analyze(void); 
extern clock_t clock(void) ; 
extern void fginit(void); 
extern void session_name(void); 


WINDOWPTR wn; 

ine i: 

int watrib,batrib; 
int rv = 0; 

int rerr 0; 

int row,col; 


clock t start; e 
static struct pmenu smenu = ( 
Op FALSE. 0. 
dae UBS! St 
Doe Sere be nse 0; 
ae SEM Extraction Program“, 
4, 3, “Naval Postgraduate School", 
6, 6, "1. setup SEMEX AN 
7, 6, "2. Acquire Image Mig oe 
8, 6, “3. Clip Image eg es 
9, 6, "4. Tag Features ea 
10, 6, "5. Size Features al 
11, 6, "6. Analyze Features”, 6, 
13, 3, “Use cursor keys to select”, 
14, 3, "Press [Enter] to execute 
ISR, "Press: [Esc] to quit 
99, 99, "",99 ) 


start = clock()}; 
if(wns_mtflg =* 7) blue = BLACK; 


printt(’ Initializing. 
fginit(); 
strcpy(filename," 
session_name(); 


wWatrib = v_setatr(WHITE, BLACK,0, BOLD); /* 


for(i=0; i<25: i++) { 
v_locate(0,i,0); 
v_wca(0, Oxb0, watrib, 80); 
) 
v_hidec(); 
wn_init(); 


pa 
/* 
/* 
/* 
ma 
1" 
EX 
p* 
/* 


/* 
/* 
pe 
/* 
/* 
/* 


Please wait.. 


ads 


setup function */ 

Acquire images from camera */ 
Clip Image */ 

Feature Tagging */ 

Feature Sizing */ 

Analyse Feature data */ 

Returns Number of clock ticks */ 
Frame Grabber Initialization */ 
Use date/time as session name */ 


One Window */ 

scratch integer */ 

scratch atributes */ 

for menu choice */ 

return value from functions */ 
for positioning windows */ 
Time variable */ 


/* define main menu */ 
/* page 0, window open, lndx */ 
/* accept menu field 3 thru 8 */ 
/* field 0 ~ info *f 
0, /* field 1 = info */ 
O /* field 2 = info */ 
fe field 3 - choice 1 */ 
/* field 4 - choice 2 */ 
/* field 5 ~ choice 3 */ 
/* field 6 — choice 4 */ 
Iž field 7 ~ choice 5 %/ 
/* field 8 ~ choice 6 */ 
O, /* field 9 - info */ 
sO. J7* field 10 = into, */ 
O/A feld II c mni T 
/* menu terminator */ 
/* start timing */ 
/* remap if mono */ 
eos 
/* set up Frame Grabber */ 
/* blank out filename */ 
/* use date/time for session */ 


window attribute */ 
build the back drop */ 
position cursor */ 

the fast way */ 


p* 
[s 
p* 


hide the cursor */ 
save entry screen */ 


/* 
/* 
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/* 
* Popup Menu 


xy 
do { 
watrib = WHITE<<4|BLACK; /* window colors */ 
batrib = blue<<4|WHITE; /* border colors */ 
if( (!DO SEQ) || (rv < 3) || rerr ) /* popup menu to get user choice */ 
rv = wn_popup(0,3,23,31,17,watrib,batrib,&smenu, FALSE) ; 
else { /* auto sequence */ 
Eve /* get next menu list */ 
if Cry >. 5)srv = 0; /* end of menu list */ 
} 
Crystal /* set to internal sync */ 
switch (rv) { /* do user command */ 
case 1; /* setup configuration */ 
rerr = setup(); 
break; 
case 2: 
rerr = acquire(); /* Acquire Image */ 
break; 
case 3: 
rerr = clipmain(); /* Clip Image */ 
break; 
case 4; 
rerr = tagmain(); /* Feature Tagging */ 
break; 
case 5: 
rerr = sizemain(); /* Feature Sizing */ 
break; 
case 6: /* Analyse Feature data */ 
rerr = analyze(); 
break; 
case 99: /* exror or ESC key */ 
default: 
rerr =; 
break; 


} /* end switch */ 
} while(rv !=99); 
wn_exit(); /* restore entry screen */ 
_clearscreení GCLEARSCREEN); /* clear screen */ 
printf("ininininitSession's Activities have been recorded in Zs", session); 
printf("\n\n\tSEMEX was on for Z.1f minutes.”, 

(float) (clock()=start) / (float) (CLK ICK *i(ciock t) 6560), ); 
printf("\n\n\tDeveloped by ECE Dept, Naval Postgraduate School"); 
printf("\n\n\tHave A nice day! \n\n\n"); 
if(( fp = fopení(session,"a")) == NULL) 

printf("a\n\n\Unable to open session file ZsinitSEMEX Aborted.", session); 
else { 
fprintf(fp, "inSEMEX was on for Z.1f minutes.", 
(float) (clock()-start) / (float) (CLK_TCK * (clock_t) 60) ); 
fclose(fp); 
y) /* if-else */ 
exit(0); /* Successful termination */ 
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p* 
* Initial PCVISIONplus Frame Grabber setup (Refer to global.h for defines) 


ee 

void 

fginit(void) 

{ 
sethdw(REGBASE , MEMBASE , MEMORY) ; /* Set hardware definitions */ 
setdim(XSIZE,YSIZE. DEPTH); /* set frame dimensions */ 
feon(); /* turn on frame grabber */ 
initialize(); /* set up registers and LUTs */ 
setcamera(0); /* Select Camera */ 
extsync(); /* Use Camera video sync */ 
waitvb(); /* syne to vertical blanking */ 
select_mem(MEM A); [* Select frame A */ 
display_mem(MEM A); /* display frame A */ 
sclear (GRAY); /* clear whole TV screen */ 
setlut (INLUT,0O): /* Use INPUT LUT, bank 0 */ 
setlut(GRNLUT,O); /* Use GREEN OUTPUT LUT, bank O */ 

) 


/* session _name() gets the current date and time and creates a session 
* filename 


Ry 
void 
session_name(void) 
{ 
struct dosdate t date; /* dos.h date structure */ 
struct dostime_t time; /* dos.h time structure */ 
static char *month[12] = { "Jan, “Feb”, “Mar”, "Apr", “May”, “Jun” 
“Juls AUR e Sepi, Oct. ANO” on Dec >: 
_dos_getdate (&date); /* get current date */ 
_dos_ gettime (&time); /*™ get. current time */ 
if (date.day < 10) /* single digit date, pad zero */ 
sprintf(session,"2Z3s0Z1ld.ses " month(date.month - 1), date.day); 
else /* double digit date */ 
sprintf(session, "%3sZ2d.ses " month[date.month - 1], date.day); 


if ( (fp = fopen(session,"a")) == NULL) { 
printf("a\n\n\Unable to open session file 2%s\nSEMEX Aborted.”, session); 
exit; /* terminate */ 
} 
else { 
if (time.minute < 10) /* single digit minute, pad zero */ 
fprint£(fp, "inOpening Session on 22d 23s 24d at 22d:02Z1d.", 
date.day, month[date.month - 1], date.year, time.hour, time.minute); 
else /* double digit minute */ 
fprintf(fp,”\nOpening Session on 22d 23s 24d at 22d:22d.", 
date.day, month[date.month - 1], date.year, time.hour, time.minute); 
fprint EC Ep == == A O 
fclose(fp); /* close session file */ 


/* End ot file SEMEX.C */ 
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FILENAME: SETUP.C 
CALLED BY: semex() and calls clip() 
LAST MODIFIED: 12 Mar 91 by LEE YEAW-LIP 
PURPOSE: Setup is used to set up the SEMEX program to configure it 
to suit the user’s needs. Various prompts can be turned on 
or off to stream-line the processing of SEMs. 
It changes the set of global variables and flags which the 
various functions in SEMEX make use of, or test for. 
Latest modification incorporates check_equipment which 
facilitates the setup of the light table and camera. 
measure_line() has also been added to determine the 


magnification of the image. 
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r 
#include "global.h" 
#include <string.h> 


setup() 

{ 
extern void bool2YN(int, char *); /* prototype declaration */ 
extern void check equipment(void) ; [* Gittow, 
WINDOWPTR wn; /* one window pointer */ 
WIFORM frm; /* form pointer “/ 
int watrib, batrib; /* window, border and */ 
unsigned fatrib; /* field attributes */ 
Static char *emsgl = "Only range of 0 to 100 allowed. Press any key"; 


Static char *emsg2 = "Only range 2 to 256 allowed. Press any key"; 
Static char *emsg3 = "Only range 1 to 255 allowed. Press any key”; 
Static char *emsg4 = "Only range 0 to 255 allowed. Press any key”; 

static char *emsg5 = "Only range 256 to 512 allowed. Press any key”; 


char DUALES /* string buffers */ 
char bZ(5)]. 53131. [* ditto "7 
char b4[5), bS5[3); S /* ditto */ 
char b6[5), b7[3); Ia ditto */ 
char b8[9], b9[3]; /* ditto */ 
char b10[5], b11[5], b12[3); /* ditto */ 
char b13[6), b14[3]: /* unused at the moment */ 
char b15[3), b16[3], b17[3], b18[3); /* Y/N input buffers */ 
char b19[3), b20[3); /* ditto */ 
Check_Equipment(); /* check camera and lights */ 
fatrib = (BLUE<<4) | WHITE | BOLD; /* field color */ 
watrib = v_setatr(WHITE, BLUE,0,0); /* window color */ 
batrib = v_setatr(RED,WHITE,0, BOLD); /* border color */ 
wn = wn_open(500,1,5,60,19,watrib,batrib); 
if (twn) { 

printf("\a\n Unable to open window. Aborting..." ); 

exit(1); 
) 
wn_title(wn," SETUP DEFAULTS "); 
frm = wn_frmopn(22); /* open 21 + 1 fields */ 
if (!frm) { 

printf("\a\n Unable to open form. Aborting..."); 

exit(1); 
) 
itoa(GAIN_LVL,b0,10); 
wn_gint(SET,frm, 0,wn, 2, 1,”°GAIN LEVEL 5 "GLatrib,*_”, 


&GAIN_LVL,3,0,100,b0,emsg1,emsgl); 
bool2YN(DF_GAIN,b1); 
wn_gbool(SET,frm, 1,wn, 2,30,"Use Default[Y,N] : ",fatrib,’ ', 
&DF_GAIN,b1,NSTR,NSTR); a 
itoa(OFFSET_LVL,b2,10); 
wn 'gint(SEI, frm, 2,wn, 3, 2, OFFSET LEVEL o tatribo os, 
SOFFSET_LVL,3,0,100,b2,emsgl,emsg1); 
bool2YN(DF_OFFSET ,b3) ; 
wn_gbool(SET,frm, 3,wn, 3,30,"Use Default[Y,N] : ",fatrib,’_’, 
£DF_OFFSET,b3,NSTR,NSTR); 
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) 


/* 


itoa(LT_MARGIN,b4,10); 
wn_gint(SET,frm, 4,wn, 4, 1,"LEFT MARGIN Atari, 
£LT_MARGIN,3,0,511,b4,emsg4,emsg4); 
bool 2YN(DF_1M,b5); 
wn_gbool(SET,frm, 5,wn, 4,30,"Use Default[Y,N) : ",fatrib,’_’, 
&DF_LM,b5S,NSTR,NSTR) ; 
itoa(RT_MARGIN,b6,10); 
wn _gint(SET,frm, 6,wn, 5, 1,”°RIGHT MARGIN E ENERO 5 
&RT_MARGIN,3,1,512,b6, emsg5,emsg5) ; 
bool2YN(DF_RM,b7); 
oone POOLSEL frm, 7/,wn, 3,30, Use Default(Y,N]. 3 " fatrib,” ”, 
£DF_RM,b7,NSTR,NSTR); 
SpEmob(De. 27. .3£”,VSCALE); 
An RELOat(SET,frm,8,wn, 6, 1, Y-SCALE FACTOR ”,fatrib,*_', 
£VSCALE,7,3,0.0,999.0,b8,emsg81l,emsgl); 
bool2YN(DF_VSCALE,b9); 
Wo epool (SEI, frm, 9,n, 6,30,"Use Default [Y,N] : “,Eatrib,*_”, 
£DF_VSCALE,b9,NSTR,NSTR); 
itoa(OVERSIZE,b10,10); 
wn gint (SEIT, frm,10,wm, 7, 1, "Max Feature Size 2 ”,fátrib,*_”, 
&OVERSIZE,3,2,256,b10,emsg2,emsg2); 
itoa(UNDERSIZE,b11,10); 
wn gint(SET,frm,11,wn, 8, 1,"Min Feature Size : %”,fatrib,’_’, 
&UNDERSIZE,3,1,255,b11,emsg3,emsg3); 
bool2YN(DF_SIZE,b12); 
wn_gbool(SET, frm,12,wn, 8,30,"Use Defaults[Y,N]): ",fatrib,’_’, 
&DF_SIZE,b12,NSTR,NSTR); 
bool2YN(HELP_LVL,b15); 
wn_gbool(SET,frm,13,wn,10, 1,"ALL: Enable HELP screens [Y,N] 
fatrib,’—’ -GHELP LVL, b15,NSIR,NSIR); 
bool2YN(DO_SEQ,b16); 
wn_gbool (SET, frm,14,wn,11, 1, “SEMEX: CLIP, TAG and SIZE without asking[Y,N]): 


fatrib,' ’,&DO SEQ,b16,NSTR,NSTR); 
bool2YN(DF_INVERT,b17); 
wn_gbool(SET, frm,15,wn,12, 1, "ACQUIRE: Complement Image without asking [Y,N): 


fatrib,* *,£DF_INVERT,b17,NSTR,NSTR); 

bool2YN(LOAD RAW,b18); 

wn_gbool(SET, frm, 16,wn,13, 1,"CLIP: Load RAW Image without asking[Y,N) 
fatrib,’ ’,&LOAD RAW,b18,NSTR,NSTR); 

bool2YN(LOADCLIP,b19); 

wn_gbool (SET, frm,17,wn,14, 1, TAG: Load CLIPPED Image without asking[Y,N] 
fatrib,’ ’,&LOADCLIP,b19,NSTR,NSTR) ; 

bool2YN(LOADTAG,b20); 

wn_gbool (SET, frm,18,wn,15, ESI ZE: Load TAGGED Image without asking[Y,N] 


fatrib,’_’ ,&LOADTAG,b20,NSTR,NSTR); 

wn_gtext(SET,frm,19,wn,16, 1,"Session Filename: ”,fatrib,’_',18,session,NSTR,NSTR); 
wn_dtext(SET,frm,20,wn,18,1,"Press [Esc} to accept Existing defaults"); 
if(!wn_frmget(frm)) { /* read form */ 

printf("\a\nMemory corrupted. Aborting wn_frmget()"); 

exit (1); 
} 
wn_frmcls(frm); /* close form */ 
wn_close(wn) ; /* close window */ 
return(0); 


* Convert 1’s to Y’s and 0’s to N’s 


= 
void 


bool2YN(int n, char *s) 


{ 


if (n == TRUE) 
strcepy(s,"Y”); 

if (n == FALSE) 
strcpy(s,”N"); 
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FUNCTION NAME: check_equipment() 

PURPOSE : This function allows setting up of the light table and 
camera prior to actual acquisition of images. 
Afer acquiring a test image, it can be inverted and 
clipped to ascertain the uniformity of illumination. 
A graphic cursor facility is also provided to measure 


various features. 
ve ve ve ve ve Fe ve Fe Fe ve Fe Pe Fe We Se ve ve Fe Se Xe Xe ve te He de WY Ve se Pe WY We Ve Ae Pr We oe Wh oe oe Fe WY oe oe ove oe oe ve He ve Fe He te se oe te ve te ve ve ve oe oe oe ve Wr Oe 


/ 


void check_equipment(void) 


{ 


extern int measure line(WINDOWPTR w); /* function prototype */ 
extern int clip(WINDOWPTR w, int t); /* ditto */ 

WINDOWPTR wn; /* window handle */ 

int watrib, batrib; /* Window, border and */ 
unsigned fatrib; /* field attributes */ 

char c; /* scratch for user response */ 
int i: /* scratch index */ 

int done; /* scratch flag */ 

int gain, offset; /* gain and offset setting */ 
watrib = v_setatr(WHITE,BLUE,0,0); /* set window attribute */ 
batrib = v_setatr(RED,WHITE,0,BOLD); /* set border attribute */ 


wn = wn_open(500,8,13,50,14,watrib,batrib); /* open dialog window */ 
if (!wn) { 

printf("\a\n\n Unable to open window"); 

exit(l): 
} 
wn_title(wn,” EQUIPMENT SETUP "); 
wn_printf(wn,"\n\n\tSet up camera and lights[Y} ?"); 
v kflush(); /* empty keyboard buffer */ 
c = getch(); /* get user response */ 
while ( (c != 'n') €£ (c != 'N’) ){ /* repeat until No */ 

/* Turn on camera */ 

wn_clr(wn); /* clear window */ 

n_title(wn," SET UP CAMERA AND LIGHTS ”); 

n_puts(w,2,5,"Turning on GRAB Mode...”); 


setcamera(0); /* connect camera 0 */ 
extsync(); /* external sync */ 
grab(NO WAIT); /* ensure grab mode ON */ 


grab(NO WAIT); 

/* Fine tuning the board’s gain and offset */ 

wn_puts(wn,1,5," GAIN AND OFFSET SETTING ”); 

gain = GAIN _LVL; 

setgain(gain); /* initial gain setting */ 
v_kflush(); /* prevent spurious input */ 
wn_puts(wn,7,1,"Use [+] and [-] keys to adjust"); 
n_puts(wn,8,1,"Press [ENTER] to continue”); 
wn_puts(wn,2,1,"Gain (0 highest, 100 lowest): ”); 
n_print£(wn,"23d"”, gain); 

while ( (c=getch()) != ENTER ) 


{ 
wn_locate(wn,2,31); 
if ( (c == '+') ££ (gain < 100) ) 
gain = gain + 5; /* increment in steps of 5 */ 
if ( (c == >") && (gain > 0) ) 
gain = gain — 5; /* decrement in steps of 5 */ 
setgain(gain); /* change the camera gain */ 
wn_printf(wn,"23d",gain); 
} 
GAIN_LVL = gain; /* Update new gain level */ 
DF_GAIN = TRUE; /* Use default gain */ 
offset = OFFSET_LVL; 
setoffset(offset); 


wn_puts(wn,4,1,”"Offset (0 darkest, 100 lightest): "); 
wn_printf(wn,"Z3d" offset); 

while ( (c=getch()) != ENTER) 

{ 
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wn_locate(wn,4,35); 


if ( (c == '+") && (offset < 100) ) 
offset = offset + 5; /* increment in steps of 5 */ 
if ( (c == '-') ES (offset > 0) ) 
offset = offset - 5; /* decrement in steps of 5 */ 
setoffset(offset); /* change the camera offset */ 
n_print£íwn,"23d" offset); 
) 
OFFSET LVL = offset: /* Update new offset */ 
DETOPFSET = TRUE, /* Use default offset */ 


/* Snap a frame (this will stop further acquisition) */ 
wn_title(wn,” ACQUIRE IMAGE - SNAP MODE ”); 


done = FALSE; /* check whether image snapped */ 
wn_clr(wn), 
Vekt lush ():; /* prevent spurious input */ 


wn_printf(wn,”\n\n\tPress [SPACEBAR] to snap an image”); 
wn_printf(wn,”"\n\t repeat until satisfied"); 
while ( (c = getch()) != ENTER) 


{ 
wn_puts(wn,6,4,"Wait..."); 
waitvb(); /* Wait for vertical blanking */ 
snap(WAIT), /* acquire a frame */ 
wn_puts(wn,6,4,"Done! "); 
wn_puts(wn,12,4,"When done, press [ENTER]"); 
done = TRUE; /* flag that image taken */ 
} 
if (!done) { 
waitvb(); /* Wait for vertical blanking */ 
snap(WAIT); /* turn grab off */ 
} 
crystal(); /* internal sync for stability */ 
setcamera(1); /* disconnect camera 0 */ 
c = ?’Y’; 
while ((c == 'y') |] (c == 'Y') ) { /* repeat */ 
wn_clr(wn); /* clear screen */ 
if (measure_line(wn) ) /* measure line lengths */ 
return; /* error detected. Abort */ 
wn_printf(wn,"\n\n\tMeasure another feature [N]?"); 
v_kflush(); /* empty keyboard buffer */ 
c = getch(); 
} 


wn_clr(wn); 
/* Complement Image in frame memory */ 


v_kflush(); /* empty keyboard buffer */ 
wn_printf(wn,"\n\n\tComplement Image [Y]?"); 
c = getch(); 
if ( (c!='N’) Es (c!='n’) ) { 
DF_INVERT = TRUE; /* Set auto-invert flag */ 
wn_printf(wn,"\n\tComplementing Image..."); 
complement (IXS, 1YS,NCOL, NROW) ; 
) 
else 
DF_INVERT = FALSE; /* Reset auto-invert flag */ 
/* Threshold image */ 
wn_printf(wn,"\n\n\tUpdating frame memory..."); 
maplut(GRNLUT,0,IXS,IYS,NCOL,NROW); /* update frame memory */ 
wn_clr(wn); /* clear window */ 


threshold(GRNLUT,0,HIGHEST, HIGH) ; /* Initial threshold */ 
wn_printf(wn,"\n\tReduce to guage Lighting Uniformity”); 


clip(wn, HIGH); /* threshold image */ 
wn_clr(wn); /* clear window */ 
wn_printf(wn,"\n\n\tAdjust Lighting and Try again[Y) ?"); 
v_kflush(); /* clear keyboard buffer */ 
c = getch(); /* get user response */ 
linlut(GRNLUT,O); /* restore LUT */ 

) /* end while */ 

wn_close(wn) ; /* ell done "/ 

return; JA return *f 
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* FUNCTION NAME: measure_line() 
* CALLED BY: check_equipment() 
* LAST MODIFIED: 7 Mar 91 by LEE YEAW-LIP 


we ree ee ee ee ee ee re ee ee ee a eS Se 


PURPOSE: 


+ 


a. 
= 


This provides a set of functions f 


Presently step size is 10 pixels 
It also calculates the magnificati 
AUNAR ve ve ve vr ve ve ve ve ve ve ve ve ve oe ve ve ve ve ve ve oe ve ve ir ve ve ve air vy vr ve We vr air ty vy vr Ve oy Ve ve k 
nf 

/* #include "global.h" */ 


> 


We ve de ve ty We e A e e e A or 


or generating a graphics 
und with the directional 
zontal lines. The 


pixel at a time 
ne pixel at a time 


ne pixel at a time 
one pixel at a time 
n steps 

in steps 

steps 


down in steps 


* cursor (cross—hair), moving it aro 
* keys and drawing vertical and hori 
w directional keys function as follows: 
w Up Moves cross-hair up one 
N Down Moves cross-hair down o 
k Left Moves cross-hair left o 
* Right Moves cross-hair right 
s Home Moves cross~hair left i 
= End Moves cross-hair right 
5 PgUp Moves cross-hair up in 
* PgDn Moves cross-—hair 

Y 


on of a 5 micron line. 
¥e ve ve ve ve ve ve ve ve ve ve ve ve ve ve ve ve He ve ve ve vr vy vr 


#tundef UARROW /* undefine to prevent conflict */ 

#undef DARROW JA ditto */ 

#tundef LARROW /* ditto */ 

undef RARROW /* ditto */ 

#undef PAGEUP /* ditto */ 

#tundef PAGEDN [" aitco ~/ 

#undef HOME /* ditto */ 

#undef END /* ditto */ 

#define UARROW 0x48 /* redefine key constants */ 

#define DARROW 0x50 /* ditto */ 

#define LARROW 0x4b Edit ton 

#define RARROW Ox4d fer aitto *7 

#define PAGEUP 0x49 /* ditto */ 

#define PAGEDN 0x51 /* ditto */ 

#define HOME 0x47 /* ditto */ 

#define END Ox4f 

#define STEP 10 /* step size for fast move */ 

static int xpixval(9], ypixval[9]; /* pixal value under cursor */ 

static int line(NROW]; /* allocate memory for a line */ 

/* 

* measure a line marked by two graphic cursors 

x 

measure_line(WINDOWPTR w) 

{ 
extern void put_cursor(int x,int y); /* function prototype */ 
extern void unput_cursor(int x,int y); [* ditto */ 
extern void put_line(int xl,int yl,int x2,int y2); /* ditto */ 
extern void unput_line(int x1,int y1,int x2,int y2); /* ditto */ 
extern void chkkey(char c, int *x, int *y); /* ditto */ 
extern float calc_line(int xl,int yl,int x2,int y2); /* ditto */ 
char Cc; /* scratch for key pressed */ 
int x = 470, y = 65; /* initial cursor location */ 
int x1, y /* store 1st point on line */ 
float len; /* length of line */ 
int color: f* pixel color */ 
v_kflush(); /* prevent spurious inputs */ 


wn_printf(w,"\n\n\tLine Measure Feature"); 


wn_printf(w,"\n\n\tPosition cursor at first point"); 


wn_printf(w,"\n\tusing ARROW keys. Press ( 
put_cursor(x,y); e 
wn_printf(w,"\r\tCoordinates: X1: %3d 

while ( (c = getch()) == ’\0’) [S 


ENTER) when done\n"); 

put cursor at the center */ 
AI 9 x.y): 

if cursor direction key */ 


98 


c = getch(); /* check direction */ 
unput_cursor(x,y); /* remove old cursor */ 
chkkey(c, &, &y); 
put_cursor(x,y); /* put cursor back on */ 
wn_printf(w,"\r\tCoordinates: X1: 26a TIAS y X y): 

) /* end while */ 

xl = xX; /* save Ist point "/ 

yi y; 

put_line(x1,y1,x,y); {/* put out a line */ 


/* get second point on the line */ 
wn_printf(w,"\n\n\tStretch line to second point"); 
wn_printf(w,"\n\tusing ARROW keys. Press [ENTER] when donen"); 


while ( (c = getch()) == '10') /* if cursor direction key */ 

{ 
c = getch(); /* check direction */ 
ubput line(x1,y1,x.y); /* remove old line */ 
unput_cursor(x,y); /* remove. cursor */ 
chkkey(c, &x, &y); /* get new cursor position */ 
put_cursor(x,y); {/* put cursor back on */ 
put_line(x1,y1,x,y); /* put new line */ 
wn_printf(w,"\r\tCoordinates: A2: 238 LL ADO Y) 

) /* end while */ 

len = calc line(xl,yl,Xx.y); /* calculate distance */ 

wn_printf(w,"\n\n\tFeature is 2.1f pixels long.", len); 

if (len > 0.0) ( /* compare with 5 micron ref */ 
VSCALE = len/5.0; /* to give Pixel Conversion Factor */ 


wn_printf(w,"\n\tFor 5 micron vertical line,"); 
wn_printf(w,"\n\tVertical Scale factor is 2.4f pixels/micron”,VSCALE); 
if (( fp = fopen(session,"a") ) == NULL) ( 
wn_printf(w,"\a\n\n\tCANNOT open session file 2Zs", session); 
wn_printf(w,"\n\n\tPress any key to continue"); 


getch(); 

return(1); /* premature termination */ 
} 
else { 


fprintf(fp, "\nSETUP:\tMeasured length is 2.1f pixels long",len); 
fprintf(fp,"\n\tVertical scale factor is 2.4f pixels/micron",VSCALE) ; 
fclose(fp); /* close session file */ 

) /* end if-else */ 

wn_printf(w,"\n\n\tPress Any key when done "); 


v_kflush(); /* empty keyboard buffer */ 
getch(); 
) 
unput_line(x1,y1,x,y); /* remove line */ 
unput_cursor(x,y); /* remove cursor */ 
return(0); /* everything OK */ 
) 
/* 


* Put out a cross-hair graphic cursor at location x,y, saving the 

* pixels under the graphic cursor. The intensity of the cursor 

* will be peak-white or peak-black depending on the pixel intensity 
* if the center point. 


s 

void 

put _cursor(int x,int y) 

{ 
int i: /* scratch counter */ 
inex yl; /* cursor pixel counters */ 
int color; /* pixel color */ 
ce ox 4: /* get the left position */ 
yl= y~ 4; /* get the top position */ 


/* save pixels under the cursor */ 

for ( iè 0; 1 < 9; i++) { 
xpixval[i] = rpixel(xlti,y); /* save pixel value */ 
ypixval({i] = rpixel(x,ylti); 


J39 


if (rpixe 
color 
else 
color 
/* put ou 
for(i 


wp 
“Pp 


) 


/* 


l(x,y) > GRAY) 
= BLACK LEVEL; 


= WHITE LEVEL; 
t the cursor */ 
0; i < 9; itt) { 
ixel(xl+i,y,color); 
ixel(x, ylti, color); 


/* enhance visibility */ 
/* of cursor */ 


draw horizontal */ 
draw vertical */ 


/* 
/* 


* Restore image to original state without cursor 


a 
void 
unput_cursor( 


{ 


int i> 
int: XL oy 
xl = x = 
yi ey = 
for (i= 
wpixe 
wpixe 
) 
} 
/* 
* check key 
ny 
void 


chkkey(char c 
{ 


switch (c 
case 


case 


case 


case 


case 


case 


case 


case 


defau 


int x, int.» 


a 

4; 

0; i < 9; i++ ) { 
1(x1+i,y xpixval(i]); 
l(x,yl+i,ypixval[i}); 


pressed for direction 


Ant A Int ty) 
Al 
UARROW: 
Cys: 
To Cty) 
break; 
DARROW: 
e ; 
LL) 
break; 
LARROW: 
CXE: 
TTC x) 
break; 
RARROW: 
Cex) I; 
if ((*x) 
break; 
PAGEUP: 
(*y) == STEP; 


< 


IYS) 


if ((*y) < IYS) (*y) = LASTROW; 


break; 
PAGEDN: 
(*y) += STEP; 


if ((*y) > LASTROW) (*y) = IYS; 


break; 
HOME : 
(*x) == STEP: 


if ((*x) < IXS) (*x) = LASTROW; 


break; 

END: 

(*x) += STEP; 

if ((*Xx) > LASTCOL) (*x) 
break; 

It: 

break; 


(*y) = LASTROW; 


> LASTROW) (*y) STIIS; 


IXS) (*x) = LASTROW; 


LASTCOL) (*x) = IXS; 


IXS; 


scratch counter */ 
cursor pixel counters */ 


/* 
p* 


/* 
/* 


get left position */ 
get top position */ 


restore horizontal stroke */ 
restore vertical stroke */ 


/* 
/* 


/* move up */ 


[£ round */ 


wrap 


/* down */ 


move 


/* wrap round */ 


as left */ 


move 


/* round */ 


wrap 


/* move right */ 


/* round */ 


wrap 


/* move up fast */ 


/* round */ 


wrap 


/* down fast */ 


move 


5 round */ 


wrap 


/* move left fast */ 


/* round */ 


wr ap 


pe right fast */ 


move 


/* round */ 


wrap 


100 


) /* end switch */ 


) 
/* 
* calculate the distance between two points 
wy 
float 
calce linetint xl, int yl, int x2,int y2) 
{ 
if ( abs(x2 - x1) > abs(y2 - yl) ) /* horizontal line */ 
return( abs(x2-xl) ); 
else /* vertical line */ 
returní( abs(y2-y1) ); 
} 
/* 


* draw a vertical or horizontal line 
* saving the image pixels underneath 


xy 
void 
puts tine(int xi, int yl, int x2, int y2) 
{ 
int i; /* scratch counter */ 
int color; /* line color */ 
if (rpixel(xl,yl) > GRAY) /* enhance visibility */ 
color = BLACK_LEVEL; /* of line */ 
else 
color = WHITE LEVEL; 
if ( abs(x2 - x1) > abs(y2 - yl) ) { /* draw horizontal line */ 
if (x2 > x1) { 
for (i = x1; i <= x2; itt) { /* draw left to right */ 
line[i-x1] = rpixel(i,yl); /* save pixel values */ 
wpixel(i,yl,color); /* before overwriting */ 
} 
} 
else { /* (x2 < x1) */ 
for (i = x2; i <= xl; itt) { /* draw right to left */ 
line[i-x2] = rpixel(i,yl); /* save pixel values */ 
wpixel(i,yl,color); /* before overwriting */ 
} 
) /* end if(x2 > x1)-else */ 
) 
else { /* draw vertical line */ 
if (y2 > yl) { 
for (i= yl; i <= y2; itt) { /* draw top down */ 
line[i-y1] = rpixel(x1,i); /* save pixel values */ 
wpixel(x1,i,color); /* before overwriting */ 
) 
} 
else { Lye = 1) 7 
for (i = y2; i <® yl; i++) { /* draw bottom up */ 
line[i-y2] = rpixel(x1,i); /* save pixel values */ 
wpixel(x1,i,color); /* before overwriting */ 
} 


} /* end if(Y2 > yl)-else */ 
) /* end if(abs(x2 -x1))-else */ 
) 


/* 
* erase a vertical or horizontal line 
* restoring the image pixels underneath 


ey 
void 
unput_line(int xl, int yl, int x2, int y2) 
{ 
int i: /* scratch counter */ 


if ( abs(x2 - x1) > absiy2 - yl) ) ( /* erase horizontal line */ 
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if (x2 > x1) ( 
for (1 = xl: 1 <= X2: 445) 
wpixel(i,y1,linefi-x1]); 
} 
else { 
for (i = x2; i <= xl; i++) 
wpixel(i,yl,line{i-x2]); 
) /* end if(x2>x1)-else */ 
) 
else ( 
if (y2 > yl) i 
for (i = yl: i <= y2; mitt) 
wpixel(x1,i,lineli-y1]); 
) 
else ( 
for (11 =y2: 1 <= IRIT) 
wpixel(x1,i,lineli-y2])); 
) /* end if(y2>y1)-else */ 
) /* end if(abs(x2-x1))-else */ 


f* End of tilessrlur.€ */ 


/* 
p* 


p* 
ES 
p* 


/* 


/* 
px 


[* 
ix 
/* 


erase left to right */ 
by restoring original values 


(x2 < xl) */ 


erase right to left */ 
by restoring original values 


erase vertical line */ 


erase top down */ 
by restoring original values 


(yl Y LJ 


erase bottom up */ 
by restoring original values 
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FILENAME > ACQUIRE.C 
CALLED BY: SEMEX main program 
LAST MODIFIED: 17 Mar 91 by LEE YEAW-LIP 
PURPOSE : This function allows images to be acquired from the 
PCVISIONplus framegrabber board. Afer acquiring a frame, 
the image can be cropped and inverted if desired before 
it is saved. Values applied to the image are saved in 
a session file as well as in the comment line of the image 


header. 
te te te te te oe He te te He te te Xe ete Xe He Xe te te Xe He He He He He He He He Xe He He Xe He Xe Xe Xe ¥e He ¥e Xe He Fe He He He Ye He He He Ye ec He Xe He Xe Fe Xe He Ye Fe Fe Xe Xe Xe Xe Ve te 





Se te ee SS ee ee ee ee a. 








#include "global.h” /* Necessary defines */ 
#include <string.h> /* string prototypes */ 


/ 


{ 


* 
Y 


* 


Acquire image from camera, invert, crop and then save it 


/ 

acquire() 
WINDOWPTR wn; /* window handle */ 
int watrib, batrib, /* Window, border and */ 
unsigned fatrib; /* field attributes */ 
char €; /* scratch for user response */ 
char s>[5]: /* scratch string */ 
int i; /* scratch index */ 
int nocam; /* flag indicating camera off */ 
int Em, rm; /* left and right margin setting */ 
int gain, offset; /* gain and offset setting */ 
watrib = v_setatr(WHITE,BLUE,0,0); /* set window attribute */ 
batrib = v_setatr(RED,WHITE,0,BOLD); /* set border attribute */ 


wn = wn_open(500,8,13,50,14,watrib,batrib); /* open dialog window */ 
if (!wn) { 
print£("lalnin Unable to open window”); 
exit(1); 
) 
pe 
* Read in an image or turn on camera 
= 
wn_title(wn,' IMAGE ACQUISITION "); 
n_print£(wn,"ininitRead Image from file [N}?"); 


v k£flusht); /* empty keyboard buffer first */ 
c = getch(); /* get user response */ 
if ( (c == 'Y’) || (c == 'y') ) ( 
getim(wn,0); /* read image from file */ 
nocam = TRUE; /* Note camera is off */ 
} 
else { 
wn puts(wn,2,3,"Turning on GRAB Mode... ”); 
setcamera(0); /* connect camera 0 */ 
extsync(); /* external sync */ 
grab(NO WAIT); /* ensure grab mode ON */ 
grab(NO WAIT); 
nocam = FALSE; /* camera on flag */ 
ae 
* Fine tuning the board's gain and offset. 
i 


wn_title(wn," GAIN AND OFFSET SETTING "); 
gain = GAIN LVL, 
setgain(gain); /* initial gain setting */ 
wn_puts(wn,7,1,"Use (+) and [-] keys to adjust"); 
wn_puts(wn,8,1,"Press [Enter] to continue”); 
if (!DF_GAIN) { 
wn_puts(wn,2,1,"Gain (0 highest, 100 lowest): ”); 
wn_printf(wn,"23d",gain); 
while ( (c=getch()) != ENTER ) 
{ 


103 


wn_locate(wn,2,31); 


if ( (c == ’+') && (gain < 100) ) 
gain = gain + 5; /* increment in steps of 5 */ 
if ( (c == ’=") && (gain > 0) ) 
gain = gain - 5; /* decrement in steps of 5 */ 
setgain(gain); /* change the camera gain */ 
wn_printf(wn,"43d", gain), 
} 
GAIN LVL = gain; /* update global gain value */ 
} 
offset = OFFSET LVL; 
setoffset(offset); 


if (1DF OFFSET) { 
wn_puts(wn,4,1,"Offset (0 darkest, 100 lightest): "); 
wn_printf(wn,"23d" ,offset) ; /* display current offset */ 
while ( (c=getch()) != ENTER) { 
wn_locate(wn,4,35); 
if ( (c == '+") && (offset < 100) ) 


offset = offset + 5; /* increment in steps of 5 */ 
if ( (c == ’~’) && (offset > 0) ) 
offset = offset - 5; /* decrement in steps of 5 */ 
setoffset(offset); /* change the camera offset */ 
wn_printf(wn,"%3d",offset); /* display new offset */ 
) 
OFFSEI_LVL = offset; /* update global offset value */ 
} 
} /* end if-else */ 
/* 
* Snap a frame (this will stop further acquisition). 
*/ 


if (!nocam) 
wn_title(wn," ACQUIRE IMAGE - SNAP MODE "); 


while(1) ( 
int done = FALSE; /* check whether image snapped */ 
wn_clr(wn); 
if (!nocam) ( /* using camera */ 
v_kflush(); /* clear keyboard buffer */ 


wn_printf(wn,"\n\n\tPress {SPACEBAR) to snap an image"); 

wn _printf(wn,"\n\t repeat until satisfied"); 

while ( (c = getch()) != ENTER) ( 
wn_puts(wn,6,4,"Wait..."); 


waitvb(); /* Wait for vertical blanking */ 
snap(WAIT) ; /* acquire a frame */ 
wn _puts(wn,6,4,"Done! "); 
wn_puts(wn,12,4,"When Satisfied, press {Enter)"); 
done = TRUE; /* flag that image taken */ 
) /* end while */ 
if (!done) { 
waitvb(); /* Wait for vertical blanking */ 
snap(WAIT); /* turn grab off */ 
) /* end if !done */ 
crystal(); /* internal sync for stability */ 
setcamera(1); /* disconnect camera 0 */ 


) /* end if !nocam */ 

/* Clear unwanted areas */ 

wn_clr(wn) ; 

wn_printf(wn,"\n\n\tPress {SPACEBAR} to crop image");/* line 2 */ 
wn_printf(wn,"\n\t one vertical line at a time"); 


if (DF _ LM) { 
if ( LT_MARGIN > IXS ) { /* auto crop left margin */ 
wn_printf(wn,"\n\tCropping left edge..."); 


for (i = IXS; i < LT_MARGIN; i++) 
viclear(i,IYS,NROW, BLACK LEVEL); 


} 
lm = LT_MARGIN; 

) 

else { /* let user crop left margin */ 
lm = IXS; /* start from left edge */ 
wn_puts(wn,5,2,"Left margin: "); /* lines */ 
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v kfiushi); /* clear keyboard buffer */ 
wn printf(wn,"\n\n\tWhen done, press [Enter]"); /* line 7? */ 
while ( (c = getch()) != ENTER) { 
vlclear(lm, IYS,NROW, BLACK LEVEL); /* clear one vertical line */ 


wn_locate(wn,5,15); /* line 5 */ 
wn_printf(wn, "23d", 1m); /* display current left margin */ 
if (lm >= RT_MARGIN) /* right edge reached */ 
break; /* > stop */ 
else 
lm++: /* do next line */ 
} /* end while */ 
LT_MARGIN = 1m; /* update global */ 
} /* end if-else */ 
Ti CDE ANA 
if(RT MARGIN < NCOL) { /* auto crop right margin */ 
wn printf(wn,"\n\tCropping right edge..."); 


for (i = NCOL; i > RT MARGIN; i--) 
vlelear(i, IYS,NROW, BLACK_LEVEL); 
rm = RT_MARGIN; 


} 

} 

else { /* let user crop right margin */ 
rm = NCOL; 
v_kflush(); /* clear keyboard buffer */ 
wn_puts(wn,5,2,"Right margin: db E /* line 5 */ 


wn_printf(wn,"\n\n\tWhen done, press [Enter]"); /* line 7 */ 
while ( (c = getch()) != ENTER) { 
viclear(rm,0,NROW, BLACK LEVEL); /* clear one vertical line */ 


wn_locate(wn,5,16) ; /* line 5 */ 
wn_printf(wn,"%3d",rm); /* display right margin */ 
if (rm <= LT MARGIN) /* left edge reached */ 
break; Te Step, Vi 
else 
mn; /* do next line */ 
} /* end while */ 
RT MARGIN = rm; /* update global */ 


) /* end if-else */ 
if ( (!DF_LM) && (!DF_RM) ) { 
wn _printf(wn,"\n\n\tSatisfied with result [Y]?");/* line 9 */ 


v_kflush(); /* clear keyboard buffer */ 
c = getch(); /* get user response */ 
if ( c !@°"°N’ &6& ¢ !* *n’ ) J* not no */ 

break; /* done, get out of loop */ 


/* prepared to redo */ 
if (nocam) 


initluts(): /* restore image */ 
else { 
setcamera(0); /* connect camera 0 */ 
extsync(); /* external sync */ 
waitvb(); 
grab(NO WAIT) ; /* ensure grab mode ON */ 
grab(NO WAIT); /* restore image and start over */ 
} /* end if-else */ 
} 
else 


break; /* get out of while loop */ 
} /* end while(1) */ 
wn_clr(wn); 
/* Complement Image in frame memory */ 
if (DF_INVERT) { 
wn_printf(wn,”\n\tComplementing Image..."); 
complement (IXS, 1TYS,NCOL, NROW) ; 


} 
else { 
wn_printf(wn,"\n\n\tComplement Image [Y}?"); 
v_kflush(); /* clear keyboard buffer */ 
c = getch(); 
if ( (ct='N’) & (c!='n') ) { 
wn_printf(wn,"\n\tComplementing Image..."); 
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complement(IXS, IYS,NCOL,NROW) ; 


) 

} /* end if-else */ 

/* 

* Save image 

Sf 
wn_printf(wn,"\n\n\tUpdating frame memory..."); 
maplut(GRNLUT,0,IXS,IYS,NCOL,NROW); /* update frame memory */ 
v_kflush(); /* clear keyboard buffer */ 


wn_printf(wn,"\n\n\tSave image[Y} ?"); 
c = getch(); 


if C (c =n ES (c P= Ne 
putim(wn,0); /* Save image with .img extension */ 
io 
* Update session file 
*/ 
wn_printf(wn,"\n\n\tUpdating session file..."); 


if (( fp = fopen(session,"a") ) == NULL) { 
wn_printf(wn,"\a\n\n\tUnable to open session file %s",session) ; 
wn_printf(wn,"\n\n\tPress any key to continue"); 
getch(); 
} 
else { 
if (nocam) 
fprintf(fp,"\nmACQUIRE: Image read from Filename: %s", filename) ; 
else { 
if ( c !=—"n’ && c != “ND 
fprintf(fp,"\nACQUIRE: Image saved to Filename: žs",filename); 
} /* end if-else */ 
fprintf(fp,"\n Comment: %s",comline) ; 
fclose(fp); /* close session file */ 
} /* end if-else */ 
wn_close(wn); 
return(0); 


} 


/* end of file ACQUIRE.C ”/ 
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RC RRARRRAARARARARARARARARARARARRRRRARARRRARRARRARARARARARARR RRA RARA RA AAA 
PRETEENAME ¿“CLIPC 


* CALLED BY: semex main() and calls clipmain() 
LAST MODIFIED : 12 Mar 91 by LEE YEAW-LIP 


cr me mm nt we re ce ee ee re ee ee ee i a a A A A A A E e a A A A A A A A Ie 


+ 


+ 


* PURPOSE : This program clips (thresholds) the image on screen by 

Y taking an operator input value and forcing all image pixel 
* values above the threshold value to BLACK_LEVEL and those 
A below the threshold value to WHITE_LEVEL. 

x Latest version has an automatic threshold capability 

” called autoclip(). It determines the background value 

di for predefined 7 regions in the image and uses the 

2 darkest value. 

k ( Background == WHITE_LEVEL ; Feature == BLACK_LEVEL ) 


Fe te te We de te de Fe de te te Fe te te de de de de de de de de de de de de de de de de ve Fe ve ve Ye de de de de de de te de de de te de de We de de de he de he de dek de he de teh tehe teth 


+ 


nT 

#include “global.h" 

clipmain() 

{ 
extern int clip(WINDOWPTR w,int t); /* function prototype */ 
WINDOWPTR wn; /* window pointer */ 
int watrib, batrib; /* window, border and */ 
unsigned fatrib; /* field attributes */ 
int tval; /* threshold value */ 
char c; /* scratch */ 


watrib = v_setatr(WHITE,BLUE,0,0); /* window color */ 

batrib = v_setatr(RED,WHITE,0,BOLD);/* border color */ 

wn = wn_open(500,8,13,50,10,watrib,batrib); /* open a dialog window */ 

if (!wn) { 
printf("\a\n\tUnable to open window. Aborting..."); 
exit(1); 

ł 

wn_title(wn," CLIPPING IMAGE "); 

if ( (fp = fopen(session,"a")) == NULL) { 
wmm_printf(wn,"\a\n\n\tCANNOT open session file Zs", session); 
wn_printf(wn,"\n\n\tPress any key to continue”); 


getch(); 
wn_close(wn); /* close dialog box */ 
return(1); /* premature termination */ 
} 
if (LOAD_RAW) { /* auto loading of RAW image */ 
wn_printf(wn,"\n\n\tLoading RAW image"); 
if(getim(wn,0)) /* default ext is .img */ 
goto err; /* error detected, terminate */ 
else 
fprintf(fp, "\nCLIP:\tImage from Filename: %s", filename) ; 
} 
else { /* ask before loading */ 


v_kflush(); /* empty keyboard buffer first */ 
wn _printf(wn,"\n\n\tLoad Image from disk [N]? "); 
wn_printf(wn,"\n\tPress [ESC] to QUIT”); 


c=getch(); 
if (c == ESC) goto err; IA guitas 
aie =n] cm ty?) 
if(getim(wn,0)) /* default ext is .img */ 
goto err; /* error detected, terminate */ 
else 


fprintf(fp,"\nCLIP:\tImage from Filename: %s", filename); 
} /* end if-else */ 


/* clip the image */ 

wn_clr(wn); 

wn_printf(wn,"\n\n\tDetermining threshold..."); 

tval = autoclip(); /* Use Auto Threshold */ 
fprintf(fp,"“\n\tAuto Threshold: %3d",tval); 

tval = clip(wn,tval); /* Modify Threshold */ 
fprintf(fp,"\tUser Threshold: 23d", tval); 
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wn_clr(wn); /* clear window */ 
wn_printf(wn,"\n\n\t1l: RESTORE image to original and abort"), 
wn_printf(wn,"\n\t2: SAVE the modified image”), 

wn_printf(wn,"\n\t3: EXIT without saving”), 

wn_printf(wn,"\n\n Select option by NUMBER [3]: "); 

v_kflush(); /* empty keyboard buffer first */ 
switch (c=getch()) 

{ 


case ‘1’: 
goto errl; /* restore and quit */ 
case '2': /* update and save */ 


wn_printf(wn,"\n\n\tUpdating frame memory”); 
maplut(GRNLUT,0,IXS,1YS,NCOL,NROW) ; 


putim(wn,1); /* save image with default ext .iml */ 
break; 

case ’'3’; /* update without saving */ 

default: 


wn_printf(wn,"\n\n\tUpdating frame memory"); 
maplut (GRNLUT,0,IXS,IYS,NCOL ,NROW) ; 


break; 

) /* end switch */ 
linlut(GRNLUT,O); /* restore GREEN LUT, bank 0 */ 
fclose(fp); /* close session file */ 
wn_close(wn); /* close window */ 
return(0); /* terminated properly */ 
err: /* premature termination sequence */ 
fprintf(fp,"\n*** CLIP ABORTED **+:); 
errl: 
linlut(GRNLUT,O); /* restore GREEN LUT, bank 0 */ 
fclose(fp); /* close session file */ 
wn_close(wn) ; /* close window */ 
return(1); 

} ` 

/* 


* Clip (Threshold) routine creates a binary image. Pixel values below ’val’ 
* are changed to BLACK LEVEL while all others are changed to WHITE LEVEL 

* Returns selected threshold. 

* NOTE: The change is not permanent as frame memory is not updated. 


x The calling program needs to call maplut() to do this. 
nl 
clip(WINDOWPTR w, int val) /* passes a window handle */ 
{ 
char c; /* scratch */ 


wn_puts(w,3,5,"THRESHOLD LEVEL: “aye 
wn_printf(w, "23d" ,val); 
wn_puts(w,6,5,"Use [+] and [-] keys to adjust”); 
wn_printf(w,"\n\n\tPress [ENTER] when done"); 
v_kflush(); /* empty keyboard buffer first */ 
while ((c=getch()) != ENTER) 
{ 
wn_locate(w,3,24); 
if ( (c=='+’) && (val < HIGH) ) 


val += 2; /* increase threshold */ 
if ( (c==’-’) && (val > LOW) ) 
val -= 2; /* decrease threshold */ 


wn_print£(w,"23d",val); 
threshold(GRNLUT,0,HIGHEST,val);/* threshold GREEN LUT */ 
) /* end while */ 
return val; /* all OK and done */ 
} /* end clip */ 


pe 

* Autoclip function — This function samples the background over 7 regions 
* of the image and determines the darkest value which 
* is returned as the threshold value. The regions are 
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/ 
autoclip(void) 
{ 
int wol ALL. ASe AEC. xAsr’,AEr: j™ x positions */ 
int Ya = 1; /* y position. for region 1 “/ 
int Yb = 80; /* y position for regions 2 and 3 */ 
int Yc = 240; [* y position for region 4 */ 
int Yd = 478; /* y position for regions 5, 6 and 7 */ 
int width = 20; ¡* width ot regions */ 
int’ CENTER = 245: /* center of image */ 
int tval;: /* threshold value */ 
intii.: /* scratch */ 
mete LT MARGIN +91; /* left edge of regions 1 and 2 */ 
KEL = XS1 + width; /* end of regions 1 and 2 */ 
noc = CENTER; /* left edge of central regions */ 
XEc = CENTER + width; /* end of central regions */ 
XSr = RT_MARGIN - 1 - width; /* left edge of regions 6 and 7 */ 
XEr = XSr + width; /* end of regions 6 and 7 */ 
/* determine background intensity for each region */ 
Eval = WHITE LEVEL; /* initial threshold value */ 
tval = findthd(XSc,XEc,Ya,tval); /* find threshold for region 1 */ 
tval = findthd(XS1,XEl,Yb,tval); /* iind threshold for region 2 */ 
tval = findthd(XSr,XEr,Yb,tval); /* find threshold for region 3 */ 
tval = findthd(XSc,XEc,Yc,tval); /* find threshold for region 4 */ 
tval = findthd(XS1,XEl,Yd,tval); /* find threshold for region 5 */ 
tval = findthd(XSc,XEc,Yd,tval); /* find threshold tor region 4 */ 
tvale= findthd(xSr,AEr, Yd, tval): /* find threshold for region 5 */ 
threshold(GRNLUT,O,HIGHEST,tval); /* threshold GREEN LUT */ 
return tval; /* return threshold */ 
} 
/* 
* determine background of region defined by the parameters passed 
* returns the darkest background level 
s 
findthd(int xs, int xe, int y, int T) 
{ 
int pixval; /* pixel value */ 
int. i; /* scratch counter */ 
for (i=xs; i<xe; i++) { 
pixval = rpixel(i,y); /* read pixel value */ 
if( pixval > GRAY ) /* light pixels assumed to be */ 
T = min(T, pixval) ; /* background, take darkest */ 
} 
return T; /* return new threshold */ 


} 


/* end of file CLIP.C */ 
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* FILENAME : TAG.C 
* DEPENDENCIES : called by SEMEX main() and calls semio functions 


LAST MODIFIED : 12 Mar 91 by LEE YEAW-LIP 


== === == —_—_ — o — A A ee ee ae a ei ea e e a e e A e e e 


+ 


> 


* PURPOSE: Labels and identifies each feature in the image. Reads pixel- 
* by-pixel, left to right, top to bottom, and assigns a unique 
* ID (fid) number to each feature (agglomeration of pixels) so 
X that they can be processed by size() or saved for later 
e processing. If there are more than 254 features, a group ID 
w (gid) is also given. The tag() function requires a binary 
* image obtained with the clip() function where 
z Background == WHITE_LEVEL and Feature == BLACK_LEVEL 
Ye ve ve ve ve te Ye Ye Ae Ye Ye we e Ye vr Ye Ye Ye Ye Pe Ye Ye Ye Ye Ye ve Ye Fe we Ye Yr Ye Ye Ye we ve Ye de Ye we Ye Wie de Ye Fe Ye Fe Ye Ye Ye Ye Ye vr Ye Ye Ye vr Ye Ye Fe Fe Ye vr ve Fe fe ve vr ve ve ve 
ef 
#include "global.h" /* required for all SEMEX files */ 
#include <time.h> /* required for time functions */ 
static int fid; /* feature indices */ 
static int gid; /* gid = fid\HIGH */ 
static int maxfl; /* max feature size */ 
tagmain( ) 
{ 
extern int tag(WINDOWPTR w); /* prototype declaration */ 
WINDOWPTR wn; /* window pointer */ 
int watrib, batrib; /* window, border and */ 
char c; 
watrib = v_setatr(WHITE,BLUE,0,0); /* window color */ 
batrib = v_setatr(RED,WHITE,0,BOLD); /* border color */ 
wn = wn_open(500,8,13,50,10,watrib,batrib); /* open a dialog window */ 
if (!wn) { 
printf("\n\tUnable to open window. Aborting..."); 
exit(1); 
) 


wn_title(wn," TAGGING IMAGE ”); 

if ((fp = fopen(session,"a")) == NULL) { 
wn_printf(wn,"\a\n\n\tUnable to open session file Zs", session); 
wn_printf(wn,"\n\n\tPress any key to continue"); 


getch(); 
wn_close(wn); /* close dialog box */ 
return(1); /* premature termination */ 
} 
if (LOADCLIP) { /* load CLIPPED image without asking */ 
wn_printf(wn,"\n\n\tLoading CLIPPED image”); 
if(getim(wn,1)) /* default ext .iml */ 
goto err; /* error, don’t continue */ 
else 
fprintf(fp, "\nTAG:\tImage from Filename: Zsin",filename); 
} 
else if (!DO_SEQ) { /* not processing frame memory */ 
v_kflush(); /* empty keyboard buffer first */ 
wn_printf(wn,”\n\n\tLoad CLIPPED image from disk [N]?"); 
wn_printf(wn,"\n\tPress (ESC) to QUIT"); 
c=getch(); 
if( c == ESC) goto err; /* quit */ 
1£ A | y e) 
if(getim(wn,1)) /* default ext .iml */ 
goto err; /* error, don’t continue */ 
else 
fprintf(fp,"\nTAG: \tImage from Filename: %s\n", filename): 
} 
else fprintf(fp,”"\nTAG:"); 
| ja 
* perform feature extraction and tagging 
my 
if(tag(wn)) { /* if error, don’t save */ 


wn_printf(wn,"\n\n\tPress Any Key to continue"); 
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getch(); 


goto err; /* terminate */ 
} 
v_kflush(); /* empty keyboard buffer first */ 
wn_printf(wn,"\n\n\tSave image to Disk File [N]?"); 
c=getch(); 
if (c== ESC) goto err; e Quit sy 
Ae e y j] (cc == y 09 

putim(wn,2); /* save with default ext .im2 */ 
fclose(fp); /* close session file */ 
wn_close(wn); /* close window */ 


return(0); 


err: /* error condition */ 
Perantte( fp,” \n*** TAG ABORTED ***"); 
fclose(fp); 
wn_close(wn); /* close window */ 
return(1); /* signal back error */ 

} 

/* 

* image feature identification and tagging algorithm 

ay 

tag (WINDOWPTR w) /* window handle */ 

{ 
extern void tagrowO(WINDOWPTR w); /* function prototype */ 
extern void tagrows(WINDOWPTR w); /* declarations */ 
extern void checkmerge(WINDOWPTR w), {* ditto */ 
extern int tagmerge(WINDOWPTR w); }* ditto */ 
extern void step(void); PUEBLO n 
extern clock _t clock(void); /* ditto 5/ 
clock_t start; /* timing variable */ 
float elapsed; ¡+ aitto */ 
unsigned fatrib; /* field attribute */ 
char ibuf(10); /* scratch string */ 
Static char *emsg = "must be between 2 and 255. Press any key", 
fatrib = (BLUE<<4) | WHITE | BOLD; /* field color */ 
wn_clr(w); /* clear window */ 
DEICIDE SIZE) { 

itoa(OVERSIZE,ibuf,10); /* convert to string */ 
wn_gint(XEQ,NFRM,NFLD,w,2,1,"Eliminate features larger than: ", 
fatrib,’_’',&maxfl,3,2,255,ibuf,NSTR,emsg); 
) 
else /* use default max */ 
maxfl = OVERSIZE; /* feature size */ 

start ™=' clock(); /* start timer */ 
wn_printf£(w,"\n\n\tTAGGING FEATURES in Progress\n"); 
fid = LOW; /* initialize feature count */ 
gid = Q; /* initialize group count */ 
wn_wrap(w, TRUE) ; /* word wrap ON */ 
tagrow0(w); /* tag top row */ 
tagrows(w); /* tag subsequent rows */ 
fidget eid*HIGH ~ 1; /* get total features */ 


if (fid <= LOW): { 
wn_printf(w,"\a\n\n\tOverflow has occurred. Fid= Zd", fid); 
wn_printf(w,"\n\tAborting Feature Extraction”); 


return( 1): 
) 
checkmerge(w); /* check merging window size */ 
if (tagmerge(w) ) /* merge joined features */ 
return(1); /* signal back error */ 


wn_printf(w,”\n\tFEATURE COUNT: %d", fid); 

elapsed = (float) (clock()-start) / (float) CLK_TCK; 
wn_printf(w,”\n\n\tElapsed Time: %.1f seconds”, elapsed); 
fprintf(fp,"\tTagged Zd features in 2.1f seconds”",fid,elapsed) ; 
maplut (GRNLUT,0,IXS,IYS,NCOL,NROW) ; /* update frame memory */ 
return(0); /* signal OK back */ 


1r 


) 


os 
* Tag the top row of the image 
* Note that fid cannot exceed HIGH. This is taken care of by step() 


*/ 
void 
tagrow0(WINDOWPTR w) 
{ 
extern void step(void); /* prototype definition */ 
register x, rp; 7* scratch */ 
for( x = LT_MARGIN ; x < RT_MARGIN ; xt+ ) /* start from left edge */ 
{ 
if( rpixel(x,1YS) == WHITE LEVEL ) 
continue; /* skip background */ 
/* from here on, pixel is not background */ 
if ( (x != LT MARGIN) && ((rp=rpixel(x-1,IYS)) != WHITE_LEVEL) ) { 
wpixel(x,IYS,rp); /* west pixel occupied, adopt ID */ 
continue; 
} 
wpixel(x,IYS,fid); /* no neighbor, give new ID */ 
wn printf(w,"\r\tz5a",fid+tgid*HIGH); 
step(); /* increment feature counter */ 
} /* end for x */ 
} 
/* 


* Tag subsequent rows in the image 
* Note that fid cannot exceed HIGH. 
oy 

void 

tagrows(WINDOWPTR w) 

{ 


This is taken care of by step() 


extern void step(void); 


register xl; (~S scratch =/ 


register x, y; /* current pixel location */ 
INC LB: /* pixel value */ 
for( y = 1 ; y < NRON ; y++ ) /* start from top edge + 1 */ 
{ 
for( x = LT_MARGIN ; x < RT_MARGIN ; xt+ ) /* start from left edge */ 


{ 

if( rpixel(x,y) == WHIIE LEVEC) 
continue; /* skip background */ 

/* from here on, pixel is a feature */ 

if( (x != LT_MARGIN) && ((rp=rpixel(x-1,y)) != WHITE_LEVEL) ) { 
wpixel(x,y,rp); /* west pixel occupied, adopt ID */ 
continue; 

} 

if( (rp=rpixel(x,y-1)) != WHITE LEVEL ) { 
wpixel(x,y,rp); /* north pixel occupied, adopt ID */ 
continue; 

) 

wpixel(x,y,fid); /* no N/W neighbors, give new ID */ 

/* Check rest of row for connectivity */ 


x1. =x. + 1: /* start with east neighbor */ 
while(1) 
{ 


if( rpixel(xl,y) == WHITE_LEVEL ) { /* end found */ 
wpixel(x,y,fid); /* assign ID */ 
wn_printf(w, "\r\t%5d", fid+gid*HIGH); 


step(); /* increment feature counters */ 
break; /* done */ 

) 

if( (rp=rpixel(x1,y-1)) != WHITE_LEVEL ) { /* merge */ 
wpixel(x,y,rp); /* use north pixel’s ID */ 
break; 

} 
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xt /* move further east */ 
) /* end while */ 
} ei Mend for x7 
) /* end for y */ 
} 


pe 


* Checks reliability of merge algorithm by comparing the maximum feature 


* length, maxfl, and the density of the features 


ay 
void 
checkmerge(WINDOWPTR w) 
{ 
int safesize; /* safe max window size */ 
safesize = (int) (NRON /( (int)(fid/(HIGH+1)) + 2)); 
if (maxfl > safesize) { 
wn_clr(w); /* declutter the screen */ 
wn_printf(w, "lalnitSizing window cannot support Zd features”, fid); 
wn_printf(w,"\n\tShrianking window from Zd to 2d pixels", maxfl,safesize); 
wn_printf(w,"\n\tAny larger features will be removed”); 
wn_printf(w,"\n\tIf this happens, Clip the image again"), 
wn_printf(w,"\n\tat a lower threshold so as to reduce”); 
wn_printf(w,"\n\tthe number of features."); 
wn_printf(w,"\n\n\tPress any key to continue”); 
maxfl = safesize; 
OVERSIZE = safesize; 
getch(); 
} 
else wn_printf(w,"\n\tLargest permissible feature is %d pixels",safesize); 
) 
p* 


* The following algorithm to merge joined features will fail if 
* maxfl is set too large and there are many (>254) small features 
* bunched together. i 


vy 

tagmerge(WINDOWPTR w) 

{ 
register xl,yl; /* in-window variables */ 
int x.y: /* current pixel indices */ 
int id; /* current pixel ID */ 
int nid; /* ID of pixel above */ 
int xleft,xright,ytop,ybot; /* sizing window limits */ 
int found, merge_done, nmerged = 0; /* tracks merging */ 


wn_printf(w,"“\n\tCombining joined features...\n"); 
for( y = 1 ; y < NRON ; y++ ) /* start from row 1 */ 
{ 
ytop = y ~ maxfl; /* set row search limits */ 
ybot = y + maxfl; 
if( ytop < 1) ytop = 1; 
if( ybot > NROW ) ybot = NROW; 


for( x = LT_MARGIN ; x < RT_MARGIN ; xt?) /* left edge to right */ 


{ 

id = rpixel(x,y); /* get current pixel ID */ 

nid = rpixel(x,y~-1); /* get ID of pixel above */ 

if( (id == WHITE_LEVEL) || (nid == WHITE_LEVEL) | | 
(id == nid) ) /* skip if background or */ 
continue; /* part of same feature */ 

/* Joined features exists. Merge */ 

found = FALSE; /* exist but not found */ 

xleft = x - maxfl; /* set column search limits */ 


xright = x + maxfl; 
if( xleft < LT_MARGIN) xleft = LT_MARGIN; 
if( xright > RT_MARGIN ) xright = RT_MARGIN; 


for ( yl = ytop + yl < ybot. ; yit+ ) 1 /* scan vertically */ 


merge_done = TRUE; /* assume merged */ 


ELS 


for ( x1 = xleft: xl -< xright; x1++*+y /* scanshorizontall, 7, 


if ( rpixel(xl,yl) == id tt /* same feature */ 
found = TRUE; /* found the culprit */ 
wpixel(xl,yl,nid); /* change its ID */ 
merge done = FALSE; /* assumption wrong */ 
} 


if (found && merge done) break; /* culprit found and merged */ 
) £* end tor yl */ 
/* one feature has been merged */ 
wn_printf(w,"\r\tZ%Z5d", ++nmerged); 
fidas. /* decrement feature count */ 
if (fid < LOW) { /* check */ 
wn_printf(w,"\a\n\n\tUnderflow has occurred. Fid = Zd", fid); 
wn_printf(w,"\n\tAborting Feature Extraction"); 
return(1); /* signal error back */ 
} 
} /* end for */ 
) /* end for y */ 


TOTAL = (long)fid; /* save the grand total */ 
return(0); /* signal OK back */ 
) 
/* 
* step() allocates a unique pair of numbers to each feature, namely 
* a feature ID number (fid) and a group ID number (gid). Allocation 
* is based on the limitation that fid must not exceed 8 bits. Therefore 
* Increment fid until it reaches HIGH 
* then reset it to LOW and increment gid 
* This is necessary as pixel value is 8 bits only 
a 
void 
step(void) 
{ 
tig: 
if( fid > HIGH ) { 
fid = LOW; 
eidir: 
} 
} 


/* End of file TAG.C */ 
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* FILENAME : SIZE.C 


CALLED BY: semex main() and calls semio.c functions 


even) MODIFIED : 12 Mar 91 by LEE YEAW-LIP 
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PURPOSE : This routine uses an existing or saved image that has 
been clipped with clip() and tagged with tag(). The output 
of this program is a tabular output of the calculated 
area, X-Chord, and Y-Chord that is suitable for 
input to a statistical analysis program. 

NOTES: : TOTAL (no of features) is automatically set by TAG. 
Background == WHITE_LEVEL 
Feature == BLACK LEVEL 
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ny 
#include "global.h" /* required by all SEMEX files */ 
#include <time.h> /* required for time functions */ 
#include <string.h> /* required for string functions */ 
static float yscale; /* vertical scale factor */ 
static int fid: /* track number of features */ 
static int maxfl, minfl; /* tracks max and min features */ 
static int xmax, ymax; /* max feature lengths */ 

static int xmin, ymin; /* min feature lengths */ 

static long minarea, maxarea; /* smallest and largest areas */ 
static long *“aptr; /* pointer to area store */ 
static int *xptr, *yptr; /* pointer to x,y stores */ 


sizemain( ) 


{ 


extern int size(WINDOWPTR w); /* prototype declaration */ 
WINDOWPTR wn; /* window pointer */ 

int watrib, batrib; /* window, border and */ 
unsigned fatrib; /* field attributes */ 

char c; ; 

watrib = v_setatr(WHITE,BLUE,0,0); /* window color */ 

batrib = v_setatr(RED,WHITE,0,BOLD); /* border color */ 


wn = wn_open(500,8,12,50,10,watrib,batrib); /* open a dialog window */ 
if (!wn) 
{ 
printf(’\n\tUnable to open window. Aborting..."); 
exit(1); 
} 
wn_title(wn," SIZING FEATURES "); 
if ((fp = fopení(session,"a")) == NULL) { 
wn_printf(wn,"\a\n\n\tUnable to open session file 2Zs" session); 
wn_printf(wn,”"\n\n\tPress any key to continue”); 


getch(); 
wn_close(wn) ; /* close dialog box */ 
return(1); /* premature termination */ 
) 
if (LOADTAG) ( /* load TAGGED image w/o asking */ 
wn_printf(wn,"\n\n\tLoading TAGGED image"); 
if (getim(wn,2)) /* iff error don’t size */ 
goto err; /* terminate */ 
else 
fprintf(fp,"\nSIZE:\tImage from Filename: Zs\n",filename); 
) 
else 7if (IDO SEQ) 1 /* not processing frame memory */ 


wn_printf(wn,"\n\n\tLoad TAGGED image from disk [N]?”); 
wn_printf(wn,”\n\tPress [ESC] to QUIT"); 


v_kflush(); /* empty keyboard buffer first */ 
c = getch(); 
if( c == ESC) goto err; It quit */ 
if(c == 'Y'|| c = 'y”) /* read image with ext .im2 */ 
if (getim(wn,2)) /* if error don't size */ 
goto err; /* terminate */ 
else 
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fprintf(fp,"\nSIZE:\tImage from Filename: 2Zsin", filename); 


} 
else fprintf(fp,”\nSIZE:"); /* process frame memory contents */ 
/* feature sizing */ 
if(size(wn)) { /* if error don't save */ 
wn_printf(wn,"\n\n\tPress any Key to continue”); 
getch(); 
goto err; /* terminate */ 


) 
n_printf(wn,"ininitSave data [Y]?"); 


v kflush(); /* empty keyboard buffer first */ 

c= getch(); /* get user response */ 

if ( (ct=’N’) ££ (c!="n') ) outdata(); /* display data */ 

fclose(fp); /* close session file */ 

/* free dynamically allocated memory */ 

free(xptr); 

free(yptr); 

free(aptr); 

v_kflush(); /* empty keyboard buffer first */ 

wn_printf(wn,"\n\n\tSave image to Disk File [N]?"); 

c=getch(); /* get user response */ 

ADS O 
wn_printf(wn,"\n\n\tUpdating frame memory..."); 
maplut(GRNLUT,0,IXS,IYS,NCOL,NROW); /* update frame memory */ 
putim(wn,3); /* save image with ext .im3 */ 

) 

wn_close(wn); /* done, close window */ 


return(0); 


err: /* terminate phase */ 
fprintt (fp, \n*** SIZE ABORTED “77 
fclose(fp); /* close session file */ 
wn_close(wn); /* close window */ 
return(1); 
} 
/* feature sizing algorithm */ 
size(WINDOWPTR w) /* pass window pointer */ 
{ 
extern clock_t clock(void); /* function prototypes */ 
extern int outdata(void); ¡"TALES 
extern void pixelsize(WINDOWPTR w); TARGEO 
clock_t start; /* timing variable */ 
float elapsed; (ARAI CCOMS/ 
unsigned fatrib; /* field attributes */ 
char c; /* scratch */ 
char ubuff[20]; /* scratch string buffers */ 


/* help message strings */ 

static char *hlp1 = "Sets vertical pixel scale factor [1]"; 

Static char *hlp2 = "Features larger than this will be discarded [100]”; 
static char *hlp3 = "Features smaller than this will be discarded [1]"; 
/* error message strings */ 

static char *emsgl = "must be between O and 100. Press any key"; 

static char *emsg2 = "must be between 2 and 256. Press any key"; 

Static char *emsg3 = "must be between 1 and 255. Press any key"; 

static char *emsg4 = "must be between 100 and 9999. Press any key"; 


wn_clr(w); /* clear window */ 
fatrib = (BLUE<<4) | WHITE | BOLD; /* field color */ 
if (DF _VSCALE) 
yscale = VSCALE; /* use default conversion */ 
else { 
sprintf(ubuff,"%27.2£",VSCALE) ; ¡"convert to string */ 
wn_gfloat(XEQ,NFRM,NFLD,w,2,2,"Vertical Scale factor: ", 
fatrib,' ”, £yscale,7,2,0.0,100.0,ubuff,hlpl,emsg1); 
) 
if (DF SIZE) { /* set feature limits */ 


minfl = UNDERSIZE; 
maxfl = OVERSIZE; 
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} 


else { 
itoa(UNDERSIZE,ubuff,10); /* convert to string */ 
wn_gint(XEQ,NFRM,NFLD,w,4,2, "Discard Features SMALLER than. `, 
fatrib, “,émintl, 31,255, ubutt nlp, emsg5); 
itoaCOVERSIZE ,ubuff ,10); /* convert to string */ 


wn_gint(XEQ,NFRM,NFLD,w,6,2, "Discard Features GREATER than: ", 
fatrib,' ,amaxtl,3,2,256, ubutt hlp2 emsgZ):; 
} 
/* Dynamically allocate memory for the TOTAL number of features 
* generated by TAG. If this 1s not available, user option is allowed. 
s 
if (TOTAL <= 0) { 
wn_printf(w,"\n\n\tImage not recently tagged. Continue [N] ?"); 


v_kflush(); /* prevent spurious input */ 
c = getch(); /* get user response */ 
if ((e =="Y’) || (e == *y*)) £ 


itoa(2000,ubuff,10); 
wn_gint(XEQ,NFRM,NFLD,w,B,2, "Max Number of Features Expected: ", 
fatrib,’ ° ,&TOTAL,4,100,9999, ubuff,hlp2, emsg4) ; 
} 
else return(1); /* don’t size */ 
} /* end if */ 
xptr = (int. ”) calloc(TOTAL ,sizeof(int)); 
pb = (int f) calloc(TOIAL,sizeof(int)): 
aptr (long *) calloc(TOTAL,sizeof(long)); 
/* Check for successfull memory allocation ey 
exper || typtr [| raptr ) { 
wn_printf(w,"\a\n\tNot enough Memory to allocate!"); 
wn_printf(w,"\n\tReduce TOTAL = 2d and try again", TOTAL); 
n_printf(w,”"ininlitPress any key to continue”); 
getch(); 
return(1); 


} 
fprintf(fp,”\tVertical scale: 2f pixels/unit length”,yscale); 


fprintf(fp,”"\n\tMin Length spec: 24d\tMax Length spec: 74d" ,minfl,maxfl),; 
/* Begin Sizing Routine X 


start = clock(); /* start timing */ 
wn_wrap(w, TRUE) ; /* turn wordwrap ON */ 
wn_printf(w,’\n\n\tSizing features\n”); 

pixelsize(w); /* size image in pixels */ 


elapsed = (float) (clock()-start) / (float) CLK_TCK; 
wn_printf(w,”\n\n\tElapsed time: 2%.1f£ seconds”, elapsed ); 
fprintf(fp,"\n\tSizing took 2.1f£ seconds”,elapsed); 


return(0); /* signal OK back */ 

} 

p* 

* Size features by pixel count. No scaling factors taken into account. 

e 

void 

pixelsize(WINDOWPTR w) 

{ 
register int xl, yl; /* scratch indices */ 
int y: /* current pixel location */ 
int currval; /* current pixel value */ 
int xleft, ytop, xright, ybot; /* sizing box */ 
int xlen=0, ylen=0; /* zeros feature dimensions */ 
long pixarea=0; /* zeros feature areas */ 
int reject=0; /* zeros rejects */ 
int too_large=0; /* zeros oversized features */ 
int too_small=0; /* zeros undersized features */ 
fid = 0; /* Initialize Feature counter */ 
xmax = 0; /* Max X-Chord reset to zero */ 
ymax = 0; /* Max Y-Chord reset to zero */ 
xmin = NCOL; /* Min X-Chord reset to max */ 
ymin = NROW; /* Min Y-Chord reset to max */ 
minarea = NCOL * NROW; /* Min area reset to max */ 
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maxarea = 0; /* Max area reset to zero */ 


for( y = IYS ; y < NROW ; yt+ ) /* do for all rows */ 

{ 
if( fid == TOTAL ) break; /* Quit when all features sized */ 
ytop = y; /* set up y coordinates */ 
ybot = y + maxfl; /* for sizing box */ 


if( ytop < IYS) ytop = IYXS; 
if( ybot > NROW ) ybot = NROw; 
for( x = LI MARGIN; x < RI MARGIN; xt) /* do for all columns */ 


{ 


currval = rpixel(x,y); /* get current pixel value */ 

if( fid == TOTAL )break; /* all features sized */ 

if( (currval == WHITE LEVEL) || (currval == BLACK_LEVEL) ) 
continue; 


/* pixel is part of feature yet to be sized */ 

xleft = x - maxfl; /* Set up x coordinates */ 

xright = x + maxfl; /* for sizing pox 7/7 

if( xleft < LT_MARGIN) xleft = LT_MARGIN; 

if( xright > RT_MARGIN ) xright = RT_MARGIN; 

for( yl = ytop ; yl < ybot ; yl++ ) /* do for all rows in box */ 


{ 
for( xl = xleft ; xl < xright ; xl++ ) /* do for columns */ 
{ 
if( rpixel(xl,yl) ‘t= currval ) continue; /* skip */ 
pixareat+; /* increment pixel area */ 
xlen++; /* increment x length */ 
wpixel(xl,yl,BLACK LEVEL); /* mark off as counted */ 
) /* end for xl */ 
if( xlen==0 ) break; /* passed bottom edge */ 
if( xlen > xptr[fid) ) /* update max x length */ 
xptr(fid) = xlen; 
xlen = 0; /* reset x length */ 
ylen++; /* increment y length */ 
) /* end for yl */ 
yptr(fid) = ylen; /* store max y length */ 
aptr(fid] = pixarea; /* store max area */ 
ylen = 0; /* ready for next feature */ 


pixarea = 0; 

/* Collect statistics for rejects */ 

if ( xptr[fid] < minfl || yptr[(fid}] < minfl ) { 
too_small++; 


rejectt+, 

xptr(fid)] = 0; /* re-initialize x element */ 
yptr[fid}] = 0; /* re-initialize y element */ 
continue; /* don’t increment fid */ 


if( xptr[fid}) >= maxfl || yptr[fid}) >= maxfl) { 
too_larget+; 


reject++; 

xptr[fid})] = 0; /* re-initialize x element */ 
yptr[fid) = 0; /* re-initialize y element */ 
continue; /* don’t increment fid */ 


) 


/* Calculate Min/Max values */ 
xmax = max( xmax, xptr[fid] ); /* widest feature */ 
ymax = max( ymax, yptr(fid) ); /* tallest feature */ 
xmin min( xmin, xptr(fid) : /* narrowest feature */ 
ymin min( ymin, yptr({fid) ) ; /* shortest feature */ 
maxarea = max( maxarea, aptr[(fid) ) ; /* largest feature */ 
minarea = min( minarea, aptr[{fid] ) ; /* smallest feature */ 
fidit: /* get next feature */ 
wn _printfi(w,"\r\t25da", fid); 
) /* end for x */ 
) /* end for y */ 
TOTAL = (long)fid; /* Record new total */ 
wn_clr(w); /* clear dialog box */ 
wn_printf(w,"\n\n\tSized Zld Features within specifications", TOTAL); 
fprintf(fp,"\n\t2%ld Features were within specifications" ,TOTAL); 


li 
wwe w 
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) 


/* 


Y 


Y 


vr 


if (too_small > 0) { 
wn_printf(w,"\n\tZd Features were less than Zd pixels",too_small,minfl); 
fprintf(fp,"\n\t%d Features were less than 2d pixels",too small,minfl); 
} 
if (too_large > 0) { 
wn_printf(w,"\n\t%d Features were greater than 2d pixels”,too large,maxfl); 
fprintf(fp,"\n\t%d Features were greater than Zd pixels",too_large,maxfl); 
} 
if (reject > 0) { 
wn_printf(w,"\n\tZd Features were REJECTED",reject); 
fprintf(fp,"\n\tZd Features were REJECTED", reject); 


Output routine to display x and y dimensions, and 2 areas 
AREA] assumes feature is elliptical and uses PI*xlen*ylen 
AREA2 converts directly from pixel area to feature area. 


"y 
outdata(void) 
{ 
const float PI4 = 0.785398; /* define pi/4 */ 
FILE “fdata; /* data file pointer */ 
WINDOWPTR w; /* window pointer */ 
int watrib, batrib; /* window, border and */ 
unsigned fatrib; /* field attributes */ 
intaj: /* scratch index */ 
float Cx, Cy, Ca; /* Conversion constants */ 
float fxmax, fymax, fxmin, fymin; /* scaled statistics */ 
float xlen, ylen; /* scaled dimensions */ 
float area, fmaxarea, fminarea; /* scaled areas */ 
char ubuff(20], datafile(20); /* scratch string buffers */ 
char c; 
watrib = v_setatr(WHITE,BLUE,0,0); /* window color */ 
batrib = v_setatr(RED,WHITE,0,BOLD),; /* border color */ 
w = wn_open(800,0,13,50,23,watrib,batrib); /* open a dialog window */ 
if Ciw) 
{ 
printf("\n\tUnable to open window. Aborting...”); 
exit(1); 
) 


wn_title(w,” TABLE OF FEATURE DATA "); 
/* Calculate Conversion Constants required to properly scale and convert 
* pixels to dimensioned units. Depends on VSCALE and ASPECT_RATIO. 
* ASPECT RATIO is defined in SEMEX.C 
ah 
Cx = ASPECT RATIO/yscale; /* x conversion constant */ 
Cy = 1.0/yscale; /* y conversion constant */ 
Ca = Cx*Cy; /* area conversion constant */ 
fprintf(fp,"\n\tConversion constants: Cx=Z%f Cy=Zf Ca=Z%f", Cx, Cy, Ca); 
wn_printf(w,"\n\n\tImage filename: 7s", filename); 
chgext(datafile,filename,".dat”); /* change extension to .dat */ 
fatrib = (BLUE<<4) | WHITE | BOLD; /* field color */ 
wn_gtext(XEQ, NFRM,NFLD,w,4,1,"Save DATA in filename: ", 
fatrib,' ',18,datafile,NSTR,NSTR); 
wn_printf(w,"\n\n\tSending output to FILE 2sin”,datafile); 
if( ( fdata = fopen(datafile,"w") ) == NULL > ( 
wn_printf(w,"\a\n\n\tCANNOT Open output file Zs\n", datafile); 
wn_printf(w,”\n\tPress Any Key to continue”); 


getch(); 

wn_close(w) ; 

return(1); 
) 
forl j 0 j < TOTAL ; jt+ ) { 

xlen = xptr[j}) * Cx; /* true x length */ 

ylen = yptr[3] * Cy; /* true y length */ 

fprintrífdata 260 210.91 210,.3t 268,3 ¿8.3tin , j+1, PI4*xlen*ylen,*(aptr+3)*Ca,xlen,ylen); 
) 
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wn_print£(w, "in ID NO AREA C AREA M X-Chord Y-Chord\n"); 
fort 3 =D 5) 5 TOTAL A 


xlen = xptr{j] * Cx; /* true x length */ 

ylen = yptr{j] * Cy; /* true y length */ 

wn_printf(w, “\n26d 210.3f 210.3f 28.3f 28.3£”, 3+1, Pl4*xlen*ylen, “(aptr+j)*Ca, xlen, ylener 
ff 0 (3220) LOA /* pause on full page */ 


wn_printf(w,"\n\tPress [ENTER] for MORE or [ESC] to QUIT"); 
if ((c = getch()) == ESC) break;/* quit if ESC pressed */ 


it (9005, ClALe = a1) /* not end yet */ 
wn_printf(w,"\n ID NO AREA _C AREA M X-Chord Y-Chordin"); 
} 
} 
fclose(fdata); /* close data file */ 


/* display statistics */ 
xlen = xmax*Cx; 

ylen = ymax*Cy; 

area = maxarea*Ca; 


wn_printf(w,"\n\n AREA M X-Chord YChord y, 
wn_printf(w, "in Max OE ALOE 210.3f", area, xlen, ylen); 
fprintf (fp NTN E AREA M X-Chord Y=chora )- 
fprinti(fp “\n\tMax 210.3€£ ADE 210.3f", area, xlen, ylen); 


xlen = xmin*Cx; 

ylen = ymin*Cy; 

area = minarea*Ca; 

wn_printf(w, "Ain Min IO E ALO Se 210.3f", area, xlen, ylen); 
fprintt( fo, "\n\tMin AOS E £10,305 210.3fin", area, xlen, ylen); 
wn_printf(w,"\n\n\tPress Any Key to continue"); 

v_Ktlush<); 


getch(); 
wn_close(w); . 
return(0); /* signal healthy end */ 


} 


/* End of file SIZE.C */ 
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fp" 
* PILENAME: ANALYZE.C 
* CALLED BY: SEMEX main() 


> 


+ 


+ + + + + + & F 
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LAST MODIFIED: 14 Mar 91 by LEE, YEAW-LIP 

PURPOSE: This set of routines analyzes the data files put out by 
the size() functions. It first allows the user to specify 
the data files and then merges the data from them into an 
array. The volume and data are then calculated and a 
histogram built using histo_vol(). This can be plotted 
using SEM.M inside MATLAB. 


ty Fe ve ve ve ve ve ve vr tr ty ty ve vy ty vy ve ty ty ty ve Hy oy ty ty oe He HY HY oy oy WY hy oy ee Ve ve ve ve oy Ye oy Hy vy vy Wy Wye ve ve ve ty ty ty ve WY ve HY ve ve vy ve ty ty ty Wy vy 


#include "global.h" /* required by all SEMEX modules */ 
include <dos.h> /* dos prototype definitions ”/ 
include <math.h> /* math prototype definitions */ 
#include <string.h> /* string handling prototype defn */ 
typedef struct { /* structured list of file pointers */ 
char name[ 13]; f* string to contain data file */ 

) SLIST; 
typedef struct ( /* Matlab MAT-file structure */ 

long type; /* type */ 

long mrows; /* row dimension */ 

long ncols; /* column dimension */ 

long imagf; /* flag indicating imag part */ 

long namlen; /* name length (including NULL) */ 
} Fmatrix; 
static SLIST “list: /* Define pointer to list */ 
Static char sdate[13]; /* date string */ 
static int fatrib: /* field attribute */ 
static float *cptr; /* pointer to area store */ 
static size_t memsize; /* size of memory allocated */ 
Static int Nfiles; /* number of data files */ 


analyze() 


{ 


extern int merge data(WINDOWPTR wn); /* prototype definition */ 

extern void histo _vol(WINDOWPTR wn); /* ditto */ 

WINDOWPTR wn; /* window handle */ 

int watrib, batrib; /* window and border attributes */ 
char c; /* scratch for user response */ 
watrib = v_setatr(WHITE,BLUE,0,0); /* window attribute */ 

batrib = v_setatr(RED,WHITE,0,BOLD); /* border attribute */ 


fatrib = (BLUE << 4 ) | WHITE | BOLD; /* field attribute defined */ 
wn = wn_open(500,8,13,52,13,watrib,batrib); /* open dialog window */ 
if(!wn) { 
printf("\a\nUnable to open window. Aborting..."); 
exit(1); 
) 
wn_title(wn,” ANALYSE FEATURES "); 
if ( (fp = fopení(session,"a")) == NULL) { 
n_print£(w,”"lalninitUnable to open session file 2s”, session); 
wn_printf(wn,"\n\n\tPress any key"); 


getch(); 
wn_close(wn); 
return(1); 
) 
else /* session file opened */ 
fprintf(fp, "\nANALYZE:Merging data files”); 
if (merge_data(wn)) /* merge data files */ 
goto err; /* terminate if error */ 
histo_vol(wn); /* histogram the data */ 
fclose(fp); /* close session file */ 
wn_close(wn); /* close window */ 


Ln 


return(0); 


err: /* error trapped */ 
fprint£(fp, “Iin*** AÑALYZE ABORIEDOA ; 

fclose(fp); /* close session file */ 
wn_close(wn); /* close window */ 
return(0); 


) 


/* merge data() prompts the user with all files with extension .dat from 
* which to merged data. Selected files are written into a list. 


rp 

merge data(WINDOWPTR wn) 

{ 
extern void get dates(unsigned date); /* prototype definition */ 
extern int extract _data(WINDOWPTR wn, int Ncol); /* ditto */ 
char c; /* scratch for user response */ 
int Area type; /* defines type of area to use */ 
struct find_t d_file; /* structure of data files */ 


n_title(wn," SELECT DATA FILES "); 


Nfiles = 0; /* initialize file counter */ 
wn_clr(wn); /* clear dialog box */ 
_dos_findfirst("*.dat", A NORMAL, &d_ file); /* get first occurrance */ 
get_dates(d_file.wr_date); /* returns date string to sdate */ 
wn_printf(wn,"\n\n\t DATA FILENAME DATE CREATED"), 

n_locate(wn,4,4); /* place cursor at row, col */ 
wn_printf(wn," 212s Zs \n\n\tInclude[Y) ?",d_file.name,sdate), 
v_kflush(); /* empty keyboard buffer */ 

c = getch(); /* get user response */ 


if ( (c==’Y') || (c=="y’) || (c==ENTER) ) ( 
if ((list = (SLIST *)calloc(1,sizeof(SLIST))) == NULL) { /* allocate memory */ 
wn_printf(wn,"\a\n\n\tNo memory to allocate”); 
wn_printf(wn,"\n\n\tPress any key to continue”); 


getch(); - 

return(1); /* signal error back */ 
) 
strcpy(list[Nfiles).name,d_file.name); /* save filename into list */ 
Nfilest+; /* increment file count */ 


} 
/* find the rest of the data files */ 
while (_dos_ findnext(&d_file) == 0) { 


get_dates(d_file.wr_date); /* returns date string to sdate */ 
wn_locate(wn,4,4); /* place cursor at row,col */ 

wn_printf(wn," 212s Zs \n\n\tFile Count: Z3dltInclude[Y] ?",d file.name,sdate,Nfiles); 
v_kflush(); /* empty keyboard buffer */ 

c = getch(); /* get user response */ 


if ( (c=="Y*) || ce=='y') || (c==ENTER)) ( 
if ((list = (SLIST *)realloc((void *)list,(Nfiles+t1)*sizeof(SLIST))) == NULL) ( 
n_print£(wn,"lalninitNo memory to reallocate"); 
wn_printf(wn,"\n\n\tPress any key"); 


getch(); 
goto errl; /* abort */ 
) 
strcpy(list[Nfiles)].name,d_file.name); /* save filename into list */ 
Nfiles++; /* increment file counter */ 
) 
if ( c == ESC ) break; /* enough already */ 


) /* end while */ 
wn_printf(wn,"\n\n 2d datafiles selected in this directory", ,Nfiles); 
wn_printf(wn,"\n Press (Enter) to Histogram data or [Esc] to Abort"); 


v_kflush(); /* flush keyboard buffer */ 
if ( getch() == ESC) /* get user response */ 
goto errl; /* abort */ 


/* allocate memory for first element */ 

memsize = sizeof(float) ; 

if ((cptr = (float *)calloc(1,memsize)) == NULL) { 
wn_printf£(wn,"\a\n\n\tNo memory to allocate”); 
n_print£(wn,"ininitPress any key"); 
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/* 
* 
* 
Y 
ve 
ve 
vr 
= 


vol 


zetcht): 


goto errl; /* signal error back */ 
) 
wn _printf(wn,"ininitUse C)alculated or M)easured area [C] ?"); 
v_kflush(); /* clear keyboard buffer ”/ 
c = getch(); /* get user response */ 
iE (c == M || Ke == om’ y 31 
Area type = 2; /* use measured area AREA M */ 
fprintf(fp,"\n\tExtracting MEASURED Area from"); 
} 
else { 
Area type = 1; /* use calculated area AREA C */ 
fprintf( fp,’ \n\tExtracting CALCULATED Area from"); 
} 
extract _data(wn,Area type); /* extract area from col 2 */ 
Freturn( 0): /* signal back OK */ 
errl: /* termination sequence */ 
free(list); /* deallocate memory */ 
return(1); /* signal back abort */ 


Get date string given a dos date code 


returns date string to sdate declared static 
date format is: Bits 0 - 4 : Day of the month (value between 1 and 31) 


5 - 8 : Month (value between 1 and 12) 
9 - 15: Year since 1980 


d 


get_dates(unsigned date) 


{ 


char “mth; /* month string. */ 

unsigned mm,dd,yy; /* date variables */ 

Static char *month( 12] (“Jan “Feb. Mar’: Apr’, May . Jun, 
"Jul, “Aug” “Sep”, “Oct, NOV”, Dec” ): 


dd = date & Ox001f; /* get bits O to 4 */ 
m = (date & 0x01e0) >> 5; /* get bits 5 to 8 */ 
yy = (date & Oxfe00) >> 9; /* get bits 9 to 15 */ 


sprintf(sdate,"%2u 273s %4u",dd,month(mm-1],yy+1980); 
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/* extract _data() extracts a single column of data from a specified file 


ah 
extract _data(WINDOWPTR wn, int Ncol) 
( 


FILE *dfp; /* data file pointer */ 

fleat, coll >]: /* scratch array */ 

int fid: /* scratch feature counter */ 
imtas /* scratch counter */ 


wn_clr(wn); 
wn_title(wn," 
TOTAL = 0; 
for (i=0; i < Nfiles; it+) { 

if ((dfp = fopen(list[i] .name,‘r’ )) 


EXTRACTING DATA a); 


/* 
/* 


initialize total no of features */ 
read a data file */ 
NULL) { /* read data within */ 


wn_printf(wn,"\a\n\n\tCANNOT Open datafile Z%s",list[i].name) ; 
wn_printf(wn,”\n\n\tPress any key to continue”); 


getch(); 
return(1); 


) 


/* 


error recovery */ 


wn_printf(wn,"\n\t%s", list[i].name); /* display data file */ 
fprintf(fp,”"\n\t\tZs”, list[i].name); 


while (fscanf(dfp,” Zd Zf Zf Zf Zf ", &fid,&col[l)],&col[2],&col[3],&col[4]) 


/* record date file name */ 
!= EOF) 


{ /* read in data */ 
cptr[TOTAL] = col[Ncol]; /* use requested column */ 
TOTAL++: /* increment total # */ 


/* allocate memory for 1 more element */ 


memsize t= sizeof(float); 
if (( cptr = 


wn_printf(wn,"lalninitReallocation Failed. 


/* 


increase memory counter */ 


(float *) realloc((void *) cptr, memsize )) == NULL) { 


Out of memory.”"); 


wn_printf(wn,"\n\n\tPress any key"); 


getch(); 

fclose(dfp); /* close data file */ 

return(1); /* error recovery */ 

} 

} /* end while fscan dfp */ 
fclose(dfp); 
wn_printf(wn,"\n210d extracted. Running Total is Zld", fid, TOTAL); 
fprintf(fp,"210d extracted. Running Total is Zld", fid, TOTAL); 


} /* end for loop */ 
v_kflush); 


/* 


empty keyboard buffer */ 


wn_printf(wn,”\n\n\tPress [Enter] to Histogram or [Esc] to Quit"); 


free(list); /* deallocate memory */ 
if ( getch() == ESC ) return(1); /* quit if ESC pressed */ 
return(0); /* successful */ 


) 


/* histo _vol() uses the merged data from merge_data(), converts it to volume 
* and collates the data into 38 bins so that a histogram is obtained. 
me 


#define BINS 38 /* No of bins */ 
void 

histo_vol(WINDOWPTR wn) 

{ 


extern void savemat(FILE *fptr, int type, char *pname, /* prototype definition */ 
int mrows, int ncols, int imagf, /* for saving data */ 
double *preal, double *pimag); /* in MATLAB form */ 


char histfile[20); /* histogram file */ 


FILES* fn: /* file pointer */ 

double PI = 3.141593; /* define pi */ 

double UPPER = 180.0; /* upper limit of histogram */ 
double STEP = 0.8264; /* same as Malvern MasterSizer */ 
static double volume[BINS] ; /* bins to contain volumes */ 
Static long count[BINS]; /* bins to contain counts */ 

double limit[BINS+1]; /* array of limit values */ 

double tot_vol = 0.0; /* total volume */ 

double toosmall = 0.0; /* bin to contain small particles */ 
int reject = 0; /* track rejects */ 

double diam, vol; /* diameter and volume variables */ 
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double C; e 
double R2; pr 
int 1: {= 
long j; ea 
double Nf, T; pe 
char c; e 


wn_clr(wn); 


wn_title(wn," HISTOGRAM VOLUMES "); 


Cee PI * 4.0 /-3.0; pe 
limit[0] = UPPER, pe 
for (i0; i < BINS; itt) { 
limitit] = limitli] * STEP; i 
volume[i] = 0.0; ee 
count[i] = 0; E 
} 


wn_printf(wn,'"\n\n\tCalculating equivalent 
tor (320: 3 < TOTAL: 3++) { LE 


constant = 4/3pi */ 
scratch = Area/PI */ 
scratch index counters */ 


scratch particle counter */ 
scratch for type conversion */ 
scratch for user response */ 


constant of proportionality */ 
assign upper limit to element 0 */ 


generate limits for each bin */ 
initialize volume bin */ 
initialize count. bin */ 


volumes..."); 
repeat for all particles */ 


/* calculate the equivalent volume assuming a sphere, given area */ 


R2 = cptr[j]/PI:; e 
diam = 2.0 * sqrt(R2); qe 
vol =C * pow(R2,1.5); no 
tot voltei vol: pe 


square of the radius */ 
diameter of particle */ 
volume of particle */ 
accumulate total volume */ 


/* sieve the volumes and collate into the correct bins */ 


if (diam > UPPER) { i 
wn_printf(wn,"”\n\tMassive particle 


too large */ 
Volume=%g, Size=2f",vol,diam); 


fprintf(fp,”\n\tMassive particle Volume=%g, Size=2f", vol, diam); 


continue; 
} 
i = 0; 


while (TRUE) { ee 


if (diam > limit (it+1]) 4 {= 
volume[i] += vol; pe 
countla]++: ye 
break; pa 
) 
else ( 
itt; pr 
if (i >= BINS) { pa 
toosmall += vol; /* 
reject++; {2 
break; pe 


) 
) /* end if-else */ 
} /* end while */ 
} /* end for j */ 


fprintf(fp,”\n\t%d particles < Zf rejected. 


toosmall*100/tot_vol); 
/* display results */ 
wn_printf(wn,"\n BIN UPPER LOWER 
forto; i < BINS; itt) { 


if (limit[i] > 100.0) qn 


mm print£in, “in 23d 23.1£ 25.1£ 210,2£ 


Volume 


repeat until right bin found */ 
found */ 

accumulate volume in bın */ 
increment count */ 

get out of while loop */ 


check next smaller bin */ 
no more bins */ 

put these into toosmall */ 
increment reject */ 

get out of while loop */ 


Volume= Zg or 26.3£22",reject,limit[BINS],toosmall, 


ae Total Vol Counti): 


use one decimal point only */ 
240.2É Ala”, 


} 


¿1+1,limit(i] ,limit(i+1), volume(i), volume(i)”"100.0/tot_vol,count(i]); 
/* use two decimal points */ 
wn_printf(wn,"\n 73d 25.2f 25.2f 210.2f ZO. ok Aol. 
itl,limit(i],limit(it+1) ,volume[i] ,volume(i]*100.0/tot_vol,count[i]); 
if (1210 == 9) { /* 10 bins at a time */ 
wn_printf(wn,"\n\n\tPress (Enter] for More”); 


else 


v_kflush(); /* empty keyboard buffer */ 
getch(); /* wait for user response */ 
if (i =I BINS = 1) { /* not last bin yet */ 


/* generate new screenful */ 
UPPER LOWER Volume 22 Total Vol Count”); 


wn_clr(wn); 
wn_printf(wn,”\n BIN 


} 


/* Print out results */ 


chgext(histfile,session,”.his”); 


/* form histogram file name */ 
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wn_print£(wn,"ininitPrinting results to Zs",histfile); 

fprintf(fp,"\n\tPrinting results to 2s",histfile); 

if ((fh = fopen(histfile,"a")) == NULL) { 
wn_printf(wn,”\a\n\n\tCANNOT Open file %s",histfile); 
wn_printf(wn,"\n\n\tPress any key”); 


getch(); 
goto err3; /* premature termination */ 
} 
fprintf(fh,"\n Zs: Data from 2d data files. Total Particle count: Z1d",histfile,Nfiles, TOTAL); 
fprint£(fh, "in BIN UPPER LOWER Volume 22% Total Vol Count 22 Total”); 
for (1=0: 1 < BINS; i++) /* write to histogram file */ 
if (limit[i]) > 100.0) /* use one decimal point only */ 
Eprint£( th, an 3d 23 ES AO ETONE ZO LE Z51ld 46.2f" ,i+1, LIMA" 
limit[it+1] ,volume[i],volume[i]*100.0/tot_vol,count[i],count[i]*100.0/TOTAL) ; 
else /* use two decimal points */ 
forintfi (fh, \n <30° 25.25 25-22 1052 26520 Asi 46.2f", itl Limita» 
limit[i+1], volume[i], volume[i]*100.0/tot_vol,count[i],count[i]*100,0/TOTAL); 
fclose(fh); /* close histfile */ 


wn_printf(wn,"\n\t%d Files extracted; Particle Count = Zld",Nfiles, TOTAL); 
wn_printf(wn,"\n\tGenerate a MATLAB MAT-file [Y] ?"); 


v_kflush(); /* prevent spurious inputs */ 
c = getch(); /* get user response */ 
if ( (c != 'N’) && (c !=’n’)) { /* generate MAT-file */ 


chgext(histfile,session,".mat"); 
wn_clr(wn); 
wn_gtext(XEQ, NFRM,NFLD,wn,4,4,”MATLAB MAT-filename: ",fatrib, ’_’, 18, histfile,NSTR,NSTR); 
wn_print£(wn,”ininitPrinting MAT-file Zs"”,histfile); 
fprintf(fp,"\n\tPrinting MAT-file 7s",histfile); 
if ((fh = fopen(histfile,"wtb")) == NULL) { 
wn_printf(wn,"\a\n\n\tCANNOT Open file %s",histfile); 
wn_printf(wn,"”\n\n\tPress any key”); 


getch(); 
goto err3; /* premature termination */ 
} 
Nf = (double)Nfiles; /* for proper type casting */ 
T = (double)TOTAL; /* ditto */ 
/* MATLAB variables will be: 
* *Nfiles’,’total’,’limit’, *totalv' and *volume' 
5] 


savemat(fh, 0, "Nfiles”, 1, 1, O, &Nf, (double *)0.0); 
savemat(fh, 0, "total", 1, 1, 0, &T, (double *)0.0); 
savemat(fh, O, "limit", 1, BINS+1, O, limit, (double *)0.0); 
savemat(fh, 0, "volume", 1, BINS, 0, volume, (double *)0.0); 
savemat(fh, 0, "totalv”, 1, 1, 0, £tot_ vol, (double *)0.0); 
savemat(fh, 0, “count”, 1, BINS, 0, (double *)count, (double *)0.0); 
fclose(fh); /* close MAT-file */ 
) 
err3: /* terminate stage */ 
free(cptr); /* deallocate dynamic memory */ 
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1E 


* savemat - C language routine to save a matrix in a MAT-file. 


*MAuUtChOr J.N: Little 11=3-86 

* Adapted by Y.L. Lee 25 Feb 91 
oi 
void 
savemat(fptr, 
PIES =fpctr: 
int type; 


mrows, ncols, 


type, pname, 


int mrows; 
int ncols:; 
int imagf, 
char *pname; 
double *preal; 
double *pimag; 
{ 
Fmatrix x; 
int mn: 


.type = (long) type; 
.mrows = (long) mrows; 
ncols (long) ncols:; 
.imagf (long) imagf; 
-namlen = 
m = x.mrows * x.ncols; 


A 4 A RK OM 


sizeof(Fmatrix), 1, 


sizeof(char), 


fwrite(&x, 
fwrite(pname, 
fwrite(preal, 
if (imagf) { 
fwrite(pimag, 


} 


/* End of file ANALYZE.C */ 


TRUE) 
(ant)x.namlen, 
sizeof(double), mn, fptr); 


imagf, preal, pimag) 


7 
Pe 
/* 
/* 
p* 
Ta 
/* 
qa 
/* 
/* 


(long) (strlen(pname) + 1); 


File pointer */ 

Type flag: O for PC */ 

Add 1 for text variables. 
See LOAD for more info. */ 
row dimension */ 

colum dimension */ 
imaginary flag */ 

pointer to matrix name */ 
pointer to real data */ 
pointer to imag data */ 


ni 


fpEr); 


sizeof(double), m, fptr); 
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* FILENAME =: SEMIC.¢ 


* CALLED BY: Various SEMEX functions 
* LAST MODIFIED : 7 Mar 91 by LEE YEAW-LIP 


* PURPOSE : getim(w,n) — reads in an image from disk 
putim(w,n) - saves an image to disk 
Ye where n = 0 appends .img to filename (raw image) 


+ + + + + F 


1 appends .iml to filename (clipped) 

2 appends .im2 to filename (tagged) 

3 appends .im3 to filename (sized) 
chgext(*fd,*fs,*ext) ~- replaces the extension in the source 


and copies it to the destination. 


NOTE: All these functions have been declared globally 


We ove ve ve ve ve dy ve ve ve eeek We ve ve ve ve ve ve ye ye He He ar 


x 
#include "global.h" 
#include <string.h> 
#include <math.h> 


Static char *digit = "1234567890 
/* 

* Read an image from disk 

* passes a window handle and th 


* OSAMA eee 
a 

getim(WINDOWPTR w,int n) 

{ 


int errval; 

char c, *line2; 

char *margins; 

char *vs; 

char *lm, *rm; 

static char fbuf(MAXFLEN]; 
unsigned fatrib; 


wn_clr(w); 


/* 
/* 
[x 


i: jz 


e default file 


.im2, 3: .im3 


/* 
/* 
/* 
/* 
/* 
/* 
E 


e 


wn_printf(w,"\n\tREAD IMAGE FROM FILE "); 


strcpy(fbuf, filename); 
while(1) 
{ 
switch (n) 
{ 
case 1: 
chgext(filename, 
break; 
case 2: 
chgext (filename, 
break; 
case 3: 
chgext(filename, 
break; 
case 0: 
default: 
chgext (filename, 
break; 
) /* end switch */ 


i 

/* 

/* 
thule amis 
fbuf, .imz ): 
£fbut, .im3"): 


fbuf,".img"); 


fatrib = (BLUE<<4) | WHITE | BOLD; /* 


wn_gtext(XEQ,NFRM, NFLD,w,2,1,"Filename: 


’_', 18, filename 


»NSTR,NSTR); 


required by all SEMEX files */ 
string handling prototypes */ 
math function prototypes */ 


define a digit */ 


extension where 


scratch for errors */ 
scratch */ 

loc of margins in comline */ 
loc of VSCALE in comline */ 
left and right margin loc */ 
filename buffer */ 

field attribute */ 

clear screen */ 

make copy of filename */ 
decide which default ext */ 


to use */ 


/* clipped image */ 


/* tagged image */ 


/* sized image */ 


/* raw video image */ 


field color */ 
” fatrib, 


errval = readim(IXS,IYS,NCOL,NROW,filename,comline); 


if(errval == 0) { 
wn_wrap(w, TRUE) ; 


/* 
/* 


image successfully read */ 
allow for word wrap */ 


wn_printf(w,"\n COMMENTS: \n %s",comline); 
/* check for 2nd comment line */ 


if ((line2 = strrehr(comline,’\n’)) 


t= NULL) ( 


/* Locate the start of the keyword "Margins=" */ 
if ((margins = strstr(line2,"Margins=”")) != NULL) 
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{ /* margins exists, extract them */ 


margins = strpbrk(margins,digit), f* locate starting dizit */ 
LT MARGIN = atoi(margins), /* extract left margin ”/ 
margins = strpbrk(margins,”","); /* locate separator */ 

/* skip over comma and extract right margin “/ 

RT_MARGIN = atoi(+tmargins); /* extract right margın */ 


) 
/* Locate the start of the keywork "VSCALE" */ 


if ((vs = strstr(line2, "VSCALE=")) != NULL) 

{ PY VSCALE. exists. extract it */ 
vs = strpbrk(vs,digit); {* locate starting digit */ 
VSCALE = (float) atof(vs); /* vertical scale factor */ 


} 
) 
wn_printf(w,"\n\n\tPress any key to continue”); 
V KE Lush ():; /* prevent premature keystroke */ 
getch(); 
return(0); /* all is well */ 
} 
else { /* problem */ 
switch(errval) 
{ 
case FILE ERROR: 
Wn _printi(w,”\n\tError opening fileln”); 
break; 
case FORMAT ERROR: 
wn_printf(w,"\n\tUnknown file format\n"); 
break; 
case READ ERROR: 
wn_printf(w,"\n\tError Reading file\n"), 
break; 
default: 
wn_printf(w,"\n\tUnknown Error Z%d\n",errval); 
break; 
) /* end switch */ f 
wn_printf(w,"\a\n\tTry Again [Y]}?\n"), 
Vikflush(); /* empty keyboard buffer first */ 
c = getch(); /* get response */ 
if(c == 'N* || c == 'n') return(1); /* signal problem back */ 
} /* end if-else */ 
wn_clr(w); /* clear window and try again */ 
} /* end while */ 
) 
/* 


* Save an image to disk 
* passes a window handle and file extension where 


* DE l: „iml, 2: a MZ... 33 ime 
= 
putim(WINDOWPTR w,int n) 
{ 
unsigned fatrib; {/* field attribute */ 
int errval; /* scratch for error handling */ 
char c; /* scratch */ 
static char fbuf[MAXFLEN], ext[5]; /* scratch strings */ 
wn_clr(w); /* clear window */ 
wn_printf(w,"\n\tSAVE IMAGE TO DISK "); 
strcpy(fbuf, filename); /* make copy of filename */ 
while(1) 
{ 
switch (n) /* decide which ext to use */ 
{ 
case l: 
strcpyiext, iml ); /* append default ext .iml */ 
chgext (filename, fbuf,ext); 
break; 
case 2: 
strepy(ext,”.im2”): /* append default ext .im2 */ 
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chgext(filename,fbuf,ext); 
break; 
case 3: 
strcpy(ext,".im3"); /* append default ext .im3 */ 
chgext (filename, fbuf, ext); 
break; 
case 0; 
default: 
strcpy(ext,'".img"); /* append default ext .im4 */ 
chgext. (filename, fbhuf, ext); 
break; 
} /* end switch */ 
fatrib = (BLUE<<4) | WHITE | BOLD; /* field color */ 
wn_gtext(XEQ,NFRM,NFLD,w,2,1,"FILENAME: ", 
fatrib,' ',18,filename,NSTR,NSTR); 
/* append image defaults into second comment line */ 
sprintf(comline2," Gain= 23d; Offset= 23d; Margins= 23d, 23d; VSCALE= 27.3f", 
GAIN_LVL,OFFSET_LVL,LT_MARGIN,RT_MARGIN, VSCALE) ; 
wn_wrap(w, TRUE); /* allow for word wrap */ 
wn_printf(w,"\n COMMENTS: \n\nZs",comline2); 
wn_gtext (XEQ,NFRM,NFLD,w,4,1,NSTR,fatrib,’_’,48,comlinel,NSTR,NSTR) ; 
strcpy(comline,comlinel); /* copy first line of comment */ 
strcat(comline,"in"); 
strcat(comline,comline2); /* append additional comments */ 
n_printf(w,"ininininitStoring Image..."); 
errval = saveim(IXS,IYS,NCOL,NROW, COMPRESSION, filename,comline); 
if(errval == 0) ( 
wn_printf(w,"\n\n\tImage successfully SAVED"); 


return(0); /* all is well */ 
} 
else { 
wn_printf(w,"\a\n\tError saving file!!"); 
if(errval == ALLOCATION ERROR) 
wn_printf(w, "\n\tInsufficient Disk Space"); 
if(errval == WRITE_ERROR) 
wn_printf(w,"\n\tError writing file or values"); 
wn_printf(w,"\nin\tTry Again [Y])?"); 
v_kflush(); /* empty keyboard buffer first */ 
c = getch(); 
if(c == 'N’ || c == ’n’) return(1); /* signal problem back */ 
wn_clr(w); 
} /* end if-else */ 
wn_clr(w); /* clear window and try again */ 


} /* end while */ 
} /* end putim */ 


/* Takes the source filename, fs, strips it of its extension, copies it to 
* the destination filename, fd, and appends the new extension, ext, to the end. 


a 
void 
chgext (char *fd,char *fs,char *ext) 
{ 
char *period; /* indicates position of the period */ 
int len; /* gives length of filename less extension */ 
period = strrchr(fs,'.'); 


len = strlen(fs) ~ strlen(period); 
strncpy(fd,fs,len); 
fd{len) = ’\0’; 
strcat(fd,ext); 
} 


/* End of file SEMIO.C */ 
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2 MATLAB Seript file SEM.M used for plotting 


k Hıstogram of Particle Volume 
& Input variables are: 
Z total : Total Number of Particles 
Z totalv : Total Volume of Partıcles 
Z Nfiles : Number of data files merged 
Z limit : array containing the limits of the bins 
2 volume : array containing the volume in each bin 
clear 
clg 
cre 
BINS = 38; Z No of bins 
fprint£(*ininPLOT HISTOGRAM FUNCTION\n\n’ ); 
file = input(*MAT-fLilename to plot: ‘,"s’); 
eval(['load ' file)); 
zs = 0; 
ans = input('Do you want to input Malvern data [Y/N]? ’,’s’),; 
if (ans == 'Y” | ans == 'y') 
format compact 
ms = zeros(1i:BINS); Z allocate BINS 
for i = 1231 Z Malvern has 31 bins only 
ms(i) = input([num2str(1) *. Range ” num2str(limit(1)) * -— ’ 
nun2str(limit(i+1)) * Percent Vol = J); 
end 
for 1 = BINS:-=1:1 Z form histogram for malvern data 
zi = ms(i),; 
zs = (zs 2i 2i 0; 
end 
end 
xs = limit(BINS+1), 
ys ero: 
formi = BINS:-1:1 % form histogram for SEMEX data 


xh = limit(i): 


Ree ama ati): 
yi = volume(1)*100/totalv; 
xs = {xs xl xh xh]; 
SIS yi yi C); 

end 


ymax = max({ys zs]); 
ypos = ymax/20; 
axis({-1, 2.3, 0, max(ymax)typos) ); 
if (ans == 'Y’ | ans == 'y?) 
semi logx(xs,ys,'—-' |, xS,285, =="): 
bext (30, ymax-ypos*7, °—— Malvern’); 
text(30,ymax-ypos*8,’ = SEMEX’ ); 
else 
semilogx(xs,ys); 
end 
title(’HISTOGRAM OF PARTICLE VOLUME’ ); 
xlabel('Particle Size in Microns (Log scale)’); 
ylabel(’Percent of Total Volume’); 
text (30, ymax-ypos*2, [’Merged from ’ num2str(Nfiles) ’ images’]) 
text(30,ymax-ypos*3, ['Filename: * file '.mat']) 
text(30,ymax-ypos*4, ['Total Vol: ” num2str(totalv) * um3']) 
text(30,ymax-ypos*5, [’Partacle Count: ’ num2str(total))) 
axis; 


# End of tile SEM.M 
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